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Expressions very similar to the “Breit-Wigner” formulas are derived for the scattering and 


absorption cross section of nuclear particles, e.g. neutrons, protons, etc. The resonance energies, 
particle widths and other characteristic features can be expressed in terms of the logarithmic 
derivative of the particle wave function at the nuclear surface. Simple assumptions about the 
behavior of the incident particle inside the nucleus are used to find approximate expressions 
for (1) the particle widths in terms of the level spacing; (2) the neutron scattering cross section 
near and between resonances; (3) the total cross section and the absorption cross section 
averaged over neighboring levels for neutron energies up to 1 Mev. These results are then com- 


pared with the experimental data and are found to be in qualitative agreement. 





1. INTRODUCTION 


HE scattering and absorption of neutrons, 
protons, or other particles by nuclei is 
currently described by a picture that makes use 
of the existence of certain resonance states of the 
compound nucleus. The cross section of these 
processes, and its dependence, on the energy of 
the incident particles is given by the so-called 
“Breit-Wigner’’ formula’? which is built in close 
analogy to the corresponding expressions in 
optical dispersion. It has the following form for 
an incident particle of an angular momentum /h, 
if the spins of the particle and of the nucleus are 
disregarded :* 
1G, Breit and E. Wigner, Phys. Rev. 49, 519 (1936). 
*H. Bethe and G. Placzek, Phys. Rev. 51, 450 (1937). 
_* Consideration of the spins of the nucleus and the par- 
ticle introduces a factor (2J+1)/[(2i+1)(2s+1)(2/+1)] 
=£(J, i, s, 1) to the expressions (1) and (2), where ¢ is the 
spin of the bombarded nucleus, s of the particle and J of 
1e compound nucleus; J can take on the values J = |i+j|, 
rg -++|¢—j|, and j the values |/+s|, |J+s—1], 
-++|1—s|. It can be shown that 2, &(J, i, s, 1) =1, if the 


sum taken over all possible values of J. All expressions, 
in which averages are taken over many resonance levels, 


are therefore identical for the case s=0, i=0 and for the 
general case. 
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Ose is the cross section for elastic scattering, cabs 
is the cross section for absorption or any other 
process that changes the initial state of the 
bombarded nucleus. It includes radiative capture, 
fission, and inelastic scattering. Here k is the 
wave number of the incident particle, Z is its 
energy, E, are the resonance energy values for 
the compound nucleus, I“ is the total width 
of the resonance, I, is the neutron width, 
corresponding to the reemission of the neutron 
with its original energy, I’. is the absorption 
width and ¢, is a phase. I,‘ contains the radi- 
ation width and also the width for any other 
process that does not lead to the elastic reemis- 
sion of the incident particles, such as fission or 
inelastic scattering. The sum is taken over all 
resonance states of the compound nucleus. The 
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amplitude p.., contained in the formula for ox, 
describes the so-called “potential scattering”’ 
which is not caused by resonance effects. 

The derivation of these expressions and the 
theoretical determination of the constants ap- 
pearing in them have been given by several 
authors.” *-* The difficulty lies in the fact that 
the strong interaction between the impinging 
particles and the nucleus prevents a perturba- 
tion calculation in analogy with the dispersion 
theory of light. Especially the definition and the 
properties of the levels far off resonance have led 
to great complications and the theories were not 
able to give an unambiguous result. These levels 
contribute, mainly, to the background between 
resonances and are of special importance for the 
scattering cross section, since they are respon- 
sible, partly or fully, for the potential scattering. 

We propose in this paper a different way of 
approach to a derivation of expressions for the 
absorption and scattering cross sections. The 
result is, of course, very similar to (1) and (2), 
but it does not-contain a summation over far 
away levels. The cross sections are determined 
by a function f(£), which is equal to the ratio 
‘of the derivative to the value of the eigenfunction 
of the incoming particle at the surface of the 
nucleus. It is possible to express some charac- 
teristic quantities, as the neutron width I,” and 
the resonance energies E,, in terms of some 
general properties of this function f. In order to 
obtain numerical values for the cross sections, 
certain simple assumptions will be made as to 
the general behavior of the function f(Z). These 
assumptions cannot be exact and the results 
derived from it must be interpreted as qualitative 
indications only. : 

The fundamental idea of this derivation can 
be expressed in qualitative terms in the following 
way. Let us consider an incident beam of neu- 
trons of low energy. It is represented outside the 
nucleus by a wave of a wave number k. Near the 
nuclear surface this wave is joined smoothly with 
equal value and derivative to the wave function 
inside the nucleus, which is the solution of a 
complicated many-body problem. We expect 
that the neutron will have, in the average, a very 

‘*P. L. Kapur and R. Peierls, Proc. Roy. Soc. A166, 
277 (1938). 


SA, { Siegert, Phys. Rev. 56, 750 (1939). 
*G. Breit, Phys. Rev. 69, 472 (1946). 
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Fic. 1. Schematic representation of neutron wave func. 
tions at the nuclear surface. The wave functions are 
indicated as functions of the distance 7 from the center 
of the nucleus. r =a is the nuclear radius. Case a corresponds 
to a neutron energy between resonances, Case 6 is near 
resonance, Case ¢ is in resonance. 











high kinetic energy inside the nucleus. The de. 
pendence of the wave function inside on the 
coordinate of the incoming neutron will therefore 
be qualitatively represented by a wave with a 
high wave number K>k. We thus have to join, 
at the nuclear surface, a wave of low wave 
number k with a wave of high wave number K. 
In general, this can only be done if the-amplitude 


‘A inside is very much smaller than the ampli- 


tude outside: which we conveniently normalize 
to unity..A is of the order k/K. In the exceptional 
cases, however, when the derivative of the inside 
wave function is near zero at the surface, the 
two waves can be joined with about equal ampli- 
tude: A~w1. (See Fig. 1.) There are certain 
narrow intervals in the energy for which the 
wave function inside will have this exceptional 
property. These are the excitation energies of the 
compound nucleus for which the neutron can get 
into the nucleus, since, for these values, the 
amplitude is large inside. These energies are the 
observed resonances, and the widths of the inter- 
vals are the resonance widths. If we assume that 
the phase of the inside wave function at the 
nuclear surface changes smoothly with the 
energy, we will expect that the energy intervals, 
in which the derivative is small enough to give 


rise to large A, are proportional to the distance. 


D between resonances: [=aD where a is a 
small number of the order k/K. Thus the width 
of the resonance is expected to be proportional 
to D. (This consideration does not consider any 
absorption inside the nucleus so that this con- 
clusion applies to the neutron width only.) This 
simple picture also gives some indication of the 
scattering cross section. The outside wave func- 
tion assumes a very small value for all energies 
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at the nuclear surface except the ones near the 
resonances. Therefore, the situation is similar to 
the scattering of an impenetrable sphere of the 
same size as the nucleus, since in the latter case 
the wave function vanishes at the surface. We 
thus expect, in addition to a resonance scattering, 
a potential scattering between resonances whose 
cross section is similar to the one of an impene- 


trable sphere. 
2. THE CASE OF NEUTRONS WITH /=0 
(a) Derivation of the General Expressions 


The wave function y of the incident particle 
outside the nucleus is a linear combination of an 
incoming and an outgoing wave: 


g= ry =o ikr 4 ne**", (3) 


where k is the wave number of the incoming 
neutron. If there is no absorption, the absolute 
square of 7 is unity; in case of absorption, 
|n|?<1. The absorption cross section is given by 


Cavs = (4 /k*)(1—|n|*), (4) 


which is easily understood, if one remembers that 
x/k? is the maximum possible absorption cross 
section which is reached when there is no out- 
going wave e+** at all. 

The scattering cross section is given by 


Oac= (x /R*)|1+n]?*. (5) 


This can be seen by observing that, in a plane 
wave the /=0 part of the wave function has the 
form ry=e—*"—e**': this differs from (3) by 
(1+n)e+*r which, therefore is the scattered 
wave. 

We now express the two cross sections in terms 
of the logarithmic derivative of ¢ at the nuclear 
radius a. We define a function f(Z) by 





f(E) =al(de/dr)/ 9] ,=c. (6) 

The value of 7 is connected with f by 
- =") : 
n=e lay ; (7) 


where x =ka. The nuclear radius a is the smallest 
distance from the center of the nucleus at which 
the incoming particle is no longer under the 
influence of the nuclear forces. 
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We obtain an expression for the cross sections 
by inserting (7) into (4) and (5) and by putting 
f=fo—ith, where fo and h are real functions; we 
note that h/ is necessarily non-negative in order 
to insure | »|*<1. 























4nxh| 1 |? 4x xh 
Tabs = —] ee, (8) 
k? |x+if k? (x+h)?+fo? 
4x|x cosx — f sinx|? 
Bl. ef 
4a x 3 
= —|———-+e* sinx| . (9) 
k? |i(x+h) —fo 








It is seen that resonances occur always when 
fo=0. We therefore call the values EZ for which 
fo(Z) =0, the resonance energies EZ, of the nu- 
cleus. In the neighborhood of the resonances we 


ut 
fo(E) = (E—E,) (dfo/dE)e=8, (10) 


and introduce the following abbreviations (the 
negative sign in (11) and (12) will be explained 
later and is caused by the fact that dfo/dE can 
be shown to be negative) : 





r,, = —2x/(dfo/dE)&=£,, (11) 
r,”) = —2h/(df/dE) =e, (12) 
Then 
T r,r,” 
Tabs 


~ ke (E—E,)?+(P, +1PQ)?2/4 
EnearE,, (13) 
An ir,” oe 2 
xc =— +e sinx] , 
k* |E—E,+(t/2)(T,. +1.) 
Enear£E,. (14) 





These expressions show the familiar forms (1) 
and (2) although they contain only contributions 
from one level. They are valid only as long as 
(10) is a good approximation for fo(Z). Formulae 
(9) and (14) contain the characteristic ‘‘poten- 
tial scattering” term which, if it were present 
without the other term, would give rise to a 
scattering cross section : 


Osc = (44/k?) sin*ka. 


This is equivalent to the scattering of an impene- 














trable sphere of radius a. It must be said, how- 
ever, that the split of the scattering amplitude 
into a “potential scattering” term: e** sinx and 
a “resonance term”’: 


4°, /(E—E,+i(T,+T,) /2) 


is so far purely formal since no use was made of 
the fact that a is the nuclear radius. Actually a 
could have been any radius r for which the wave 
function of the incoming particle is given by (3) ; 
this means a could assume any value larger than 
or equal to the nuclear radius. It is evident that 
any change in value of a entails a corresponding 
change in f(Z) which makes the values (8) and 
(9) of the cross sections invariant to the choice of 
a. It will be shown in the next section, however, 
that f(Z) is expected to have certain simple 
properties which depend only on the nuclear 
structure inside the nucleus, in case a is chosen 
to be the nuclear radius. Specifically, relation 
(10) is then supposed to be valid over an energy 
region much larger than the line width [, +1”. 
In this case expression (14) is valid also for 
energies at which the resonance term becomes 
smaller than the potential scattering term, so 
that the potential scattering assumes real sig- 
nificance. 


(b) Definition and Determination of the Func- 
tion f(£) 


It is possible to make qualitative statements 
as to the behavior of the function f(Z), if the 
wave-length X of the wave function outside of 
the nucleus is much longer than the average 
distance d between particles inside the nucleus: 


k=1/k>d. (15) 


The actual value of X is then of no importance 
for the shape of the wave functions inside since X 
is very much longer than the wave-lengths inside 
which are of the order d. The value of f(£) will 
depend on the energy £ of the incoming par- 
ticle, only by being an explicit function f(W) 
of the excitation energy W of the compound 
nucleus: W=E+B, where B is the binding 
energy of the neutron to the nucleus. Expression 
(11) shows therefore the characteristic propor- 
tionality of I,“ with k, since fo or its derivative 
should not depend explicitly on k. 
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The actual determination of f(W) itself can 
be done only by finding the solution of the 
Schrédinger equation for the compound Nucleus 
corresponding to the energy W, and by cal. 
culating the logarithmic derivative of this soly. 
tion with respect to r at the nuclear surface. 
Because of condition (15), the value of this 
derivative does not change appreciably outside 
the nucleus in distances of the order d. The func. 
tion f(W) is therefore well defined, in spite of the 
fact that the actual position of the nuclear 
surface is not sharply determined. The Schrg. 
dinger equation for the compound nucleus can 
be written in the form: 


Ay (rn: + ra) = WV (ri -- +1), (16) 


where 7;:--ra are the coordinates of all con. 
stituents. Since W is the energy of the compound 
nucleus after the neutron has entered from out- 
side, it belongs, strictly speaking, to the con- 
tinuous spectrum of the Hamiltonian Hy. We 
assume, for the sake of simplicity, that no other 
particle but the incoming neutron can be 
emitted at the energy W. The wave function ¥ 
is then completely defined by the boundary 
conditions 


W(ri---ra4)=0 for re= 0, i=2---A, 


if r; are the coordinates of the incoming neutron. 
We then obtain a definite value for 


| -o(—* / (rw) (17 


since, for r,=a, YW has assumed the form of a 
product ¥(r1)x(re---ra), where x is the eigen- 
function of the initial (or residual) nucleus. 
This function fo(W) is not yet the function 
f(W) in which we are interested, since the 
damping effect of the radiative absorption has 
not been considered yet. It is possible to express 
the effect of the radiative absorption by an 
imaginary addition to the Hamiltonian Hy: © 


AV=(Ko+H)¥, H'=—T,/2. 


The term H' would add to the eigenvalues of Ho 
an imaginary term —iI',/2, which describes the 
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characteristic exponential time dependence : 


|v| 2m | e— #4 o/2)¢ | 2—¢ Tat, 


The solution of H which is needed to determine 


f(é) is given by 
Hv=Wv 


which also can be written in the form 


Hw =(W+4T./2)¥. 


: Hence, if fo(W) was determined according to 


(17) from Eq. (16), the actual f(W) will be 

given by 

f(W) =f W+iT 2/2) 
~fo(W)+1(T/2)(dfo/dW), (17a) 

if T,/2 is very small compared to W. The 

imaginary part —h of f(W) is thus 


h= —(Te/2)(dfo/dW), 


which agrees with our definition (12) of Ta. 
The resonance value W,=B+E£, and the 
neutron widths I,‘” are defined by: 


f(W,)=0, Th = —2x/(dfo/dW) wow, (18) 


Since the imaginary part of f(W) is necessarily 
non-positive in order to fulfill in (7) the condition 
|n|2<1, we conclude from (17a) 


df,>/dW <0. 


In view of the fact that the derivative of fy is 
never positive and that fo has zero points at 
W=W,, it is convenient to write it in the form 


fo(W) = —Ka tan 2(W), (19) 


where K is a constant of the dimension of a wave 
number and 2(W) is a monotonically increasing 
function, which is equal to multiples of x for the 
resonance values W,: 


2( W,) =nr, 


The value of K can be estimated in the fol- 
lowing way: the wave function W inside the 
nucleus depends on the coordinates of all con- 
stituents. We make the qualitative assumption 
that ¥, in its dependence on r, has the character 
of a periodic function with a high wave number: 


r¥~cos(Kr+5) (21) 
from which one obtains (see (6)) 


f= —Ka tan(Ka+4). 


nm an integer. (20) 





If this is compared to (19), it seems plausible to 
assume that K is a wave number corresponding 
to the kinetic energy of a particle in the nucleus. 
We would like to emphasize that this pro- 
cedure is not equivalent to the use of a one- 
particle model as the form of (21) may perhaps 
suggest. Expression (19) uses only the fact that 
W has a periodic dependence on r with a period 
given by K. The function z is left, so far, com- 
pletely undetermined. It would have the form 
z= (2mW/h*)‘a in the one-particle model. 

The order of magnitude of K can be obtained 
in the following way: the maximum momentum 
P of N equal free Fermi-Dirac particles in their 
lowest states, if kept within a volume V, is given 
by 

P=2nh(3N/4rV)}. 


If A is the number of constituents in a nucleus, 
we put N=A/4, because of the 4 types of par- 
ticles in the nucleus and get: 


P/h=(2m/ro)(9/64x")', 


where fo is the elementary radius which deter- 
mines nuclear radii by a=79A'. If ro>=1.5 X10", 
we obtain P/hk~1.0X10" cm which corre- 
sponds to an energy of 20 Mev. The average 
kinetic energy of a neutron entering the nucleus 
will be of the order of this energy plus the binding 
energy B. The wave number K should therefore 
correspond to an energy of about 28 Mev which 
would be: 


K =1.2X10"* cm. (22) 


The function f(W) would be qualitatively 
determined if some assumption could be made 
as to the behavior of f(W) between the values 
nx which it reaches at W,. It seems natural to 
assume that z(W) increases smoothly from nx to 
(n+1)z if the energy W goes from one resonance 
to the next. This assumption can be expressed by 


(dz/dW)w,=2/D* (23) 


where D* is an energy of the order of the average 
level distance D. D itself is, of course, a function 
of W. The function fp can then be written as 
function of E by means of (19) in the form: 


fo(EZ) = —Ka tan[(x/D*)(E—E,)]. (24) 


It is evident that this expression has only very 
approximate validity. 
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Fic, 2. Function 2(Z) (Eq. (19)) plotted vs. the energy 
E, E; to Es are resonance energies. The solid curve gives 
a smooth dependence as assumed in this paper. The 
broken line gives an example of an anomalous dependence 
which would lead to low absorption cross sections. 
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The relation between D* and the average 
level distance D depends on the detailed proper- 
ties of the function z(W) which, so far, is only 
defined by (20). If z(W) is laid as smoothly as 
possible through its values given by (20), we 
obtain D* =D in case of equidistant levels. It is 
seen in Fig. 2 that fluctuations in the level 
distance tend to decrease D*/D, since the slope 
of z at resonances will be closer to the slope of 
the straight connection with the neighbor reso- 
nance which lies closest. It seems to be a good 
definition of a ‘“‘smooth curve” to assume that 
D*, as defined in (23), is the smaller of the two 
distances between the resonance W, and the 
two adjacent resonance values. If this definition 
is adopted, one can show that a random level 
distribution gives 


D*/D=}. (25) 


The assumption of a “‘smooth”’ function 2(W) is 
somewhat arbitrary, and may not be justified 
in every case. One may argue, however, that the 
many-body problem of the nucleus is in some 
way equivalent to a one-body problem with a 
coordinate r which has a much larger interval at 
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its disposal than 0<r <a, because of internal re. 
flections and exchanges. The eigenfunction inside 
the nucleus has then the form C cos[K(r+ s)] 
where s is a length much larger than a, and K the 
average wave number for the excitation energy 
W. One then expects a smooth dependence on W 


‘of the function z=K(W)-(a+s). One can, of 


course, explain any value for D*/D by a suitable 
function z(W). The dotted curve in Fig. 2 gives 


an example, which would lead to an abnormally 


small value. 


3. THE ABSORPTION AND SCATTERING CROss 
SECTION OF s-NEUTRONS 


The assumptions made in the preceding 
section give rise to a simple expression for the 
neutron width I',, which becomes according to 
(18), (19), and (23): 


2k D* 
Yr,” —— oe (26) 


Tv 


The neutron width for s-neutrons can be deter- 
mined experimentally from an analysis of the 
neutron absorption resonances. The material 
available indicates that the values of I, are all of 
the order of 10-*4/E, ev where E£, is the resonance 
energy in electron volts. This value would give 
rise to D*~10 ev which seems to be a reasonable 
order of magnitude for the resonance level 
distance. 

Let us consider the average of gaps Over an 
energy interval AE which contains many levels: 

1 7 


r,T dE 
(abs) Av — — 
D k? (E—E,)?+(Tn+Ta)?/4 
2xT, Ta 4n D* IT, 








(28) 


Here the integral is taken over the neighborhood 
of a typical resonance E£, within the interval. It 
can be shown from the more general expression 
(8) that only the immediate neighborhood of the 
resonances contributes to the integral. I’, and I’ 
are the average values of I,” and I," over the 
resonances included within AE. Expression (26) 
was used to eliminate I’, in the final expression. 

For very low energies the following inequality 


7 N. Bohr and J. A. Wheeler, Phys. Rev. 56, 426 (1939). 

















ide 
) 1 
the 


ble 
ves 


lly 


ing 
he 


26) 


8) 


on 
he 
lr. 
he 
6) 
n. 











holds: 's>T'n. Since, for heavy nuclei (A >100) 
r,~0.1 ev and D* is of the order of 10 ev, we 
fnd T.=I', according to (26) for E~7000 ev. 
We thus expect the following expression to hold : 


4e D* 250 
(Gsm w= D ~ (Eye)? 





10-**cm?, E<7000 ev 
A>100 (29) 


if Ew is the incident energy in ev and (25) is 
used. No measurement of (¢abs)w in this energy 
region is available at present for heavy nuclei 
(A >100). It is not possible to derive a similar 
expression for lighter nuclei because of the 
larger D which would restrict the limit of validity 
to even lower energies. The larger D will then 
make it impossible to take an average over 
many levels within the validity of the expres- 
sion. It should be noticed that the cross section 


(29) is an average over an energy region con- 


taining many levels and it should therefore not 
be applied to the thermal absorption cross 
section. . 

We now turn to the scattering cross section 
as a function of energy. The value of ¢,. for any 
energy £ can be calculated from (9) and (24). 
The evaluation gives rise to a line shape of the 
characteristic form given in Fig. 3. This figure is 
calculated for the indium resonance, and a value 

*=25 ev was used which would give the ob- 
served neutron width® of 3X 10-* ev in expression 
(26). The depression at the low energy side and 
the high shoulder at the high energy side can 
also be understood as destructive and con- 
structive interference of the ‘‘potential scatter- 
ing’ amplitude with the resonance term as 
shown in (9). This phenomenon has been also 
described by Bethe.’ The extensions of these 
features in the energy scale are much larger than 
the width ([,+T,,) of the line. It is of the order 
2D/(Ka), which is about } of the level distance 
in case of heavy nuclei. 

No measurements of o. as a function of energy 
have been made so far, but the value of o4- for 
thermal energies or for energies near the thermal 
value are known for several elements. We expect 
these values to be near 4ra?. However, they may 

* J. Hornbostel, H. H. Goldsmith, J. H. Manley, Phys. 
Rev. 58, 18 (1940); J. Rainwater and W. Havens, unpub- 
lished ; we are grateful to these authors for communicating 


their results before publication. 
*H. A. Bethe, Rev. Mod. Phys. 9, 95 (1937). 
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be above or below this value if there is a level 
nearby below or above the energy at which the 
scattering was measured. The distance of this 
level must be at most of the order of about a 
fifth of the average level distance, in order to 
influence the scattering cross section. The nuclear 
radius a which determines the potential scat- 
tering is roughly defined as the distance at which 
the wave number of the incoming neutron 
assumes values of the order K. It is thus defined 
with an inherent inaccuracy of +1/K. We 
therefore expect the scattering cross section away 
from resonance to be near to 4rR’, where R is 
the nuclear radius determined from other sources. 
The deviation from 47R? should be of the order 
8xR/K which is about 37 percent of the total for 
elements near Fe and 25 percent of the total at 
the upper end of the periodic table. 

Table I shows a number of scattering cross 
sections measured, together with the values of 
4rR? based upon the assumption R= 1.5 X10-*A! 
cm. Of course, there will be fluctuations away 
from this average value. This table contains 
scattering cross sections mostly of elements 
which show small absorption at low energies. 





10,900 
- Tmax + 2420 x10 em? 


p—— P*008ev 











a 





xc 
(10 24 ome) 
6 
ad 


° SO 100 150 200 250 300 
Etev) 





























or 





Fic. 3. The scattering cross section for neutrons as a 
function of neutron energy in the neighborhood of a 
resonance. The numerical values used for T and EZ, are 
those for the indium resonance. Because of a typographical 
error the value of os. for thermal energy is given wrongly 


in this figure. It should read: ¢(thermal) = 1.4710 cm’. ~ 
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This selection is attributed to experimental 
reasons. It follows therefore that in the case of 
’ heavier elements (A >100) there is no resonance 
nearby to interfere with the potential scattering. 
The cross section should therefore be 4xa?, which 
should not differ from 42rR? by more than about 
30 percent. The values of Table I bear out this 
expectation with the possible exception of Sb 
and Hg. In these cases, the actual value of R 
may be considerably different from the assumed 
value R= 1.5 X10-"A", In the case of the lighter 
elements (A <100), however, the neutron width 
is probably much larger than the absorption 
width because of the larger D, so that resonance 
scattering may occur in spite of weak absorption. 
We expect then, in some cases, to find an inter- 
ference effect of resonance and potential scatter- 
ing which may account for the abnormally high 
cross sections in Fe, Ni, Cu, Ge, Se, and Sr, 
and for the low value in Mn. 

Finally, measurements” of the change in phase 
of the scattered wave with respect to the incident 
wave has been measured for several heavy nuclei 
for energies in or near to the thermal region. In 
most cases this phase change was 180°. This 
result is in agreement with the assumption that 
the potential scattering p,. is the scattering by 
an impenetrable sphere for, in this case, the 
outgoing wave must be exactly out of phase with 
the incident wave. This phase relation will hold 
as long as the potential scattering term is the 
most important term in gs. In the resonance 
region, this phase relation holds for neutron 
energies greater than the resonance energy. It 
does not hold for energies smaller than the 
resonance energy. This is, of course, just a very 
small portion of the neutron spectrum so that for 
most energies the phase change will be near 180°. 


4. GENERALIZATION TO NEUTRONS />0 AND TO 
OTHER KINDS OF INCIDENT PARTICLES 
(E.G., CHARGED PARTICLES) 


So far, we have considered only neutrons in 
the s state, i.e., /=0. When neutrons with / >0, 
or charged particles, protons, deuterons, a-par- 
ticles, etc., are considered, the radial equation 
for the particle for r >a now will contain potential 
energy terms. In the case of charged particles, 
1=0, the potential energy term is Zze?/r where Z 


1 E, Fermi and L. Marshall, Phys. Rev. 70, 103 (1946). 
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TaBLeE I. Elastic scattering cross sections o,, for rea al | 
and epithermal neutrons. r 








Se 
— 


4nR? Energy (ev) Reference 


Element 


Ti 
Mn 
Cr 
Fe 


Co 
Ni 
Cu 


a 
2 





i) — os 
Wh & in i 
= & = 


335 
AAANAAAA EK EKKAF EKKO 


NoOKAW Ww 
wn 


— — 
oF RK AQQ--- 
in r 


— 
CRACK ASNOMPUMNDA RAR 


Rua sS YS ft 
"mm 


SCLC CLC LCCONINAAAAAAAAMMGAH PP POW 
SOMDHPWNOUNDOCAUPRWHSOORUWO A UPR SHBOD 


— 
ane, aber 
oe 








afin” Goldhaber and G. H. Briggs, Proc. Roy. Soc. A162, 127 

B J. Rainwater and W. Havens, — We thank Rainwater 
and — for communicating these results 

CN. Bradbury, F. Bloch, H. Tatel, and P. A. Ross, Phys. Rev. 
52, 1023 wgsi). 
net LL, Hanstein, Phys. Rev. 59, 489 (1941). 

The fourth column gives the energy at which os, was determined. 
In all references but reference A, the value measured is always the 
total cross section ¢s¢-+¢«. In these cases, oa was either known to be 
much smaller than og, or it was possible to subtract unambigously 
term proportional to 1/y which was considered to be ea. 
neutron energies absorbable in Cd. The velocity selector measurements 
were taken mostly in the neighborhood of 1 ev (~1) except where 
other ener; are quoted. In the latter cases, the nearest level was 
already influencing the value at 1 ev, so that it was advisable to choose 
a lower value. 

The scattering cross section of Cd was determined as an additive 
term to the best fitting Breit-Wigner expression covering the resonance. 


is the nuclear charge and z is the particle charge. 
When /+0, then for both cases, neutron or 
charged particles, one must introduce an ad- 
ditional effective potential energy equal to 
h?l(l+-1)/2mr*. In any of these cases, the par- 
ticle may be considered to be moving in the re- 
gion r >a, in a potential V(r). It is then possible 
to generalize the discussion of Sections 2 and 3 
as follows. 

Let u(r) equal the product of r and the wave 
function of the incident particle moving in the 
potential V(r), which asymptotically for larger 
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behaves as e~**": 


u(r)—e rl and similarly v(r)—ret@r—t#/2) | 
10 1 2 


The wave function ¢ for the particle can then be 
written as in (3) as: 
g=r=utm, || <1. (31) 
The absorption cross section is given by 
Oabs = (2/-+1)(4/k*)[1 — | |?) (32) 
and the scattering cross section by 
Oe = (21-+1)(4/k*) |1+7|?. (33) 
A function f(Z) may be defined in analogy to (7) 
as: 


, U 


— 
f(B) =a(¢'/¢)=2—— (34) 





where w’ and v’ denote derivatives with respect 
to kr at r=a. Solving (34) for 7 we obtain 


xu’ — fu 
q= —— t (35) 
xu’ — fu 
It is convenient to express u and v and their 
derivatives in terms of their absolute value and 
phase: 


v= |v|e*, 
v’ = |v’ |e*’. 


u=|v|e—*, 


u'= lv’ | et (36) 


Note 6’ and |v’| are not the derivatives of 6 and 
|v|. If one uses the relation 
|wv’| sin(é’— 48) =1 for all r (37) 


one may transform (32) and (33) into 


xh/|o|* 





ie ‘(x/|v|?-+h)? +[fo—x|v’/o| cos(é’— 8) ]?’ 


x/|v|? ‘ 


(38) 





Tv 
Osc ere 


Here f=fo—ih, |v|, |v’|, 6 and 8’ are evaluated 
at x=ka. The resonances occur when 
fo=x|v'/v| cos(d’—8). (40) 


In the neighborhood of the resonances, we expand 
fo(E) as in (10) and introduce the following 
abbreviations 


rtm -=(5), = (41) 


dfo 
r = — 21> =) (12) 
E=E, 


We obtain 
r, (Tr, 
, (42 
(E—E,)?+((0, +1.) /2}? 





7 
Tabs “St 1) 


4a 
alt ial 1) 


sr,” 2 
+Pe| » (43) 
(E-E,) +i/2(7, +1.) 





with re 
Po =e* sind. 


+sinée*| . (39) 


i(x/|v|?+h)—[fo—x|v'/v| cos(s’ — 8) ] 





These expressions are almost identical with the 
ones obtained for neutrons with /=0. The term 
in ¢sc, Which represents the potential scattering 
is changed, and the expression (41) for the 
neutron width contains the characteristic factor 
|v|-* which expresses the effect of the potential 
field outside the nucleus. This factor is the ratio 
|v( )/v(x)|* of the intensities of an outgoing 
wave at infinity and at the nuclear surface. It is 
a small number for repulsive fields. Its ap- 
pearance in the expression of the neutron width 
is expected, since the nuclear eigenfunction 
adjusts its form to the wave function of the 
particle at the nuclear surface, whereas the width 
is connected with the current of outgoing par- 
ticles through a sphere around the nucleus with 
a very large radius. The factor |v|~* has also 
been used in earlier papers by Kapur and 
Peierls,‘ Weisskopf and Ewing," and Bethe.® 

There is a limitation to the validity of these 
considerations, which corresponds to the con- 
dition (15) of the case without barrier. Here the 


u'V. F. Weisskopf and D. H. Ewing, Phys. Rev. 57, 
472 (1940). 
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condition is: 
1/x>d, x=k|v'/v| rma; 

where «x is the instantaneous wave number at 
the nuclear surface. This condition is necessary, 
in analogy to (15), to insure small variations of 
the wave function outside, in intervals of the 
order d. It puts a limit of about several Mev’s 
to the height of the barrier to which our expres- 
sions are applicable. 


5. APPLICATION TO NEUTRON CROSS SECTIONS 
FOR HIGHER ENERGIES 


When different / values are involved, the func- 
tion v will depend upon |. This dependence will 
be indicated by a subscript /. In the case of 
neutrons v; is given by: 


v,=1(9x/2)*Hi44 (x). 
The value” |v;|? is then 


|o2|?=$ax| Hy 4 (x) |’, 
|vo|?=1, 

|v |? = (1+-2*) /x?, 

| v2|2?=(9+3x2+24) /x4, 

| vs |? = (225+45x°+ 6x4+x°) /x°®, 


The function 6,(x) is: 
Ji+4(x) 

5:(x) = —tan— ; 

Ni+4(x) 

69 =x, 

6;=x—2/2+cot™ x 

bg =x—2+cot—!4(x?—3), 

63 =x—34/2+cot—[x(x?— 


The width I, for the emission of a neutron of 
angular momentum / is therefore given by 


2k 1 Di* 
K |v,|? T 


15) /(6x?—15)]. 


r= (46) 


if the same assumptions about fo(Z) are made as 
in (24). This is the generalization of (26). 
These expressions can be applied to the ab- 


1% The functions |:(x)/x|*, |v:’(x)/x|?, 8G), 8’(x) have 
been tabulated for /=0(1)20 and x =0(.1)10 n “Scattering 
and Radiation from Circular Cylinders pond Spheres,” 
Morse, Lowan, Feshbach, and Lax, a report issued by the 
N.D.R.C. Division (6). 


(44) 
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sorption cross section of neutrons. Let us con- 
sider the average cross section over a small 
energy integral AE containing many resonances, 
This is the cross section observed, if the neutron 
beam is not sufficiently monochromatic to dis- 
tinguish resonances. The average over (42) for 
a given / gives: 


2m(T nil ai)av 
Di((Tnddw+(T adm) 


Here D, is the level distance of resonance levels 
excited by particles with an angular momentum 
Ih; (Tni)w and (Ta:)4 are the average neutron 
width and absorption width, respectively, for 
neutrons of angular momentum /h in the energy 
region EZ. This can be transformed by means of 
(41) and (23) and (46) into: 
D* 1 


(edn = 21+-1)— 
( Ss, (1+(DPnz)av/{Tat)m) | 02]? 


(47) 
where |9;|? is given in (45). We now assume that 
the ratio (I'a:)~/D:* and the wave number K is 
independent of / and of the energy. One then can 
write: 





7 
(Cabs)? = ra +1) 





(T'nt)a (T'no) av x 
(Tat)av ~ (Peo)w | 02]? ~ xolo |?” 





where x9 is the value of x =ka for which 


a Ka 
“a >To 


(Tno)w=(T'ao)w, Xo 


The total absorption cross section can be obtained 
by summing (47) over all /’s: 


x «© D* 


(Cabs) av = aC +1) 


————.. (48 
|v: |?-+2/xo “ 


The assumption that (I'a:)w/D;* is inde- 
pendent of the energy will break down when 
inelastic scattering or fission sets in which, in _ 
this treatment, are included in (Tai)w. This 
increase in (I'a:)4 will make xo increase sharply 
at higher energies. We may assume that the 
inelastic scattering is unimportant for energies 
less than 0.4 Mev. 

The value of x9=4xKaI,/D will vary from 
element to element but it is expected to be much 














smaller for lighter elements than for heavier ones. 

This comes mainly from its dependence on the 
level distance D. We expect that D is substan- 
tially larger for light elements. The ratio T,/D 
for any width I’, corresponding to an individual 
transition is usually considered to be independent 
of D. However, Ta: is the sum of all individual 
transition widths to lower states of the compound 
nucleus. The number of possible end states is 
roughly inversely proportional to D, so that 
T'z:/D is expected to decrease considerably with 
increasing D. 

The variation of D seems to occur, however, 
only at values of A<100. The fact that the 
neutron widths in known nuclear resonance 
levels show approximately the same proportion- 
ality factor of ./Z in elements different as 
indium (A=115) rhodium (A=103) and Au 
(A =187) indicates, according to (26), that the 
level distance D at the neutron resonances does 
not change appreciably for these values of A. 
This is probably the effect of two opposing 
trends: (a) The level distance D decreases with 
increasing A for a given excitation energy of the 
compound nucleus. (b) The excitation energy at 
which D is measured is the binding energy of the 
neutron, which decreases with increasing A for 
A>100 and which is roughly constant for 
A<100. For A>100 the two effects a and b 
operate in opposite directions, giving rise to D 
and consequently to an x» which is comparatively 
constant. In the region A <100, only the trend 
given in (a) remains, which will give rise to an 
increase in D and thus to a decrease in Xo for 
decreasing A. This probably explains the marked 
decrease of the absorption cross section, ¢abs, for 
elements for which A<100 as observed by 
Hughes.” 

Expression (48) is, strictly speaking, valid 
only for nuclei with spin zero. If the spin is not 
zero, reemission of the neutron is possible with 
an / different from the incoming one, without 
exciting the nucleus. The partial widths cor- 
responding to these processes should be included 
in (I's1)~ since they represent an “‘inelastic”’ scat- 
tering inasmuch as the state of the nucleus has 
changed during the process. Thus (I'a:)4 would 
be made larger than the pure radiation width, by 


an amount which depends on the neutron energy. 


18 Hughes, Phys. Rev. 70, 106 (1946). 
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Closer investigation and application of selection 
rules reveals, however, that this effect is not very 
important. It decreases oa», by less than 10 per- 
cent if x51. 

Expression (48) for the absorption cross section 
of neutrons can be compared with experiments 
performed by Greisen, Linenberger, Miskel, and 
Segré,!4 on the elements Ag, In, I, and Au, with 
neutrons of energies between 5 and 400 kev. 
Figure 4 shows these results plotted as a function 
of x?=k*a’, which is proportional to the energy. 
An energy scale is given for In and Au. The solid 
curves are plots of expression (48) with different 
values of xo. It is hard to predict the actual value 
of xo, even in elements like In and Au, where one 
resonance level is analyzed, since the properties 
of this level are probably not typical. The best 
values® for the In resonance are: 


r',.~3 X10(E,)tev, T.~0.08 ev 
and for Au: 
r,<~6X10-*(E,)tev, T.~0.15 ev. 


This gives x»~0.05 for both elements. The actual 
value of xo is probably somewhat higher since 
these two resonances are unusually strong which 
points to an abnormally large I,. Moreover, 
there is an additional uncertainty due to the 
fluctuations in the nuclear radius from the 
assumed average value. Figure 4 shows that the 
theoretical curves give a fair representation of 
the experimental material. 

The rather sharp decrease of oap./oo for Ag 
and Au for energies greater than a few hundred 
kev is caused by the onset of inelastic scattering 
of the neutrons. The cross section for this com- 
peting process is not included in the measured 
absorption cross section. Bismuth and lead," 
however, seem to represent notable exceptions 
since they have shown no observable absorption 
in this energy region. These five nuclei (lead has 
four isotopes) fall completely out of the rules 
observed with so many other heavy nuclei. 
Within the framework of this theory, this would 
only be explained by an unusually large level 
distance D for these elements, which would point 
to an unusually low binding energy of the last 
neutron. 


i "MWe are are my to these authors for communicating 
their results before publication. 











FESHBACH, PEASLEE, AND WEISSKOPF 











AB 
2 
s 




















/j 






































é ' / 
= ° 
e— “gs 
i Say = Ya rN / 
e Ps! a be LA i) Fic. 4. The absorption crogs 
INQ ~~ Paid el k be a ng | = neutrons as a func. 
— ve +1 . ° tion of (ka)*, The experimental 
a Pie a a Lo 70! oints Acuna eo by Grei 
= inenberger, Miskel, and 
: Phe] — ” ° We are grateful to these authors 
z nN ” wl for communicating their results 
bo _ es a Lo prior to publication. 






































































The elastic cross section of neutrons can be 
calculated with the help of (43). We consider 
again the average (os.)w over many resonances, 
since this is the magnitude which can be 
measured most readily. We obtain 


1 Ert+D/2 Mr 


sc) = — ween 21 1 
(ore) (+1) 


E,—D/2 


x | (Tut) 
(E —E,) +3((Pnt)av +(Ta1)m) 


4a mT niw\ . , 
-=a+1] (1- > ) sin 5: 











us r,? 
+- =| 
— 2 DT ata) 
With the aid of approximations (46) and (47), 
we may now write 


4 
(00)m = LX Kose)” =oo9+— — — 


+1 1 
"yo: |?L 1+ (wol vs] 2/2) 


where X =Ka, 6; is defined by Eqs. (36), (45b) 
and apo is: 





—2 sin’ (49) 


4n 
oo= 7 > (21+ 1) sin?6). 


} 49 190 400 1000 —s kev in 
2 28 28 70 280 | 700 kev Au 
° i i omen 
0005 oe 068 O& c2 oo 8 2 : «. =) USS 
x (ka? 


+Pe| dE 






This is the scattering cross section of an im. 
penetrable sphere of radius a, oo is 4ra? for smal] 
energies (ka<1) and approaches 27a? asymp- 
totically for large energies (ka>1). 

The quantity usually measured is not (oy 
but rather the total scattering cross section 


(o av = (osc)mv + (Cabs) av- 









The theoretical formula for (o:),, is obtained by 
adding (49) and (48) 


4n D* x 21+1 


(0 s)m = oo+— —— Lt 


? DX |v, |? 








cos26;. (50) 





Note that (01) is independent of x». 

The total nuclear cross section is, according to 
(50) always greater than oo. For low energies, 
one need only use the /=0 term in (50). 


(o1)w =42ra*{1+D*/DXx]. 








The excess of (o1)4 over 4xa? for small x will be 
large since the 1/x term (1/v law of absorption) 
becomes very large for small x. For larger 
energies, several values of / contribute to the 
cross section. In the limit the positive and 
negative terms cancel so that (o1)m—a0. 
Measurements of the total cross section of 
several elements of neutrons in the energy region 
between 20 kev and 800 kev are made by 
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Leipunsky'® and Wattenberg. Figure 5 shows 
these experimental results in units of a’, 
together with curves for oo and (o;)» for in- 
dium, and D*/D=}4. Curves for other nuclei 
may be calculated quite easily from this one by 
making use of the fact that (¢1)w—oo, from (50), 
depends inversely on the nuclear radius. In 
comparing the theoretical and experimental 
values, one must keep in mind that a was 
assumed to be 1.5X10-%A? cm. The actual 
value of the nuclear radius will probably deviate 
from this average value. In Fig. 5, representative 
elements from among the middle and heavy 
elements have been included. It is seen that the 
qualitative features of (o.)4 are in substantial 
agreement with the theory for the elements Fe, 
Ni, W, Pb. Although the cadmium points are 
of the right order of magnitude, they do not 
show the characteristic decrease with increasing 
ka. This behavior is exhibited by other elements 
with A~100 such as Sn, I, Sb, all of which have 
about the same (1) as Cd. 


6. RELATION TO OTHER METHODS 


The derivation of the resonance formulas pre- 
sented here is in many respects similar to the 
method employed by Kapur and Peierls.‘ It may 
be useful for the understanding of the situation 
to point out the similarities. 

It is possible to write the expressions (8) and 
(9) in a form in which the familiar sum over 
resonance levels appears, if the following relation 
is used : 


1 & 
E => —————-q= r . 0 
-_ f(E) —ix S,—E ies 





where S, are the poles of the function g(Z), which 
are defined by 


f(S,) =ix. (51) 


The C, are the residues: C,=(1/(df/dE))z8<s,. 
Here x is considered as a constant parameter inde- 
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Fic. 5. The total cross section for neutrons as a function 
of ka. erorat is the total cross section curve for A = 115. 
The eo curve is the cross section for scattering by an “im- 
penetrable sphere.” The experimental points are typical 
examples taken from Wattenberg (reference 17). Similar 
results have also been obtained by Leipunski (reference 16). 


pendent of E. It follows from (17a) and (18) that 
S,=E,—il'&) /2—4T,,/2 (52) 


where £, is the solution of fo(Z)=0. If the 
imaginary parts of S, are small compared to the 
real ones, the relation C,=—TI, /2x holds 
because of (18), and we can write (8) in the form 


- (T,(T,”) 2 
Tabs — 


k? e_ E+. 4T.)/2 





(53) 


after using the relation 


h= —41,“ (df/dE)z-s,. 


Expression (9) can similarly be written in the 
form: 





Tr,” 2 





Ca =-— 
R2 





" E-E,+i(P, +1) /2 


46 A. I. Leipunsky, J. Phys. Academy of Science, U.S.S.R. 3, 231 (1940). 
1¢ J. Wattenberg, unpublished. We are grateful to the author for communicating his results before publication. 


+e sinx] . (54) 
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These expressions are the same as the ones 
derived by Kapur and Peierls. The eigenvalues 
S, are defined in their paper by the condition, 
that the wave y of the emitted particles satisfies 
the boundary condition 


dy /dr —itky =0, 


which is identical with our condition (51). 

The cross sections for a given energy E 
between resonances are represented in this form 
as the effect of contributions of a large number 
of levels E,, whereas in our case the same cross 
sections are determined by the function f at the 
energy £ in question. The two approaches are 
mathematically equivalent because of the iden- 
tity (50a). The method presented here seems to 
offer the advantage of an easier interpretation of 
the magnitudes involved. Qualitative conclu- 
sions as to the behavior of the logarithmic de- 
rivatives f(£) of the wave function at the surface 
of the nucleus can be drawn more readily than 
conclusions concerning a whole series of eigen- 
solutions of a complicated eigenvalue problem at 
energies far off the value E. Especially the levels 
far off resonance and their contribution to scat- 
tering offer some difficulties of interpretation. 
It was pointed out by Siegert,’ that they do not 
correspond to any physical state since the 
boundary condition contains the k of the incident 
particle and not the k which would correspond 
to the emission of the particle by the resonance 
level itself. Siegert’s. method avoids this difficulty 
by using a different boundary condition for the 
definition of the compound states. However, the 
physical meaning of compound states of high 
energy is questionable, when their width becomes 
larger than their distances, which is bound 
to happen at higher energies. 

Bethe” has derived expressions similar to ours 
for the average cross sections and particle widths 
by ascribing to the nucleus a strong absorption 


17H. A. Bethe, Phys. Rev. 57, 1125 (1940). 
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coefficient in form of an imaginary potentiaj 
energy. This absorption was supposed to describe 
the ‘‘amalgamation”’ of the incident particle with 
the nucleus. It therefore was assumed to be so 
strong that the wave function of the incoming 
particle decreases by a factor e~ at a distance } 
of the order of the reciprocal of our K. Bethe’s 
expression for the average neutron capture crogs 
section for / =0 (oabs)w = 22°kd [his formula (27a)] 
is identical with ours (28) if 


b=(1/K)(D*/D)(2/z). 


Bethe’s model does not exhibit any resonances 
because of the strong absorption assumed. This 
absorption is very much larger than the actual 
one due to radiative capture, and the corre. 
sponding breadth is so large that it would 
broaden any levels until they merge. The simi- 
larity of Bethe’s results with ours comes from 
the fact that his fictitious absorption causes the 
wave function of the incoming particle to 
assume a high logarithmic derivative after 
entering the nucleus of the order K~1/b. This 
fact alone determines the average cross sections, 
whereas the special form (24) of the logarithmic 
derivative is important only for the finer details 
relating to resonance phenomena. Bethe’s as- 
sumption of an absorption coefficient of the 
order b~! appears to be equivalent to our assump- 
tion (21) about the average wave numbers inside 
of the nucleus. Our way of introducing high wave 
numbers inside the nucleus seems to be more 
general since it also permits the description of 
resonance phenomena. 

The authors express their appreciation to Dr. 
H. H. Goldsmith for his invaluable help in 
analyzing the unpublished experimental material 
obtained from the recently declassified work of 
the Manhattan Project. They also wish to thank 
Dr. G. Placzek for many enlightening discus- 
sions.'® 


18 The research described in this 5 was supported 
in part by contract NSori, U. S. Navy Department. 
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The direct production of RaE and Po by deuteron bombardment of bismuth has been rein- 


vestigated. Present results indicate that as the deuteron energy is lowered from 8.8 to 7.3 Mev, 
the production ratio of RaE/Po increases 25 percent but as thé deuteron energy is dropped from 
7.3 Mev to 5.5 Mev, the ratio more than doubles. The trend is in agreement with the Oppen- 
heimer-Phillips reaction in which RaE is produced by a d, reaction and Po is produced by a 
d, n reaction. Assuming these reaction cross sections to be:respectively proportional to Bethe’s 
OP penetration factor and to the Gamow-Condon-Guerney penetration factor, best agreement 
of the variation with energy of the ratios of these cross sections is obtained with a nuclear 


radius of 2.0A?X10-* cm for bismuth. 








INTRODUCTION 


_— disintegration products of bismuth 
under bombardment by deuterons with 
energy in the region of 5 to 9 million electron 
volts has been studied experimentally and two 
different results were found.'? Since bismuth 
with its single isotope and high atomic number 
affords an excellent opportunity for the compari- 
son of two competing processes under deuteron 
bombardment—proton emission and neutron 
emission—the problem was reinvestigated experi- 
mentally at the University of Michigan Physics 
Department in 1940. Because of the war, the 
results have only just now been interpreted and 
are herewith reported. 


THE NUCLEAR REACTIONS 


The nuclear reactions have been previously 
reported and there is a d, p and a d, m reaction 
which may be described by means of three 
processes : 


I. 1d?-+-g3Bi?°—>g3RaE”+4', 
Il. 1d? +-g3Bi?—>g,* Po*! 4,3 RaE?!9+-1p!, 


III. d?+-83Bi2—4,* Po?!!—+g,P0?!°+-on!. 


The first two reactions have the same end 


anon” Lantham, and Lewis, Proc. Roy. Soc. 174, 126 
2 Cork, Halpern, and Tatel, Phys. Rev. 57, 348 (1940); 
57, 371 (1940). 
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products and so may not be distinguished by 
measurement of only the end products. Process 
I is an Oppenheimer-Phillips reaction* in which 
the deuteron splits up outside the nucleus with 
the neutron being absorbed in the final nucleus 
which is identical with the compound nucleus. 
Formally, reaction II appears the same as I but 
here an intermediate compound nucleus is formed 
and the proton emission must compete with the 
more probable neutron emission as in III. If III 
is more probable than II, only I and III will 
occur and these two are distinguishable by a 
fortunate set of circumstances. The RaE formed 
in I decays by B-ray emission with a half-life of 
5.0 days into polonium and the Po thus formed, 
as well as the Po produced in III, decays with a 
lifetime of 136 days by the emission of a-parti- 
cles. Thus, by the measuring of the a-particle 
activity as it changes, the amounts of each 
product may be determined. The number of Po 
nucleii N2 present at any time is given by 


N; = [ArW1°/(A1 — Az) JLexp( — el) —exp(—Az:) ] 
+N,2° exp(—Agf). (1) 


N° is the number of RaE nucleii formed by the 
bombardment in reaction I (and II) while N;° is 
the amount of Po nucleii present at t=0 (at the 
instant of bombardment) in reaction III. A; and 
A2 are, respectively, the decay constant of RaE 


*H. A. Bethe, Phys. Rev. 53, 39 (1938). 
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Fic. 1. Straight-line plot of alpha-particle emission as a 
function of time for the various targets with deuteron 
bombarding energy. Ordinates should be multiplied by 
1.10X 108, 1.60 108, 10.0 108, and 3.06 x 108 for the 0, 1, 
2, 2 mil curves, respectively. 


(1.6010-* second) and Po (5.89X10-® sec- 
ond-"). Ny° and N2° may be found from the 
experiment by integration or by differentiation 
of (1). For the latter method, the alpha-particle 
emission rate at any time is given by da/dt 
= —):N2 which can be combined with (1) to give 


= NC fi(t)/fo(t) ]+Ne°; 2 
= / i= »Lfa(t)/ fat) ]+N2°; (2) 


fi(t) = DAro/ Qu — a) J 
X [exp(— Ast) —exp(—Ait) ], 
fo(t) =r. exp(—r2/). 


Expression (2) is the equation of a straight line 
and so if the a-particle rate divided by fe is 
plotted as ordinate against f:/f2 as abscissa, the 
intercept N2° gives the yield in III and the slope 
N° is a measure of the yield in I (and II). 
Figures 1 and 2 show plots of the data in this 
form. The vertical lines through each point are 
a measure of the statistical fluctuation of the 
number of counts. This fluctuation is taken as 
N* where N is the total number of alpha-particles 
measured in a given interval. 
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THE OBJECT OF THE EXPERIMENT ANp \ 
GENERAL PROCEDURE 


For a given target of bismuth, it is possible to 
compare the cross sections for deuteron bom. 
bardment of the two types of reactions, emission 
of a proton and emission of a neutron. Generally, 
when a compound nucleus is formed at excita. 
tions below the barrier height, it is considerably 
more probable that a neutron is emitted than a 
proton because, for one thing, the neutron does 
not have to penetrate the potential barrier. 
However, if the mechanism proposed by Oppen- 
heimer and Phillips prevails, it is possible for the 
d, p reaction to be more probable. A test of 
whether or not the OP process occurs is to 
compare its cross section with an ordinary reac. 
tion as the deuteron bombarding energy is varied, 
the OP process diminishing relatively as the top 
of the deuteron barrier is approached. In this 
experiment, the energy of the deuterons was 
varied by interposing aluminum foils between 
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Fic. 2. Straight-line plot of alpha-particle emission as a 
function of time for the various targets with lower bom- 
barding -z. The ordinates should be multiplied by 
2.00 108, 2.00108, 10.0108, 10.0108, 263 10* for 
the 3, 4, 5, and 6 mil curves, respectively. In comparing 
Fig. 2 to Fig. 1, it is seen that there is greater slope diver- 
gence in Fig. 2 although the intercept spread is over the 
same range. This increase in slope indicates relatively more 
Ni(RaE) with lower energy. 
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RADIUM E AND POLONIUM FROM BISMUTH 


TABLE I. Experimental results. 
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° N = o 108 o 10%8 
| Mer. 10-7 xo? = “coulombetimebr. | -Nw/Nge R ox as 
0 8.7 55.2 11.3 0.012 0.58 4.9 5.0 34 6.8 4.3 4.9 
1 8.2 36.9 7.3 0.015. 0.5 5.1 5.2 18 3.4 5.0 5.0 
2 7.8 20.9 3.7 0.0165 1.0 5.6 5.7 9.2 1.6 4.9 6.0 
2 7.8 29.9 5.3 0.015» 0.5 5.6 5.7 15 2.6 
2* 78 7.8 1.3; ° 0.027; 1.5 5.8 6.0 6.3 *\ 6.3 
3 7.3 17.8 2.78 0.016. 0.83 6.2 6.3 8 1.30 8 
4 6.8 27.2 3.2 0.040,(?) 1.8 8.1 9.5 5 0.5 17 
5 6.3 7.5 0.72 0.0367 2.0 10.2 10.5 1.5 0.14 6.5 
6 6.3 11.3 0.83 0.0565 2.0 13 14 1.4 0.10 
6 5.8 0.52 0.033 0.072 4.0 13.6 16 0.5, 0.03, = 











* Target of approximately 3.2 mg/cm? of Bi on Al. 


the beam and the bismuth target. Hence, for a 
given energy, a bismuth target was bombarded 
and the rate of emission of alpha-particles was 
measured from the time the bombardment 
stopped and plotted as already explained. Then 
N,° and N;2° which are the end products of the 
d, p and d, m reaction, respectively, may be found 
and so their ratio for each energy. Since for each 
target the deuteron energy and distribution is 
identical, an excellent value of the ratio of these 
two quantities NV = N,°/N,° may be found. This 
ratio must be corrected for the time of bombard- 
ment especially for the weak samples which have 
both a long bombardment time and a large ratio. 
The reason is that the radium E is always 
decaying into polonium and so at the end of a 
bombardment there will, in general, be less 
radium E and more polonium thereby giving a 
smaller ratio. The true ratio of cross sections is 
to a first and sufficient approximation : 


N 
R= 
1—4\,7(N+1) 


where \; is the decay constant of radium E and 
T is the bombardment time. This correction 
amounts to 1 percent for an N of 5.6 and a 
bombardment of 0.5 hour; and to 9 percent for 
a bombardment of 2 hours with an N of 13. 
These values with N and R and T are shown in 
Table I. 

There are two types of targets which have been 
used in this experiment and both are “thin.” 
HLL! used evaporated layers of bismuth on thin 
aluminum foils with a total thickness of about 
1 air-cm. The second procedure*—used here 
—employs a physically thick piece of bismuth 


(3) 


















for the target but makes use of the short 
range 3.8 air-cm polonium alpha-particles com- 
pared to the long range deuterons to secure an 
essentially thin target. Since the alpha-particles 
had to travel one centimeter before entering the 
ionization chamber grid to be registered, the 
only ones counted must originate within 2.8 
air-cm of the surface of the bismuth target. At 
the deuteron energy limits of 6 and 9 Mev used 
in this experiment, the deuterons lose, respec- 
tively, 0.07 and 0.05 Mev/cm or 0.20 and 0.14 
Mev in the effective thickness (2.8 air-cm) of the 
target. Hence no correction need be made here 
for the target thickness. There is a disadvantage 
in using this thick target scheme, for the alpha- 
particles coming into the chamber have residual 
ranges extending from zero to 2.8 air-cm and 
so fluctuations in over-all counter sensitivity 
give variations in counting rates. However, a 
change in counter sensitivity of 1 percent gave 
a change of less than 1 percent in the counting 
rate. Before each count, the amplifier sensitivity 
was set to within 1 percent of a fixed value and 
generally the statistical fluctuations in counts 
was about 1 percent so that this error does not 
affect the readings more than the statistical 
fluctuations of the counting. The advantage in 
this type of thick target is that the target is 
simple to prepare and there is little danger in 
burning it out during bombardment or breaking 
it during the subsequent measurements. 

Each bombardment was monitored with a 
current integrator to obtain absolute cross sec- 
tion. However, since the cyclotron was never 
considered to give an identical bombardment for 
different runs, the absolute cross section is not 











H. TATEL AND J. 








Residual Beam Current (arbitrary units) 





aR cm.-@ 49 $2 





L | j 


°o 


' 1 ! 
uev-—* 66 68 90 


Fic. 3. Plot of residual cu 
of deuteron beam of cyclotal 
7 through different thicknesses of 


- yn RANGE aluminum foil. 











40 so 3 70 
Thickness of Aluminum (mg/cm?) 


considered to be a reliable result although the 
trend should be fairly consistent. 

The energy of the deuteron beam was esti- 
mated visually by measuring the range in air of 
the beam, the mean range being taken as the 
point along the beam at which the brilliancy of 
an irradiated fluorescent screen started to dimin- 
ish rapidly. This measurement was made as a 
check on a residual range determination in which 
the residual beam current to a Faraday cup was 
measured as the beam passed through aluminum 
foils of various thickness, Fig. 3. The air equiva- 
lent of aluminum was taken as 1.52 milligrams 
per cm? per air-cm and corrected according to 
Livingston and Bethe. The effective maximum 
deuteron range is taken as 54 air-cm or 8.8 Mev. 
Although, the beam is not monoenergetic, the 
trend of the results should not be seriously 
affected as the cross section is changing so 
rapidly that for any experiment only the high 
energy tail is efficacious in making the transmu- 
tation. The straggling in energy for a 5- and 
9-Mev deuteron beam is about 0.1 and 0.16 Mev, 
respectively, while the approximate straggling 
of the beam in the cyclotron was 0.5 Mev. 


RESULTS 


Bismuth targets were bombarded with deu- 
terons from zero to 0.006 inch of aluminum in 





steps of 0.001 inch (6.6 mg/cm’). These steps 
decreased the mean energy of the beam by 
approximately 0.5 Mev per step from a maximum 
of 8.8 Mev. During bombardment, the beam 
current was kept to within 70 percent of a norm 
and the total quantity of deuterons and the time 
of bombardment recorded. After bombardment, 
the targets were mounted in front of a parallel 
plate ionization chamber and linear pulse ampli- 
fier arrangement. Counting rates of the emitted 
alpha-particles were measured at intervals for 
from five to ten days. All targets showed an 
initial rate, an increase, and then (if followed 
long enough) a decrease. The graphs in Figs. 1 
and 2 show that each target accurately follows 
the assumed pattern of formation of RaE and 
Po. The quantities of these end products, their 
ratios, and the ratios corrected for finite bom- 
bardment time are all shown in Table I. The 
absolute cross section is roughly computed as- 
suming that all the Bi nucleii in 2.8 air cm (10.6 
mg/cm?) are effective as targets and that the 
effective solid angle is defined by the circular 
1.5-cm opening of the ionization chamber 1 cm 
from the target under observation. For com- 
parison, the ratios from the two previous meas 
urements are also tabulated Razz! and Rear. 
These results are plotted as follows: cross section 
for formation of Po (ca) in Fig. 4 and the 
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corrected ratio R for finite bombardment time, 
of RaE/Po or cap/can in Fig. 5. 


SOURCES OF ERROR 






A major source of uncertainty in the measure- 
ments is the deuteron energy. There are two 
difficulties. First, the beam had a considerable 
spread as shown in Fig. 3 and second, the energy 
varied from time to time. This may be the 
reason for the disagreement of the two measure- 
ments of R with 5 mils of Al. Because so many 
measurements were taken on each sample, it is 
thought that the measurement of R on any 

particular sample has an over-all accuracy of 
10 percent at 6 mils to 3 percent at zero mils of 
Al absorber in the beam. 


DISCUSSION 





As has been pointed out, the d, p reaction 
may formally be described by either an OP 
process, I, or by the formation of a compound 
nucleus with subsequent proton mission, II. 
The reaction III is clearly of the latter type. In 
order to compare the OP process, I, with an 
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Fic. 4. Plot of cross section of d, m reaction as a function 
of the energy. The two curves are plotted for the Gamow- 
Condon-Guerney penetration function for small and large 
nuclear radius. Both curves arbitrarily fit at about 8 Mev. 
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reaction as a function of the deuteron bombarding energy. 
The curves are the ratio of Bethe’s Oppenheimer-Phillips 
penetration factor to the Gamow-Condon-Guerney pene- 
tration factor. The best fit is with the larger nuclear radius 
of 2.0A'X 10-" cm. 


“ordinary”’ process as III, it would be necessary 
for reaction II to be extremely improbable. How- 
ever, there are no experimental data by which 
reaction II] may be ruled out compared to 
reaction III. 

It is usually considered that when a heavy 
compound nucleus is formed, the competition 
of the various processes is such that the emission 
of a neutron is considerably more probable than 
emission of a proton and so reaction III would 
be more probable than reaction II. This is 
deduced by Bethe,‘ paragraph 79, in the following 
manner: For the ordinary type of reaction in 
which a compound nucleus is formed by deuteron 
bombardment with subsequent emission of neu- 
trons and protons: 


Can= Tr EaP dT, /T 
Cap = aXEsP JET, /T 


and 

(4) 
where £4 is the deuteron “sticking function,”’ Pz 
is the deuteron penetrability, /, is the “‘effective 
orbital momentum,” and T is the total width of 


*H. A. Bethe, Rev. Mod. Phys. 9, 69 (1937). 
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all levels of the compound nucleus in the energy 
region in question. All these quantities are func- 
tions of the incoming particles leaving I,, I'p, 
the partial widths of the states with emitted 
particles or p in the final residual nucleus. 


Hence, 
Can/Cap=T'n/T >». (3) 


These partial widths depend on the total number 
of available states and on the penetrabilities of 
the emitted particles both of which in turn 
depend upon available energy. The difference in 
this available energy is 


U,—U,=e.—0.8 Mev+By,, 


which is expression 40 from Bethe,’ where e_ is 
the energy of the 6-ray of RaE (in this particular 
case 2 Mev), the 0.8 Mev is the mass difference 
of neutron and hydrogen atom, and B, is the 
proton barrier height (5 Mev) whence 


U,— Up~6 Mev. 


Hence, there is considerably more energy avail- 
able for the neutron emission than for the proton 
emission and since the partial widths depend 
upon the energy available, it is seen that I, 
should be greater than I, and, therefore, from 
(3) the cross section for emission of the neutron 
should be greater than the cross section for the 
emission of the proton. While there is not enough 
knowledge on the details of these nucleii available, 
a probable estimate of this ratio may be obtained 
from statistical considerations. Bethe (41)* shows 
that the ratio of the neutron and proton partial 
widths is given by 


r,/T',=exp(U,— U,)/T. 


The nuclear temperature T for bismuth for this 
type of excitation will be approximately 1.3 Mev 
and so the ratio of the cross sections would be a 
rather large number e* which would indicate 
that reaction III is considerably more probable 
than reaction II and thus the experiment is 
measuring only reaction I and III. Assuming 
then that this is true, it should be noted that 
Bethe’s critiques* for observing the OP process 
are fulfilled in this reaction, i.e., the deuteron 
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bombarding energy is below the barrier he 
minus the deuteron binding energy (2.2 Mev) 
The barrier height of bismuth for deuterons ig 
9.8 and 14.3 Mev for the large and small nuclear 
radius, respectively. 

The only factor used here in the theoretical 
treatment of the comparison of the Oppenheimer. 
Phillips type of reaction to the ordinary type of 
reaction are the penetrability factor. These are 
also the most rapidly varying factors with respect 
to energy and it is usually considered that the 
ratio of the cross sections (y) is given by the 
ratio of the penetrabilities. Hence 


¥ =c0r/caca~ Por/Paca, 


where the subscript OP refers to the Oppen. 
heimer-Phillips penetrability function and GCG 
refers to the Gamow-Condon-Guerney type of 
penetrability for the deuteron. 

In fitting the theoretical y to the experiment, 
the most important parameter is the nuclear 
radius. Near the barrier height, this is especially 
important as the theory predicts that y=1 at 
the barrier height minus the binding energy of 
the deuteron. For a radius of 2.0A!X10-" em 
and 1.44!x10-" cm, this occurs for Z=83 at 
7.6 and 12.1 Mev, respectively. In Fig. 5, the 
data have a more slowly varying region in the 
neighborhood of 7.6 Mev which is what the 
theory predicts for the larger nuclear radius. 
The y's for these two values are plotted also in 
Fig. 5. The values for the smaller radius are 
taken from Volkoff’s® discussion of Bethe’s OP 
theory. The theoretical curves are fitted to the 
experimental values at 7.6 Mev. The trend of 
the results are in general agreement with an 
OP process of finite nuclear radius and the radius 
best matching these results is the larger 2.04! 
x<10—* cm. Of further interest, it is seen that 
d, n follows rather closely the Pace function and 
so the effect of the other factors which enter in 
the cross sections other than the penetration 
factors varies slowly compared to them. 

This investigation was made possible by 4 
grant from the Horace H. Rackham Fund. 


5G. M. Volkoff, Phys. Rev. 57, 866 (1940). 
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The slow neutron velocity spectrum of Au, In, Ta, W, 
Pt, and Zr have been investigated. The results for Au and 
In check the values previously obtained using the ‘‘old” 
system. The results for the gold level at (4.8+0.1)ev with 
oI = 600, and for the main indium level at (1.44+0.02) ev 
with ooI!*=210 agree with the earlier results. Additional 
indium levels were found at (3.8+0.2) ev with ooI'~120 
and at (8.6+0.4) ev with oI’ ~300. The cross-section curve 
for Ta below 1 ev is well matched by o=[(7.2+0.4) 
+(3.0+0.1)E-+]. Tantalum has levels at (4.1+0.1) ev 
with ool*~44, at (10.0+0.3) ev with ooI*~25, at (13+0.5) 
ev with ooI?~3, at (22+2) ev with opI*~18, at (37243) ev 
with ooI*~400, and other dips probably indicating mul- 
tiple strong levels near 100 ev, 300 ev, and at higher 
energies. The cross section curve for W below 1 ev is well 
matched by o=[(5.7+0.2)+(2.72+0.05)E-#]. There are 
levels at (4.0+0.1) ev with ooI*~13, at (7.4+0.2) ev with 


(Received November 4, 1946) 


ool*~5, at (18.0+0.5) ev with ooI*~3000 (perhaps mul- 
tiple), at (45+2) ev with ooI*~400 (perhaps multiple), at 
(180+20) ev with ooI*~10,000 if single (probably mul- 
tiple), and a dip at 1100 ev which is probably multiple. 
The results for Pt show pronounced crystal interference 
effects in the thermal region. The 1/2 line is probably 
o =[(12.040.3)+(1.0340.06)E-4], There are levels at 
(11.5+0.4) ev with ooI*~S5, at (18.241) ev with ooI*~30, 
and broad dips near 100 ev and near 1000 ev indicating the 
presence of several strong unresolved levels. The results for 
Zr below 0.6 ev are well matched by o=[(6.8+0.3) 
+ (0.74+0.10)E-*]. There are weak levels at (1.09+0,03) 
ev with opI'~5 and at (2.30+0.07) ev with opI'~8. A 
complex dip near 7 ev probably consists of two levels at 
Eo=(5.7+0.5) ev with ooI'~10 and at (7.6+0.4) ev with 
ool*~40. 








INTRODUCTION 





PT HIS is the second of a series of papers in 

which the results of measurements using 
the new Columbia University Slow Neutron 
Velocity Spectrometer are presented. A descrip- 
tion of the first system and a detailed analysis of 
operational factors is given in two previous 
papers." The description of the new system and 
the results for Cd, Ag, Sb, Mn, and Ir are given 
in the first paper of this series.* 

For convenience in interpreting these results 
a few of the factors involved in the operation of 
the apparatus and in the analyses applied to the 
data are listed below. A complete discussion of 
these factors is given in the three previous 


papers. 


*This document is based on work performed under 
Contract W-7405-Eng-50 with the Manhattan Project and 
Columbia University. It will also be published in the 
Manhattan Project Technical Series, Division III, as part 
of the contribution of the Physics Department, Columbia 
University. Publication assisted by the Ernest Kempton 


Adams Fund for Physical Research of Columbia University. 

1L. J. Rainwater and W. W. Havens, Jr., Phys. Rev. 
70, 136 (1946). 

*W. W. Havens, Jr. and L. J. Rainwater, Phys. Rev. 
70, 154 (1946). 

*L. J. Rainwater, W. W. Havens, Jr., C. S. Wu and 
J. R. Dunning, Phys. Rev. 71, 65 (1947). 
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OPERATION 


1. A burst of neutrons is produced at the 
source between t=0 and t=r, and these neutrons 
are detected 6.2 meters from the source between 
t=t, and t=t++7 to define the neutron energy in 
terms of its time of flight ¢, over the measured 
path length. 

2. Sixteen adjacent timing intervals are meas- 
ured simultaneously for each of two independent 
detectors, thus 32 measurements are made 
simultaneously. 

3. The detection intervals are spaced in time. 
an amount r equal to the width of the individual 
intervals. The resolution function of the instru- 
ment may be considered as triangular in shape 
with a base of approximate width 2r or twice the 
spacing between adjacent points. Actually the 
width of the resolution triangle includes other 
factors such as a term due to the counter length, 
and, below 0.3 ev, a term due to the diffusion 
time of the thermal neutrons leaving the source. 
In almost all cases, the resolution width may be 
taken as 2 to 3 times the spacing between 
adjacent points. 

4. The results of the measurements are pre- 
sented as plots of the slow neutron transmission 
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Fic. 1. The slow neutron transmission of 0.210 g/cm? 
of iF This curve shows the resonance dip in transmission 
at 4.8 ev. 


of the sample vs. the slow neutron time of flight 
in usec. /meter. Scales of energy are provided for 
convenience. The relation E = (72.3/t)? ev gives 
the neutron energy for a timing ¢ usec. /meter. 

5. The cross section, ¢, is related to the trans- 
mission, J, by the relation T=e—"’ where n 
=atoms/cm? of the element studied (if only the 
element is present). From this equation a cross 
section scale indicated on each curve is computed, 
and if o is a slowly varying function of ¢t, the 
measured value of o will correspond to the true 
value within the accuracy of each measurement. 
Near a resonance the effect of the resolution 
. width will affect the measured cross section ap- 
preciably, and for resonances above 1 ev, the 
value of ¢ calculated from the transmission curve 
is usually very much smaller than the true value 
of «. In particular the true value of oo at exact 
resonance may be from 2 to 1000 times larger 
than the directly measured value ¢ max., from 
the transmission dip. 

The quantity oI can usually be determined 
to within a factor of two by a method previously 
described. * This value is usually given in addi- 
tion to the position (energy) of the resonance 
level. In some cases where the condition noo>1 
is not satisfied, a rough value of oI’ may be 
given. Although these values are only roughly 
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quantitative, any information concerning \the 
“strength” of the levels, was thought to be of 
value. 

6. For many elements the cross section curye 
for energies lower than the first resonance may 
be closely matched by the sum of a 1/v term and 
a constant term. The equation for the Cross 
section in this form is given where possible, 

7. The units of energy £, cross section ¢ and 
the oI and oI values will usually not be stated 
specifically. It is to be understood that energy js 
expressed in units of electron volts, ev, and crogg 
sections in units of X10-* cm?/atoms. T is the 
full level width (in ev) at half-maximum in the 
Breit-Wigner formula. 

Unless otherwise specified all cross sections are 
averages for the natural element and are not 
assigned to particular isotopes. 


GOLD 


The resonance in gold near 4.8 ev has previ- 
ously been studied with the first Columbia Neu- 
tron Spectrometer.? The earlier measurements 
in gold have been checked with the new system. 
In Fig. 1, the transmission curve is shown for a 
0.210 g/cm? gold foil. The position of the reso. 
nance at (4.8+0.1) ev is in agreement with the 
previous measurement. The earlier value oI* 
=600 is also in good agreement with Fig. 1. 
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Fic. 2. The slow neutron transmission of 0.193 g/cm* 
of indium. This curve illustrates the three resonances at 
8.6 ev, 3.8 ev, and 1.44 ev. 
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This level must be assigned to the Au'®’ for gold 
has only one isotope. The half-life of Au'®* 
resulting from the capture of a neutron by Au!*? 
is 2.7 days.*® If ['~0.1 ev, then oo~60,000, 
which is 20 times larger than the measured 
o max. 

INDIUM 


The indium resonance was also investigated 
with the previous system.? The main level was 
located at 1.44 ev with o oI?=210. The same 
sample has been remeasured over a wider range 
with better resolution using the new system with 
the results shown in Fig. 2. Two additional 
resonances at 3.8 ev and 8.6 ev have been located 
with ooI'~120 for the 3.8 ev level and opI'~300 
for the 8.6 ev level. For ['~0.1 ev this corre- 
sponds to o»>~1200 and 3000 respectively. There 
is no indication of other strong levels. (Since the 
sample is thin, only strong levels would be 
observed.) 

The results Ey) = 1.44 ev and ooI? = 210, for the 
main level, agree with the results of Fig. 2. The 
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Fic. 3. The slow neutron transmission of 0.03865 g/cm? 
of indium near the 1.44 ev resonance. The solid line is a 
theoretical Breit-Wigner curve with Ey = 1.44 ev, ool? =210, 
and !'=0.09 ev, corrected for the resolution width of the 
apparatus. 
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‘Amaldi, D'Agostino, Fermi, Ponticorvo, Rasetti, and 
Segré, Proc. Roy. Soc. London A149, 522 (1935). 
® McMillan, man, and Ruben, Phys. Rev. 52, 375 


(1937); Pool, Cork, and Thornton, Phys. Rev. 52, 239 
(1937); Dzelepow and Constantinov, Comptes Rendus 
U.S.S.R. 30, 701 (1941). 
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1.44 ev level was also studied using a 38.65 
mg/cm? sample with the results shown in Fig.3, 
where the resonance is more nearly resolved since 
o max. = 11,000. If [T=0.09 ev, then op = 26,000 
which is within a factor of 24 of the measured 
value of o max. 

Indium has two isotopes In"* (4.5 percent) and 
In™5 (95.5 percent). In‘ results from the capture 
of a neutron by In"™*, which has a half-life of 48 
days.* The 13-sec. and 54-min. In™* are isomers.*’ 
The level strengths for the responsible isotopes 
must be correspondingly corrected. The final 
results may be expressed : 

1. The first resonance has Ey= (1.44+0.02) ev 
ool*~210. (If T'=0.09 ev, then oo~26,000.) 
Value omax = 11,000 was obtained. 

2. The second resonance has Ey) = (3.8++0.2) ev 
ool ~120. (If [~0.1 ev, then oo>~1200.) 

3. The third resonance has Ey =(8.6+0.4) ev 
with ool ~300. (If ['~0.1 ev, then oo~3000.) 
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Fic. 4. The slow neutron transmission of 22.4 g/cm? 
of tantalum. This curve illustrates the 1/v type cross 
section in, the thermal region with (o=7.243.0F-4) and 
also shows the resonance levels at the higher energies. 


* Lawson and Cork, Phys. Rev. 52, 531 (1937); Mitchell, 
Phys. Rev. 53, 269 (1938). 

7 Lawson and Cork, Phys. Rev. 52, 531 (1937); Mitchell 
and Langer, Phys. Rev. 53, 505 (1938). 
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Fis. 5. The slow neutron transmission of tantalum. The upper curve shows the transmission for a 9.98 g/cm? sample 
with a resolution width of 3.6 ysec./meter. The lower curve shows the transmission for a 22.44 g/cm? sample with a 
resolution width of 1.8 ysec./meter. The use of the thicker sample with higher resolution increases the detail that can 


be obtained. 


TANTALUM 


The slow neutron transmission of tantalum 
was investigated using plates of tantalum metal 
obtained from the University of Chicago Metal- 
lurgical Laboratories. ‘A spectrographic analysis 
of the samples showed less than 5 percent Cb; 
less than 0.05 percent Zr and Al; less than 0.005 
percent Ti, Mo, Fe, Co, Os, Ag, Zn, Cd, Sn, Pb, 
Sb and traces of Cr, Ni, Ru, Rh, Pd, Pt, and Hg. 
For all other metals the prominent lines were not 
apparent, therefore the impurity was probably 
less than 0.05 percent for arc insensitive ma- 
terials and less than 10-5 percent for arc sensitive 
materials. Samples of 9.98 g/cm? and 22.44 
g/cm? thickness were studied. The results of 
these measurements are shown in Figs. 4 and 5. 

The general features of the curve over a wide 
energy range are shown in Fig. 4. The results 
for E<1 ev are well matched by the relation 


o=(7.2+3.0E-) and several levels are shown 
at higher energies. The high energy region was 
first investigated using a 9.98 g/cm? sample 
with about 3.2 yusec./meter resolution width. 
The results shown in the upper curve of Fig. 5 
show the presence of strong levels at 4.1 ev, 10 
ev, and 37 ev with indication of several other 
levels near 13 ev, 22 ev, and above 100 ev. 

In view of the complex structure revealed, the 
portion of the curve above 4 ev was remeasured 
using a 22.44 g/cm* sample with twice the 
resolution. The results are shown in the lower 
curve of Fig. 5. All of the structure in the upper 
curve is repeated with better resolution in the 
lower curve, thus conclusively demonstrating the 
reality of the effects obtained. 

An analysis of Fig. 5 can be made by using the 
position of the 1/v line and the constant term for 
the thick sample as determined in Fig. 4. If the 
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"curve is assumed to consist of a series of reso- 


nance dips there are strong levels at 4.1 ev, 10 ev, 
and 37 ev, weaker levels at 13 ev, 22 ev, and 

bably multiple levels near 100 ev and 300 ev, 
and in the higher energy region. Even though 
the levels below 100 ev may be multiple, an 
analysis to determine oI? has been made as- 
suming them to be single. The curve cannot be 
divided into a series of separate dips with any 
degree of accuracy near 13 ev and 22 ev. There- 
fore, the values of ooI? for these levels may be 
uncertain by more than the customary factor of 
2. 

The results for tantalum are of particular 
importance since all of the levels must be ascribed 
to the single Ta'* isotope. This is the first time 
in the knowledge of the authors that so many 
closely spaced levels have been experimentally 
demonstrated to belong to a single nucleus, 
although such effects are expected from the 
theory of level spacing.* The strength of the 
levels makes it unlikely that slight sample im- 
purities can be responsible for any of the levels 
observed. The activity induced in tantalum by 
the (, y) reaction has a half-life of 97 days.°® 

The results for tantalum may be summarized : 

1. Below 1 ev the data are well matched by 

the 1/v relation 
o=((7.2+0.4)+(3.040.1)E]. 

2. The first and main level occurs at Ey = (4.1 

+0.1) ev 
ool? ~ 44. 
(If the 1/v slope were entirely attributed to 
a level at 4.1 ev then oo’? =96 is required.) 
3. The second strong level is at Ey»=(10.0 
+0.3) ev 
ool? ~25. 
4. The third level is at E> =(13+0.5) ev 
ool? ~3: 
5. The fourth level is at Ey =(22+2) ev 
ool? ~ 18. 
6. The fifth level is at Ey) =(37+3) ev 
ool”? ~400. 
7. There are indications of probably more than 
one strong level near 100 ev, near 300 ev 
and at higher energies: 


*H. A. Bethe, Phys. Rev. 50, 332 (1936). 

*Oldenberg, Phys. Rev. 53, 35 (1938); Fomin and 
Houtermans, Physik: Zeits. Sowjetunion 9, 273 (1936); 
Houtermans, Naturwiss. 28, 578 (1940). 
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TUNGSTEN 


The slow neutron transmission of tungsten was 
studied using rolled tungsten plates supplied by 
the General Electric Company, Cleveland Wire 
Works. Six plates were used for the 29.2 g/cm? 
sample (Fig. 6) and single plates of 5.06 and 
4.90 g/cm? for the high energy region (Figs. 7 
and 8). A spectrographic analysis of the sample 
showed less than 1 percent As; 0.5 percent Bi; 
0.05 percent Zr, V, Zn, Cd, and Si; 0.005 percent 
Be, Fe, Pd, Au, Ge; traces of La, Ru, Rh, Cu, 
Hg, Al, Tl, Sn, Pb. For all other metals the 
prominent lines were not apparent. 

The entire region was first investigated with 
low resolution using the thick sample as shown 
in Fig. 6. Below 1 ev the data are well matched 
by the relation ¢=(5.7+2.72E->), and a dip is 
obtained showing the presence of one or more 
levels above 3 ev. 

The region above 2 ev was next investigated 
using the 4.90 g/cm? sample and about 3.2 
usec./meter resolution width as shown in Fig. 7. 
This shows the presence of a strong level at 18 ev, 
a strong level at 4.0 ev, and weaker levels near 
7.4 ev, 45 ev, 180 ev, and 1100 ev. Since the 
curve was so complex, it was considered of 
interest to investigate the region above 4 ev 
with higher resolution. A 5.1-g/cm? sample was 
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Fic. 6. The slow neutron transmission of 29.2 g/cm? 
of tungsten. This curve shows the 1/v type cross section 
for tungsten in the thermal region with ¢=(5.7+2.72E-4) 
and also illustrates the effect of resonances in the high 
energy region when broad resolution is used. 
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used for this purpose with the results shown in 
Fig. 8. As in the case of tantalum all of the 
effects observed with the thicker sample and 
higher resolution are repeated with greater detail 
thus giving an excellent check on the measure- 
ments. The level at 7.4 ev appears again as a 
sharp dip, a more definite dip is now obtained at 
45 ev; the level (or levels) near 180 ev is more 
pronounced; and the dip near 1100 ev is more 
definite. The shape of the curve near 10 ev may 
suggest the presence of a weak level, but the 
magnitude of the effect is so small that the 
presence of the level here is questionable. 

Tungsten has 5 stable isotopes of which 4 have 
(within a factor of 2) approximately equal abun- 
dances. The single odd isotope W'®* of 17.3 
percent abundance should be responsible for the 
strong levels. These are two well-established 
activities with half-lives of 77-d and 24.1-hr. 
induced in tungsten by (m, y) reactions.*!° 

The possibility has been suggested that the 
4.0 ev level in Fig. 7 is caused by a small tantalum 
impurity and represents the 4.1 ev tantalum 
level. A comparison of Figs. 5 and 7, however, 
indicates a small but definite difference in level 
position. In addition a calculation of the strength 


10 Minakawa, Phys. Rev. 57, 1189 (1940); Fajans and 
Sullivan, ~~ Rev. 58, 276 (1940) ; McLennan, Grimmet, 
and Read, Na 


ture 135, 147 (1935). 
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ool?~13 for the level would require that the 
sample be about 30 percent tantalum by weight. 
The results for tungsten may be summarized 
as follows: 
1. Below 1 ev the data are well matched by 
the 1/v relation 
o=[(5.7+40.2) +(2.72+0.05) E+]. 
2. The first strong level is at Ey = (4.0+0.1) ev 
ool? ~ 13. 
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Fic. 8. The slow neutron transmission of 5.06 g/cm* 
of tungsten. The resolution used in this curve is twice as 
high as that used for Fig. 7 and confirms dips at 45 ev and 
1100 ev. This data also gives greater detail for the levels 
at 7.4 ev and 180 ev. . 
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Fic. 9. The slow neutron transmission of 25.8 g/cm? of platinum above 2.2 ev. This curve shows trans- 
mission dips for platinum at 11.5 ev, 18.2 ev, 100 ev, and 1000 ev with the upper transmission dips not 
completely resolved. 





























































If the W'® isotope is responsible, then 5. The fourth level is at Ey = (45+2) ev 









_ ool®~75 for this isotope. ool? ~400. (Perhaps multiple.) 
3. The second (weaker) level is at Ey=(7.4 6. The fifth dip is at Ey =(180+20) ev 
+0.2) ev ool? ~ 10,000 if it is due to a single level. 
ool? ~S. 7. Another dip near 1100 ev is probably 
4. The third and strongest level is at Eo multiple. 
=(18.0+0.5) ev PLATINUM 
ool ~ 3000. A sample containing 25.8 g/cm? of metallic 






This may be multiple as the curve in Fig. 8 platinum was used for the slow neutron trans- 
recovers more slowly on the high energy - mission measurements..A spectrographic analysis 
side of the resonance than would be ex- of the sample showed less than 0.5 percent Ca, 
pected for a single level. Al; 0.05 percent Fe, Rh, Pd, Cu, Au; 0.005 
percent Ag and Si; traces of Cr, Ru, and In. 
For all other metals the prominent lines were 

; not apparent. 
The results of the measurements with the 
: highest resolution for energies above 2.2 ev are 
shown in Fig. 9. This shows that there is a 
* moderately strong level at 11.5 ev, a weaker 
ee level at 18.2 ev, a broad dip probably indicating 
Fic. 10. The slow neutron transmission of 25.8 g/cm? several strong levels from 70 ev to 200 ev, and 
of platinum from 0.7 ev to 3 ev. This curve shows a another broad dip near 1000 ev again suggesting 


transmission dip at about 1.25 ev which is probably caused 
by a rhodium impurity. several levels. 





























The region was extended with somewhat lower 
resolution as shown in Fig. 10. This shows a 
shallow dip at about 1.25 ev which is probably 
because of a slight rhodium impurity. 

The results of the measurements with wide 
resolution from 0.009 ev to 0.4 ev are shown in 
Fig. 11. A 1/y slope o=(12.0+1.03E-) is ob- 
tained for the first part of the curve while 
definite crystal interference effects are obtained 
for lower energies. The main discontinuities cor- 
respond to Bragg reflection limits at neutron 
wave-lengths of 1.6A and 2.4A. 

Platinum has 5 stable isotopes of which the 
isotope Pt! (35.3 percent abundance) is the 
only odd isotope and, according to theory, should 
be responsible for the strong levels.* There are 
only three activities of 18-hr. 3.3-d and 31-min. 
associated with the (m, y) reactions in plati- 
num. 

The results for platinum may be summarized 
as follows: 

1. The 1/v slope (0.04 to 0.8 ev) is 

o=((12.0+0.3) +(1.03+0.06) E+]. 

2. There are definite crystal effects for E 

<0.035 ev with transmission discontinuities 
for the neutron wave-lengths \=1.6A and 
2.4A. 
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Fic. 11. The slow neutron transmission of 25.8 g/cm?* 
< — in the thermal region. This curve illustrates 
e 1/v 


+1.03E-4) and the crystal interference effects observed at 
the lower energies. 


1% McMillan, Kaman, and Ruben, Phys. Rev. 52, 375 
(1937); Sherr, Bainbridge, and Anderson, Phys. Rev. 60, 
' 473 (1941); Pool, Cork, and Thornton, Phys. Rev. 52, 239 
(938- McLennan, Grimmet, and Read, Nature 135, 147 


pe curve in the thermal region with o=(12.0 . 
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3. There is a level at Eo=(11.5+0.4) ey ‘ 
ool? ~55. 

4. There is a level at Eyo=(18.2+1) ey 
ool? ~ 30. 

5. There are broad dips at 100 and 1000 ey 
indicating the presence of strong unresolved 
levels. 


ZIRCONIUM 


A sample consisting of 6.45 g/cm? of Zr0, 
was used to investigate the slow neutron trans. 
mission of zirconium. A spectrographic analysis 
of the sample showed less than 5 percent Ca; 
0.5 percent Zn, Si, Sr; 0.05 percent Mg, Ti, Ma, 
Fe, Cu, Al; 0.005 percent Be, B; traces of Sc, 
Y, V, Cr, Mo, Ni, Sn, Pb, Sb. For all other 
metals the prominent lines were not apparent, 


_In all cases a constant cross section of 4.1 has 


been deducted for the oxygen.§ 

The results of measurements with wide resoly- 
tion over the entire energy region are shown in 
Fig. 12. This shows that levels are present above 
1 ev and that the results for lower energies are 
well matched by the 1/v curve ¢ = (6.8+-0.74E7>), 

The region above 0.6 ev was investigated with 
better resolution as shown in Fig. 13. The dips 
indicated that there were single levels at 1.09 ey 
and 2.3 ev and perhaps a double level near 7 ey, 
The region 4 to 10 ev was then investigated 
with highest resolution as shown, and the dip 
seems to be resolved into two dips at 7.6 ev and 
4.7 ev. There are no indications of other reso- 
nances. 

Zirconium has 5 stable isotopes of which Zr" 
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Fic. 12. The slow neutron transmission of 6.45 g/cm* 
of ZrO:. This curve shows the 1/v type cross section for 
Zr in the thermal region with o=(6.8+0.74E-4) after 4.1 


has been deducted for the oxygen cross section. 
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‘= Fic. 13. The slow neutron cross section of 6.45 g/cm* of ZrO». This curve shows the transmission dips at 1.09 ev 
Ove and 2.3 ev and the odd shaped dip in the vicinity of 7 ev which might be the sum of two dips at 5.7 ev and another 
are at 7.6 ev. 
4), , 
vith (11.5 percent abundance) is the only odd isotope 3. The second level is at Eo=(2.30+0.07) ev 
dips and, therefore, should be responsible for the ool ~8. 
~ stronger levels. The weakness of the levels 4. There is bulge on the next dip indicating a 
a suggests that one or more may be owing to slight level at Eo=(5.740.5) ev with ooI'~10 
sted impurities in the sample. The induced activities (uncertain). 
dip in zirconium by (m, y) reactions are not yet 5. The next level is at Ey =(7.6+0.4) ev 
a definitely assigned. ool’ ~ 40. 
on The results for Zr may be summarized as We take pleasure in acknowledging our in- 
follows : debtedness to Professor J. R. Dunning for his 
zt 1. Below 0.6 ev the results are well matched valuable suggestions and stimulating discussions. 
by the 1/v curve We wish to thank Miss Miriam Levin, who 
o=[(6.8+0.3) +(0.74+0.10) £4]. assisted with the numerous calculations involved 
2. There is a partially resolved level at Ey in this paper, and the other members of the 
=(1.09+0.03) ev cyclotron staff who aided with these measure- 






ool ~S. ments. 
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Slow Neutron Velocity Spectrometer Studies. III. ‘ 
I, Os, Co, T1, Cb, Ge* 
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Columbia University, Pupin Physics Laboratories, New York, New York 
(Received November 8, 1946) 


The 1/v thermal line for iodine is o=[(3.8+0.2) 
+(1.12+0.05)E-4]. Crystal interference effects are ap- 
parent at low energies. There is a resonance at (20.6+0.4) 
ev with ooI*~4. A flat-bottomed dip at 37 ev suggests two 
levels at (3242) ev and (42+2) ev with approximately 
equal values of ooI*~135. The 1/v thermal line for osmium 
is o=[(15+0.4)+(2.7+0.1)E-4] with crystal interference 
effects at low energies. There are levels at (6.50.3) ev 
with ool? ~10, at (8.8+0.3) ev with oof?~35, at (20+1) ev 
with oI*~25, probable levels at (28+1.5) ev with 
ool? ~8, at (42+2) ev with ooI*~10, at 84 ev and several 
unresolved levels above 84 ev. The 1/v thermal line for 
cobalt is o=[(6.7+0.3)+(6.4+40.15)E-4]. There is an 
unusually strong dip at (115+5) ev with ooI?~200,000 


(if single) and indications of other strong unresolved levels 
in the region of 1000 to 10,000 ev. The 1/v thermal line for 
thallium is _ probably [o= (9.70.2) + (0.640.2)E-4), 
There are strong crystal interference effects at low energies, 
There is a strong dip in the transmission curve at 270 ey 
and another broad dip near 1100 ev. Columbium shows 
only a slight 1/v effect o =[(6.40.2) + (0.10+0.04) B+) 
with strong crystal interference effects at lower energies, 
Only a slight dip at 4.1 ev attributed to a slight Ta impu. 
rity was found in the higher energy region. Germanium 
shows no 1/v thermal slope because of the strong crystal 
interference effects. There is a (perhaps multiple) level at 
95 ev with ooI?*~800. 





INTRODUCTION 


HIS is the third of a series of papers in 
which the results of investigations using 
the Columbia University Slow Neutron Velocity 
Spectrometer are presented. A detailed analysis 
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Fic. 1. The slow neutron transmission of 18.86 g/cm? 
of iodine. This curve indicates the presence of one or more 
resonance levels at high energies. From 0.04 to 3 ev the 
results are well matched by the relation ¢ = (3.8+1.12E>4). 


* This document is based on work performed under 
Contract W-7405-Eng-50 with the Manhattan Project 
and Columbia srg oe It will also be published in the 
Manhattan Project Technical Series, Book III, Division 
VIII, as part of the contribution of the Physics Depart- 
ment, Columbia University. Publication assisted by the 
Ernest Kempton Adams Fund for Physical Research of 
Columbia University. 


of the factors involved in the operation of the 
instrument and the interpretation of the data is 
given in previous papers.’ In presenting the 
results of the measurements for I, Os, Co, TI, 
Cb, and Ge the conventions previously adopted*‘ 
are used. 
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Fic. 2. The slow neutron transmission of 18.86 g/cm* 
of iodine. This curve shows a strong resonance dip near 
37 ev with indication of another unresolved level near 20 ev. 


1L. J. Rainwater and W. W. Havens, Jr., Phys. Rev. 
70, 136 (1946). 

2W. W. Havens, Jr. and L. J. Rainwater, Phys. Rev. 
76, 136 (1946). 

3L. J. Rainwater, W. W. Havens. Jr., C. S. Wu, and 
J. R. Dunning, Phys. Rev. 71, 65 (1947). 

4W. W. Havens, C. S. Wu, L. J. Rainwater, and C. L. 
Meaker, Phys. Rev. 71, 165 (1947). 
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SLOW NEUTRON 


Fic. 3. The slow neutron transmis- 
sion of 18.86 g/cm? of iodine. This 
curve was taken with the highest reso- 
lution of the spectrometer. The reso- 
nance at 20.6 ev is now clearly 
separated from the higher energy dip. 
The shape of the resonance dip at 37 
ev suggests the presence of at least two 
levels probably at 32 ev and 42 ev. 


RESULTS 
Iodine 

The slow neutron transmission of iodine was 
studied using an 18.86 g/cm? sample of ana- 
lytical grade iodine packed in a thin-walled 
aluminum container. This sample was prepared 
by the Metallurgical Laboratory of the Uni- 
versity of Chicago. 

The results of the measurements using broad 
resolution are shown in Fig. 1. There is a dip at 
high energies indicating the presence of onethe 
more resonance levels. From 0.04 to 3 ev or 
results are well matched by the relation ¢ 
=[(3.8+1.12)E+]. Below 0.04 ev there is evi- 
dence of crystal interference effects with discon- 
tinuities in the transmission occurring at neutron 
wave-lengths of about 1.46A and 2.37A. 

The energy region above 2.2 ev was next 
investigated with higher resolution as shown in 
Fig. 2. The main resonance dip is seen to occur 
near 37 ev with indication of another unresolved 
level near 20 ev. The lower transmission above 
200 ev may indicate the presence of many 
unresolved levels. 

To obtain more reliable information concern- 
ing the resonances, the region above 10 ev was 
investigated using the highest resolution possible 
with the results shown in Fig. 3. The resonance 
at 20.6 ev is now clearly separated from the 
higher energy dip. The dip at 37 ev now appears 
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flat on the bottom where the transmission is’ 
considerably above zero. The shape of this reso- 
nance absorption curve would not be obtained 
if only a single resonance were present; therefore 
the presence of at least two levels, probably at 
32 ev and 42 ev, seems reasonable. In the 
analysis equal oI? values have been assumed 
for levels at 32 and 42 ev. The curve above 50 ev 
shows small irregularities that may represent 
weak levels near 80 ev, 200 ev, and 600 ev, but 
these are very uncertain. A straight line is also 
shown to indicate that these levels may be due 
to spurious effects. 

Since iodine has only one isotope I’ all the 
neutron absorption effects can be assigned to 
this one isotope. The slow neutron induced 
activity in iodine is from I*5* which has a 
half-life of 24.99 min. The results for iodine may 
be summarized as follows: 

1. The 1/9 thermal line is given by 

o=((3.8+0.2) +(1.12+0.05) £4]. 

2. There are crystal interference effects below 
0.04 ev with discontinuities in the trans- 
mission curve for neutron wave-lengths of 
about 1.46 and 2.37A. 

3. The first resonance is at Ey = (20.6+0.4) ev 
ool? ~4. 


5’ Amaldi, D'Agostino, Fermi, Ponticorvo, Rasetti, and 
Segré, Proc. Roy. Soc. London A149, 522 (1935). 
* Tape and Cork, Phys. Rev. 53, 676 (1938). 
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4. There is evidence that the next dip is due 
to two levels at (32+2) ev and (42+2) ev 
with oI? values~ 135. 
5. There are doubtful indication of weak levels 
near 80 ev, 200 ev, and 600 ev. 
Osmium 
The slow neutron transmission of osmium was 
studied using a 17.3 g/cm? sample of metal 
powder. A spectrographic analysis of the sample 
showed less than 1 percent Si; 0.5 percent Ca, 
Ti, Fe, Al, Bi; 0.05 percent Cr, Mn, Pt, Cu, Au, 
Pb; 0.005 percent Ag; traces of Zr, V, Co, Ni, 
Pa, Pd, Ir, Sn, and P. For all other metals the 
prominent lines were not apparent; therefore the 
impurity was probably less than 0.05 percent for 
arc insensitive materials and less than 10-5 
percent for arc sensitive materials. The sample 
was originally prepared at the Bureau of Stand- 
ards where it was heated in purified helium for 
two hours at 800°C. It’ was then compressed 
into a cell under a load of 20 tons. The sample 
was later opened, reduced at 120°C and recap- 
suled. The transmission of this sample was 
investigated using low resolution over the full 
timing range as shown in Fig. 4. The best 1/v 
line (0.1 ev to 5 ev) is c=(15+2.7E~). Below 
0.1 ev there are evidences of crystal interference 
effects particularly near 0.028 ev and 0.015 ev 
corresponding to neutron wave-lengths of 1.75 
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Fic. 4. The slow neutron transmission of 17.3 g/cm? 
of osmium. The best 1/v line (0.1 ev to 5 ev) is e=(15 
+2.7E~4). Crystal interference effects show below 0.1 ev. 
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and 2.36A. In addition there is a dip at highér 
energies. 

The sample was next investigated using a 
6.4 usec. /meter resolution width which indicated 
a strong dip at 8 ev and a smaller one near 80 ey, 
The region above 4 ev was then investigated 
with a 3.2 ywsec./meter resolution width. The dip 
near 8 ev was partially resolved into two levels 
near 9 ev and 6 ev as in the case of the 20 ey 
level for iodine in Fig. 2. The results of this set 
of measurements showed that it would be 
profitable to use better resolution. Therefore 
the region above 4 ev was remeasured using the 
highest resolution with the results shown ip 
Fig. 5. The peaks at 6.5 ev and 8.8 ev are here 
clearly resolved, and the peak at 20 ev is evident. 
The curve above 20 ev stays well below the line 
o=15 representing the constant term in Fig. 4, 

The dips at 84 ev and 42 ev are small, but the 
fact that they are probably real is indicated by 
the fact that they repeated on the several sepa- 
rate portions of the data which were combined 
for Fig. 5. The dip at 28 ev is less certain although 
it is also repeated in separate portions of the 
data. 

Using the line ¢=15 (from Fig. 4) as the value 
of the transmission due to scattering alone, the 
shape of the curve suggests that one or more 
levels are present at 28 ev, 42 ev, and 84 ev and 
that there are unresolved levels above 84 ey. 
The rise in transmission near zero time of flight 
is probably because of a decrease in the scattering 
cross section at higher energies. 

Osmium has 7 stable isotopes to which the 
levels may be ascribed. Os!** is the most abun- 
dant (16.1 percent) odd isotope and thus may be 
responsible for the strong levels.’ The activities 
assigned to the (m, y) reactions in osmium are 
the 32-hr. Os'™ and the 17-day Os!%*.8 

The results for Os may be summarized as 
follows: 

1. The 1/v line in the thermal region is well 
represented by the equation: 
o=((15+0.4)+(2.7+0.1) E+]. 

2. There is evidence of crystal interference 
effects at lower energies particularly for 

7H. A. Bethe, Phys. Rev. 332, 50 (1936). 

®Kurtchatow, Latychew, Nemenow, and_Selinow, 
Physik. Zeits. Sowjetunion 8, 589 (1935); Sooke 


Friedlander, Phys. Rev. 59, 400 (1941) ; Zingg, Helv. 
Acta 13, 219 (1940). 
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Fic. 5. The slow neutron transmission of 17.3 g/cm? of osmium. This curve was 
taken with the highest resolution of the spectrometer. 


neutron wave-lengths of about 1.75 and 


2.36A. 

3. There is a level at Eo = (6.5+0.3) ev 
ool? ~ 10. 

4. There is a level at Ey =(8.8+0.3) ev 
ool? ~35. 


5. There is a level at Ey = (20+1) ev ool? ~25. 

6. There is probably a level at (28+1.5) ev 
ool? ~8 (to within a factor of 10). 

7. There is a level at Ey = (4242) ev ooI®~10 
(to within a factor of 5). 

8. There is a level at Ey=(84+6) ev with 
other unresolved levels above this value. 


Cobalt 


The slow neutron transmission of cobalt was 
investigated using a 12.4 g/cm* sample of 
cobalt metal powder. Since this was too thick 
for the resonance region measurements, a thinner 
sample of 1.06 g/cm? was also employed. A 
spectrographic analysis of the sample showed 
less than 0.5 percent P, Si, Mn; 0.05 percent Mg, 
Ca, Fe, Cu, B, Al; 0.005 percent Ni, As; traces 
of Be, Cr, Mo, Ag, In, Ge, Sn, and Pb. For all 
other metals the prominent lines were not appar- 
ent. 

The results of the measurements using broad 





resolution are shown in Fig. 6. Because of the 
sample thickness the region below 0.15 ev was 
not studied. The results below 5 ev are well 
matched by the 1/2 line. ¢=(6.7+64E-), and 
a dip is seen at higher energies. 

The sample was next measured with good 
resolution above 0.5 ev as shown in Fig. 7. The 
same 1/v line (determined ,from Fig. 6) is used 


05 
0.4 


TRANSMISSION 
° 
~m 


° 


NEUTRON ENERGY IN EV ~~ 


wo 5 2 1 Q5 O38 O2 O15 





NEUTRON TIME OF FLIGHT MIC ROSECONDS/METER 
Fic. 6. The slow neutron transmission of 12.4 g/cm* 

of cobalt. The results below 5 ev are well matched by the 

1/v line ¢ = (6.7+6.4E-4). A dip is seen at higher energies. 
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Fic. 7. The slow neutron transmission of 12.4 g/cm? 
of cobalt. This figure shows a strong dip near 100 ev and 
another small dip at zero time of flight. 


with excellent agreement. This figure shows a ' 


strong dip near 100 ev and another small dip at 
zero time of flight. 

The region above 50 ev was measured using 
the thinner sample with maximum resolution as 
shown in Fig. 8. The shape of the dip at 115 ev 
is approximately the shape that would be ex- 
pected for a single strong level at this point 
although this dip may be due to multiple unre- 
solved levels. The cross section indicates that 
_ in the region from 1000 to 10,000 ev there are 
probably several strong unresolved levels. 

The level at 115 ev if single, is an exceptionally 
strong level with ooIf!?~200,000 as determined 
from the magnitude of the resonance dip. If the 
1/v slope is assumed to be entirely attributed to a 
single level at 115 ev a value of oopI?~30,000 is 
obtained. This is much smaller than the other 
value but is still extremely large for a single 
resonance level. The lower value calculated from 
the 1/v slope may indicate that the dip at 115 ev 
is multiple or may represent the interference 
effects of the strong levels above 1000 ev on the 
thermal cross section. Because of the observed 
strength of levels in cobalt they are all probably 
caused by the main Co* isotope of 99.83 percent 
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Fic. 8. The slow neutron transmission of 1.06 g/cm? 


of cobalt. This curve was taken by using a thinner sample 
and maximum resolution. 


abundance. Co* ® induced by the (, y) reaction 
has two isomeric states with half-lives of 5.3 yr. 
and 10.7 min. 

The results for cobalt may be summarized as 

follows: 

1. The results below 5 ev are well matched by 
the 1/v relation 
o=[(6.7+40.3)+(6.4+0.15)E]. 

2. There is a strong dip indicating a resonance 
at Ep=(115+5) ev oo!?~200,000 if the dip 
is because of a single strong resonance. If 
the 1/v slope is attributed to a single level 
at 115 ev then ooI?~30,000. The dip may 
be owing to several unresolved strong levels. 

3. There is evidence of other strong unresolved 
levels in the region of 1000 to 10,000 ev. 


Thallium 
The slow neutron transmission of thallium was 
studied using a sample containing 52.09 g/cm’ 
of the material. A spectrographic analysis of the 


sample showed less than 0.05 percent Ca, Pb; 
0.005 percent Fe, Cu, Al; traces of Mg, Cr, Ma, 


* Livingood and Seaborg, Phys. Rev. 60, 913 (1941); 
Deutsch and Elliott, Phys. Rev. 62, 558 (1942). 
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Cd, Si. For all other metals the prominent lines 
were not apparent. 

The, region above 0.015 ev was investigated 
using broad resolution as shown in Fig. 9. This 
shows that crystal interference effects are strong 
enough to make the determination of the 1/v 
slope difficult. The region above 0.5 ev was 
remeasured using better resolution as shown in 
Fig. 10. From this the 1/v line may be expressed 
as o=[(9.7+0.2)+(0.6+0.2)E+]. This curve 
also shows the presence of a strong dip at several 
hundred volts energy with a rise in transmission 
well above the line ¢=9.7. This rise in trans- 
mission near zero time of flight probably repre- 
sents a decrease in the scattering cross section 
at still higher energies (>50,000 ev for example). 

The high energy region above 50 ev was next 
investigated with maximum resolution as shown 
in Fig. 11. The shape of the curve suggests that 
there is a strong level at 270 ev and another 
strong level near 1100 ev although a number of 
unresolved levels may be present to give the 
observed curve. A rough analysis of the strength 
of the level at 270 ev, if single, indicates that 
ool?~ 20,000, a very strong level. If the 1/v 
slope is attributed entirely to a single level at 
270 ev, the value ooI?%~11,000 is indicated. 
Thallium has two isotopes of about 30 percent 
and 70 percent abundance to which the levels 
must be ascribed. The two slow neutron induced 
activities in thallium’? have half-lives of 4.23 
min. and 3.5 yr. 
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Fic. 9. The slow neutron transmission of 52.09 Eo 
of thallium. This curve was obtained by using broad 
resolution. The crystal interference effects are prominent. 


* Pool, Cork, and Thornton, Phys. Rev. 52, 239 (1937); 
Heyn, Nature 139, 842 (1937); Fajans and Voigt, Phys. 
Rev. 60, 619 (1941). 
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Columbium 


The slow neutron transmission of Cb was 
studied using a 19.37 g/cm? disk of Cb metal. 
A spectrographic analysis of the sample showed 
less than 1 percent Ta; 0.5 percent Zr; 0.05 
percent Ca, Y, Zn; 0.005 percent Be, Sc, Ti, Ni, 
Rh, Pd, Si; traces of V, Cr, W, Mn, Fe, Ca, 
Au, Al, Ga, Sn, Bi. For all other ‘metals the 
prominent lines were not apparent. The results 
of the measurements using broad resolution over 
a large energy range are shown in Fig. 12. This 
shows the presence of a slight dip near 4 ev, a 
slight 1/v slope above 0.05 ev and pronounced 
crystal interference effects at lower energies. In 
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Fic. 10. The slow neutron transmission of 52.09 g/cm? 
of thallium. This curve represents the region above 0.5 ev 
as measured with better resolution. The iy line may be 
expressed as ¢ = [(9.7+0.2)+(0.6+0.2)E-4]. There is also 
a strong dip at several hundred volts. 
































NOLLD3SS SSOUD 














0.1 





ie) 2 4 6 8 10 
NEUTRON TIME OF FLIGHT MICROSECONDS/METER 


" Fic. 11. The slow neutron transmission of 52.09 g/cm? 
of thallium using the maximum resolution of the spec- 
trometer. The curve shows a strong main level at 270 ev 
and another strong level near 1100 ev. Both dips may be 
owing to several levels. 
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Fic. 12. The slow neutron transmission of 19.37 g/cm? 
of columbium. Pronounced crystal interference effects are 


evident at lower energies. The best 1/v line is ¢ = (6.4+0.2) 
+ (0.10+0.04) E+]. 
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Fic. 13. The slow neutron transmission of 19.37 g/cm? 
of columbium. The transmission of neutrons above 0.5 ev 
is quite constant except the weak dip at 4.1 ev which can 
be attributed to an 0.4 percent Ta impurity. 
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Fic. 14. The slow neutron transmission of 4.70 g/cm? 


of germanium. This represents the measurements as 
obtained by using broad resolution over a wide energy 
range. The presence-of the crystal interference effects 
makes the determination of 1/v line impossible. 


particular there are sharp discontinuities in the 
transmission corresponding to neutron wave- 
lengths of about 2.05 and 2.67A. The best 1/v 
line is ¢o=[ (6.40.2) +(0.10+0.04) E+]. 

The region above 0.5 ev was then studied with 
good resolution as shown in Fig. 13. The trans- 
mission is quite constant except for the weak dip 
at 4.1 ev. This dip has been attributed to a slight 
Ta impurity.‘ From the area of the dip this 
could be accounted for by an 0.4 percent Ta 
impurity which is not unreasonable. Thus no 
levels have been found which can be attributed 
to Cb itself. Cb has only the single isotope Cb®. 
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Fic. 15. The slow neutron transmission of 4.70 g/cm? 
of germanium. This represents the region above 2.2 ev 
as studied with good resolution. 
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Fic. 16. The slow neutron transmission of 4.70 g/cm? 
of germanium. This region is explored using maximum 
resolution. The dip near 95 ev may be multiple in view of 
its peculiar shape. 


The Cb™ nucleus would be of the unstable odd 
Z even A type, and this probably accounts for 
the absence of observable levels.” 5 Nevertheless, 
an induced activity in Cb by an (mn, y) reaction 
has been identified." It has a half-life of 6.6 min. 


Germanium 


The slow neutron transmission of Ge was 
studied using a 4.70 g/cm? disk of Ge metal. 
A spectrographic analysis of the sample showed 
less than 0.1 percent of Fe; 0.05 percent Ca, Cr, 
Si; 0.005 percent of W; traces of Mg, Mo, Ma, 
and Cu. For all other metals the prominent lines 
were not apparent. The results of the measure- 
ments using broad resolution over a wide energy 
range are shown in Fig. 14. Crystal interference 


1 Pool, Cork, and Thornton, Phys. Rev. 52, 239 (1937); 
Sagane, Kojima, Miyamoto, and Ikawa, Proc. Phys. 
Math. Soc. Japan 22, 174 (1940). 
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effects in the thermal region makes the determi- 
nation of the 1/v slope impractical. The region 
above 2.2 ev was next studied using good resolu- 
tion as shown in Fig. 15. This shows a small dip 
near 90 ev. The region above 50 ev was then 
studied using maximum resolution as shown in 
Fig. 16. This shows a dip near 95 ev which may 
be multiple in view of the shape of the dip. If it 
is because of a single level at 95 ev then oI 
~800. The rise in transmission above 10,000 ev 
probably indicates a decrease in the scattering 
cross section at higher energies. There are five 
stable isotopes in germanium. The three well 
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established activities induced in germanium by 
(mn, y) reactions are 40 hr. Ge™, 89 min. Ge’®, 
and 12 hr. Ge7’. 

We take pleasure in acknowledging our in- 
debtedness to Professor J. R. Dunning for his 
valuable suggestions and stimulating discussions. 
We wish to thank Miss Miriam Levin, who 
assisted with the numerous calculations involved 
in this paper, and the other members of the 
cyclotron staff who aided with these measure- 
ments. 


2” Sagane, Phys. Rev. 55, 31 (1939); Sagane, Miyamoto, 
and Ikawa, Phys. Rev. 59, 904 (1941). 





PHYSICAL REVIEW 


VOLUME 71, 


NUMBER 3 FEBRUARY 1, 1947 


The Electrical Charge on Precipitation at Various Altitudes and Its 
Relation to Thunderstorms 
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The free electrical charges on individual precipitation 
particles were measured at various altitudes up to 26,000 
feet by an induction method that avoided touching them. 
In a weak cold front exhibiting no thunderstorm activity, 
and having the freezing level at 11,000 feet, positive 
charges averaging 0.033 e.s.u. were observed from 10,000 
to 26,000 feet. Negative charges averaging 0.04 e.s.u. were 
measured from the surface up to 20,000 feet. Positive 
charges were not observed below 10,000 feet and negative 
ones were not detected above 20,000 feet. Between these 
levels a mixture existed. The free charge on a large number 
of the particles was so great that the electric field at their 


LTHOUGH it has been two centuries since 
Benjamin Franklin flew his kite to demon- 
strate the electrical character of lightning, little 
real progress has been made in the fundamental 
understanding of the electrical charge generating 
and transferring processes in thunderclouds. 
Because every lightning discharge in a cloud 
profoundly changes the electrical state within it, 
surveys of the cloud capable of detailed inter- 
pretation must be completed in a time short 
compared to the charging time. This means that 
an active thunderhead must be explored in 
something less than 10 seconds or, if a longer 
time is necessary, a steady state cloud must be 


surface was an appreciable fraction of the breakdown field 
for air, showing that powerful electrifying agencies exist 
even under rather quiet frontal conditions. Electric field 
measurements on the airplane showed that the particle 
charges were largely neutralized by nearby charges, pre- 
sumably on the air or on the tiny cloud droplets. It is 
shown that the removal of this neutralizing charge will 
immediately produce thunderstorm electric fields and 
potentials. A time plot of the electric field on the surface of 
a flying airplane during the interval of a lightning strike 
is given. 


chosen in which the convective activity is so low 
that disruptive electrical discharge does not take 
place. The Army-Navy Atmospheric Electricity 
Research Project has had assigned to it a Flying 
Fortress that has been instrumented to investi- 
gate the electrical state of the lower atmosphere. 
This airplane has flown through almost every 
conceivable type of weather during the past two 
years to determine the character and magnitude 
of the electrical effects associated with cloud and 
storm formations.' 


‘Gunn, Hall, and Kinzer, “The Precipitation-Static 
Interference Problem and Methods for its Investigation,” 
Proc. Inst. Rad. Eng. 34, 156P April (1946). 
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Fic. 1. (a) Average electrostatic charge on precipitation 
particles at various altitudes: (b) Corrected temperature 
as measured with M1313 aircraft psychrometer. 


The present paper is concerned with the elec- 
trical characteristics of precipitation falling 
within a weak cold front. The observed changes 
in altitude ‘and speed of the airplane flying 
through and at right angles to this front showed 
that there was negligible convective activity. 
Further, the electrical fields resulting from such 
activity were measured and found to be small. 
The observations showed, too, that measurable 
electrical charges were carried by the precipita- 
tion and that these charges, through appropriate 
mechanical separation, were sufficiently large to 
produce high electrical potentials and lightning. 
Measurements taken in active thunderstorm 
areas show that the charge per droplet (or flake) 
may be as much as an order of magnitude greater 
than those encountered in a mild cold front. 
Because of the stability of the’front employed in 
this investigation, measurements made therein 
are considered more significant than would be 
the case if lightning discharges were present 
that suddenly might change the whole electrical 
configuration of the storm area. 

The weather of the present investigation was 
a weak occluded front of the cold type accom- 
panied with a low pressure trough oriented 
North-South, that moved eastward across Min- 
nesota on July 27, 1945. As the front progressed, 
it became more of a cold front in character. 
Steady rain appeared over northern Minnesota. 
The project airplane flew through the front again 
and again at various altitudes and at right angles 
to the direction of the front collecting pertinent 
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meteorological and electrical data. In addition 
to the usual ‘flight instruments and 
meteorological thermometers, the airplane Was 
equipped with two important specialized instry. 
ments designed to measure the electric field at 
the surface of the airplane and the electrica) 
charge resident on the droplets of the encoyp. 
tered precipitation. 

The electric field meters** are induction de. 
vices consisting of a highly insulated electrode 
exposed outside the airplane to the electric fielg 
and connected to the grid of a vacuum tube. 
The exposed electrode is arranged to be alter. 
nately exposed and shielded by a grounded 
metallic rotor driven by a small motor. In the 


presence of an electric field, an alternating volt. . 


age is produced on the grid of the vacuum tube 
that is proportional to the field. The output from 
the first vacuum tube is passed through an 
amplifier and thence to a direct current meter 
by way of a mechanical rectifier driven in syn- 
chronism with the grounded metallic rotor. Since 
such a system is phase sensitive, the indicating 
meter reads not only the magnitude of the 
electric field but also its direction. The electric 
field meters are mounted both on the top and on 
the bottom of the fuselage in such a way that it 
is possible to determine the probable origin of 
any measured electric field. These devices have 
been invaluable in the assessment of data and 
in determining the electrical state of the atmos- 
phere.* 

Some years ago the pioneer results of 
Gschwend,*' who simultaneously measured the 
electrical charge and size of individual raindrops, 
were reinvestigated by a new method. Two 
highly insulated inducing rings separated verti- 
cally by 75 centimeters were mounted below a 
suitable collimator and connected to a high gain 
amplifier, the output of which was connected to 
a recording oscillograph. Raindrops falling 
through the two rings induced two successive 
deflections of the oscillograph, the amplitude of 
which determined the sign and magnitude of the 


?R. Gunn, Phys. Rev. 40, 307 (1932). 

3 Waddel, Drutowski, and Blatt, Proc. Inst. Rad. Eng. 
34, 161P April Loe 

‘ Stimmel, ers, Waterfall, and Gunn, Proc. Inst. 
Rad. Eng. 34, ne P April (1946). 

5 Gschwend, Jahrb. Radioakt. u. Electronik. 17, 62, 19 
(1921). 
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charge on the raindrop, and the time interval 
between the two deflections measured the time 
it took the drop to fall the distance between the 
two inducing rings. Thus, the charge magnitude, 
sign, and the velocity of fall (or droplet size) 
could be simultaneously determined. This meth- 
od has now been adapted for use on aircraft 
with considerable success. 

A truncated sheet-metal cone is mounted on 
the airplane with the smallest circular opening 
of 9.62 cm? area presented face forward in such 
a way that precipitation coming through this 
opening will not touch the inner sides of the 
cone. Thence the particles move outward from 
the big end of the cone without being broken up 
or disturbed in any way. Inside the cone are 
mounted two inducing rings; the forward one is 
5 centimeters in diameter and 5 centimeters long, 
while the aft ring is 12 centimeters in diameter 
and 5 centimeters long. The forward ring is 
highly insulated and connected through appro- 
priate wires and shields to a pressurized cathode 
ray oscillograph equipped with a slow speed 
motion picture camera. The aft ring in the wider 
part of the cone acts as an electrostatic shield 
and is grounded to the cone and to the frame of 
the aircraft. The whole assembly is mounted on 
a streamlined strut 45 centimeters below the 
fuselage where it is moderately well shielded 
electrically from the moving propellers and from 
sharp points that might break into corona. It is 
found with such a device that very few precipi- 
tation particles ever hit the inducing ring and 
that it is easy to maintain the necessary ‘high 
insulation even in the heaviest rain. Raindrops 
striking the leading edge of the cone and breaking 
up produce charges in which the figures of the 
pulses impressed on the oscillograph are notably 
different from a charged droplet proceeding 
through the inducing ring in a normal manner. 
Thus, it is possible to discard all traces except 
those caused by droplets that have traversed 
the apparatus without touching it. The output 
from the inducing ring is fed into a standard 
high sensitivity amplifier and then to a pressur- 
ized cathode ray oscillograph. A 35-millimete; 
movie camera is arranged to take pictures of 
each single sweep. The oscillograph in general is 
calibrated by employing standard a.c. voltages 
whose amplitudes are periodically photographed 








on a calibration frame. The calibrations for 
electrical quantity and for the capacity of the 
complete system were determined by using a 


‘small air gun and putting charges of various 


magnitudes on the shot fired therefrom. These 
charges on the shot were measured after tra- 
versing the inducing ring by being captured in a 
Faraday cage connected to an electrometer. 
Supplementing this information with a knowl- 
edge of the voltage sensitivity corresponding to 
any given setting of the amplifiers and the 
capacitance (which was 357 centimeters), the 
quantity sensitivity is easily checked in flight. 
The time of sweep per traverse of the oscillograph 
was usually 0.169 second. The scale of the 
measuring equipment and sweep rate were se- 
lected so that each sweep usually contained 
from 3 to 10 impulses whose magnitude and 
direction determined the charge on each indi- 
vidual droplet traversing the inducing ring. 
Coherent sets of data were obtained using the 
above apparatus for the mild cold front described 
in an earlier paragraph and at various altitudes. 
It may be seen from Fig. la that electrical 
charges reside on the precipitation at all altitudes 
from the highest observed at 26,000 feet on down 
to the ground. In this investigated front, the 
charge on the precipitation at high altitudes is 
positive in character while at low altitudes the 
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Fic. 2. Fraction of particles having observed 
electrostatic charge. 
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Fic. 3. Vertical distribution of -the number of charged 
particles per unit volume. 


charge on the particles was negative. However, 
at intermediate altitudes, for example from 
20,000 feet on down to 10,000 feet, a mixture of 
positive and negative droplets was encountered. 
At 16,000 feet no positive droplets were meas- 
ured, but at 14,000 feet a few and at 12,000 feet 
quite a number were observed. However, at 
lower altitudes, only negatively charged particles 
were encountered. The average electrical charge 
on some 500 positive droplets above 10,000 feet 
altitude approximated 0.033 electrostatic unit. 
The average negative charge on about 700 
particles at altitudes below 20,000 feet was 
somewhat greater and approximated 0.04 electro- 
static unit. Figure 1a shows plainly the distribu- 
tion of charge on precipitation at various alti- 
tudes both for the positive and negative particles. 
It is thought to be of considerable theoretical 
importance that the freezing level is at approxi- 
mately 11,000 feet and is below the general area 
of maximum electrification of the particles. The 
corrected temperatures as determined from a 
precision ML313 aircraft psychrometer mounted 
outside the fuselage are plotted in Fig. 1b. The 
accuracy of such temperature measurements, in 
common with all such measurements on aircraft, 
is subject to some uncertainty. 

The number of droplets measured at any one 
altitude was hardly sufficient to permit a sta- 
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tistical determination of their free electrical 
charge, but the average over the whole 

is probably significant. In Fig. 2 the measured 
electrostatic charge on the intercepted d 

is plotted against the fraction of the droplets in 
any charge range. It may be seen that a 
number of the droplets (or snow flakes) 
charges less than 0.02 e.s.u. and that the number 
carrying charges larger than 0.1 e.s.u. is relatively 
small. The difference in distribution between the 
positive and negative particles is not thought to 
be significant. 

The reason why the charges on the particles 
seldom exceed 0.1 e.s.u. is easily understood, 
If a particle, for example, is assumed to be 
typical of ‘‘medium rain,” then its radius* ap. 
proximates 5X10~ centimeter, and the elec. 
tric field at its surface for a total charge of 0,1 
e.s.u. is at least as great as 12,000 volts/centj- 
meter. This value approximates the dielectric 
strength of air. It seems probable, therefore, 
that the particle electrification processes even 
in a quiet frontal rain area are so potent that 
typical droplets carry an appreciable fraction of 
the largest charge possible in.an ordinary atmos. 
phere. Since the dielectric strength of air js 
proportional to the pressure, the charge limita- 
tion imposed by corona discharge plays a rela- 
tively larger part at high altitudes than it does 
at the earth’s surface. It is worth re-emphasizing 
that these charges were measured in a weak 
cold front that showed little evidence of im- 
portant thunderstorm or electrical activity. 

It has also been found possible to determine 
the density or the number of appreciably charged 
droplets per unit volume while flying through 
normal precipitation. It should be clear that 
with a known cross-sectional area exposed ahead 
of the inducing ring and with known values for 
the corrected air speed of the airplane, it is 
possible to determine the volume of air swept 
out per unit time by the opening in the cone. 
With values of this volume and the length of 
time that a given number of impulses are meas 
ured or observed on the oscillograph, it is possible 
to determine the number of charged drops swept 
up per unit volume. It is recognized that this 
procedure will not measure the number of tiny 


®R. Gunn, “Electricity of Rain and Thunderstorm, 
Terr. Mag. Atmos. Elec. 40, 79-106 (1935). 
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droplets which bear very small charges, but 
these are of small importance in transporting 
water out of the cloud. In reducing the data, 
the number of impulses measured on the oscillo- 
graph was determined for every frame taken and 
the time of sweep for each frame noted. In Fig. 3 
is given a plot of the number of charged precipi- 
tation particles per cubic centimeter as a function 
of the altitude above the ground. It should be 
emphasized that the number of drops per cubic 
centimeter is the number of drops with sufficient 
charge on them to produce a noticeable deflection 
on the oscillograph, or greater than about 0.01 
e.s.u. Some drops were undoubtedly overlooked. 
It is interesting to notice that at 26,000 feet the 
number of drops per cubic meter is only about 
130, and the number increases more or less 
irregularly as the altitude decreases to about 
12,000 feet. At this level, which is near the 
freezing level, the number of particles per cubic 
meter approximates 750 and is a maximum. 
Below the freezing level, the density decreases 
rather rapidly to much smaller values. This 
seems to be of interest in interpreting the possible 
modes of formation of rain because if every 
raindrop in the lower atmosphere is assumed to 
come from a melting snow, or water and ice 
particle, as Bergeron has suggested,’ it seems 
obvious that there should be more particles 
per unit volume at the freezing level than there 
would be raindrops below. This is because the 
raindrops will usually fall at a greater speed than 
the snow, or the water and ice particles, and 
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Fic. 4. Vertical distribution of particle space 
charge density. 


™“On the Physics of Cloud and Precipitation,” Meteor. 
(Union Geodesique et Geophysique Internationale, Lisbon, 
1933) pp. 156-170. 
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Fic. 5. Measured electric fields on B-25 before and after 
lightning strike, August 15, 1944, 13,000 feet. 


therefore there should be a maximum in droplet 
density at the melting level. Our meager data 
suggest that this maximum may be real. A series 
of similar measurements in an active thunder- 
storm gave essentially similar results. 

It should be evident from the foregoing dis- 
cussion that with a knowledge of the charge on 
the precipitation particles and the number of 
particles per unit volume, one may calculate the 
particle space charge. Values of the particle space 
charge are plotted in Fig. 4 both for the positive 
and negative particles (solid line) and for the 
net particle space charge (broken line). Particle 
space charge must be clearly differentiated from 
the true space charge because neutralizing 
charges always surround rain droplets, and these 
must be counted in evaluating the true space 
charge. That neutralizing charges exist some- 
where outside of each droplet in the cold front 
was well demonstrated by the fact that the 
measured electric fields throughout the area 
were, in every case, orders of magnitude smaller 
than one would calculate from the free charges 
actually resident on the precipitation. In this 
connection it is well to remind the reader that if 
the region of the free charge distribution is large 
compared to the vertical height above the 
ground, then the vertical electric field E is given 
nearly enough by 


E=45 f pde-+Ey (1) 


where p is the true space charge density; Ep is 
the originally existing field produced by surface 








distributions of free charge; and z is the vertical 
coordinate. Since E» is usually negligible, E/4r 
as measured on the airplane gives the value of 
the free charge integrated throughout a prism 
of unit area and extending vertically through 
the front. 

The particle space charge is considered to 
be of great importance because, if relative mo- 
tion is brought about by any means whatever 
between the particles themselves and the air or 
cloud droplets in which they are immersed, and 
if there is a separation of the two, then free 
charge will appear in the adjacent area. In the 
presence of adequate air currents, this surplus 
charge may obviously be some appreciable frac- 
tion of the space charge appearing on the parti- 
cles alone. If the free space charge distribution 
is known, then the electric field and the potential 
distribution may be immediately calculated by 
standard electrostatic procedures. If it is assumed 
that the free space charge is equal to the space 
charge on the precipitation and that no compen- 
sating charge exists on the air outside the drop- 
lets, it is easy to calculate, using the observed 
particle space charge density shown in Fig. 4, 
that in a cloud area in which vertical convection 
takes place, that electric fields up to 20,000 volts 
per centimeter and potentials up to 10° volts 
may immediately be produced. A calculation of 
this nature has been given in an earlier paper.*® 
These calculated electric fields and potentials are 
an order of magnitude greater than those re- 
quired for a lightning discharge and thus show 
that a partial separation of the charge is adequate 
to produce observed thunderstorm electrical 
effects. For example, in our B-25 experimental 
airplane, on August 15, 1944, we flew through 
an active thundercloud and were struck by 
lightning. An electric field indicator mounted on 
the belly of the plane’ near the trailing edge of 
the main wing provided values of the electric 
field both before and after the strike even though 


some of our other measuring equipment was’ 


destroyed by the bolt. The strike occurred at 
about 13,000 feet. The electric field increased as 
shown in Fig. 5 to a value of about 3400 volts 
per centimeter in about ten seconds and re- 
covered in three seconds. These measured values 
of electric field must be corrected for distortion 
due to the presence of the airplane in the thunder 
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cloud. Other measurements show that if thie 
electric field at the time of the strike was ap. 
proximately vertical, then the correction for 
distortion was such that the electric field exist 

in areas away from the plane and at 13,000 feet, 
closely approximated 1600 volts per centimeter 
just prior to the stroke. Thus, making a reason. 
able estimate of the dimensions of the Cloud, one 
is led to the conclusion that the measured values 
of electric field and potential are smaller by 
perhaps an order of magnitude than the values 
which would be produced if all the charge op 
the precipitation and that outside of it were 
suddenly separated by a distance equal to the 
dimensions of a typical cloud. 

The significance of these results in relation to 
thunderstorm electricity is obvious. The meas. 
ures show that in a mild cold front, there is a 
definite tendency for positive electricity to reside 
on the precipitation at high altitudes; whereas, 
at relatively low altitudes, there is a marked 
tendency for a negative charge to exist on the 
precipitation. Further, simultaneous measures of 
the electric field show that the integrated true 
space charge at all levels is small and neutralizing 
charges outside the precipitation usually exist, 
Clearly therefore, if vertical air currents separate 
the associated neutralizing positive and negative 
charges to any appreciable extent, large electric 
fields and voltages sufficiently high to cause 
lightning will be produced. That a redistribution 
of electric charge will be produced in the event of 
active vertical air currents can hardly be doubted. 

Our data lead to the conclusion that even 
rain clouds are intrinsically able to produce 
lightning, and only mechanical effects are needed 
in an ordinary rainstorm to separate the already 
existing charge in the cloud from that on the 
associated precipitation. 

The reasons why positive electrical charge 
predominates on precipitation at high altitudes 
and negative at relatively lower altitudes are 
still controversial. We are continuing our meas 
urements on the charges measured in flight on 
precipitation, and it is believed that as further 
data become available, we may have, finally, 
a first-class clue as to the origin of this charging. 

The author is glad to acknowledge the assist- 
ance of Franklin E. Waterfall in securing much 
of the presented data. 
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Nuclear Excitation Functions.* II. Al*’(d; p, «)Na** 





E. T. CLARKE 
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received October 24, 1946) 


Hitherto unknown in light elements, a new reaction involving the emission of a proton and 
an alpha-particle has been observed by bombarding aluminum with 14-Mev deuterons. By an 
analysis of its excitation function, it is shown that the activity is not generated by the well- 
known (, a) reaction, but is the result of the multiple disintegration of the compound nucleus 


(AF"+4). 





INTRODUCTION 


N the course of the investigation of excitation 

functions of reactions induced in magnesium 
by deuterons, two important radioactive isotopes 
were produced and studied : the 3.0-year isotope 
Na™ produced from Mg™, and 14.8-hour Na™ 
from Mg**, both by the (d, a) reaction on these 
stable isotopes of magnesium. The excitation 
functions were found by use of the stacked foil 
technique! in which a series of aluminum foils 
coated with evaporated magnesium was bom- 
barded with 14-Mev deuterons. The reaction 
Mg*(d, a)Na*® yielded the expected curve,’ but 
the curve of Na™ vs. deuteron energy was 
anomalous in that it showed a sudden rise at high 
bombarding energies. Further investigation of 
this excitation function showed that the in- 
creased activity at the high energy end was 
being generated in the aluminum foils. 

The only manner in which Na*™ has been 
reported produced from Al?’ is through the (n, a) 
reaction, and this would not be expected to show 
a sharp decrease in excitation through a stack of 
foils only } mm in thickness. The new excitation 
curve, therefore, was studied in some detail in an 
effort to ascertain the reaction actually taking 
place. 


EXPERIMENTAL PROCEDURE 


. Anumber of aluminum foils were first stacked 
up and bombarded with the full deuteron beam 
of the M.I.T. cyclotron. Half-life measurements 
and a chemical analysis of the activity induced 


* The research described in this article was supported in 
part by Contract NSori-78, U. S. Navy Department. 
(1948) T. Clarke and J. W. Irvine, Jr., Phys. Rev. 66, 231 
: a T. Clarke and J. W. Irvine, Jr., Phys. Rev. 69, 680 


187 


in the first few foils were carried out, the results 
of which confirmed the assumption that the 
isotope being produced from the aluminum was 
14.8-hour Na*™. A stack of fifty foils, each 4.55 
mg/cm? thick, was then prepared from sections 
of aluminum sheet of known area and weight. 
Variations in thickness of the foils were found by 
this means to be less than 1 percent. After the 
deuteron beam energy had been determined,' the 
stack of foils was given a deuteron exposure of 
0.81 microampere hour and allowed to decay for 
two days to remove the shorter lived components. 

Activities in the fifty foils were measured on a 
bell-type mica-window beta-ray counter*® and a 
scale-of-128 circuit. An aluminum absorber of 
80 mg/cm? was used to bring the counting rate 
to a few thousand counts per minute, and each 
foil was placed on the absorber to minimize 
changes in geometry from sample to sample; the 
distance between counter and absorber was in- 
creased to handle the higher intensities in the 
first few foils. Readings were corrected for back- 
ground, about 5 percent of the lowest observa- 
tion, and for decay of the Na™ during measure- 
ment. Finally the stack of foils as a whole was 
compared with a 10-microgram radium standard 
on a calibrated platinum screen wall gamma-ray 
counter,‘ and was found to contain® 3.30 rd of 
Na at the time of bombardment; this quantity 
was used to determine the absolute cross section 
of the reaction. 

The activities induced in the foils were then 
plotted as a function of deuteron bombarding 


a a Good, A. Kip, and S. Brown, Rev. Sci. Inst. 17, 262 
‘A. Roberts, L. G. Elliott, J. R. Downing, W. C. 
Peacock, and M. Deutsch, Phys. . 64, 268 (1943). 
5 1 rutherford (rd) = 10° disintegrations per second; E. U 
Condon and L. F. Curtiss, Phys. Rev. 69, 672 (1946). 
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Fic. 1. Experimental excitation function of Na™ obtained by bombardment of aluminum foils. The ordinates, p= 


for three ran 


, are absolute, though in arbitrary units such that 1000= 1.74 x 10’ disintegrations/minute at bom 


ment time. Activation beyond the deuteron range is neutron induced. The Na*(d, p)Na™ activity is caused by the 


presence of 0.02 percent sodium contamination. 


energy, the result of which is presented in Fig. 1. 
The curve can be divided into three distinct 
sections: the initial part, beyond the range of 
the deuterons, showing a constant excitation 
due to neutron activation; probably both by 
AF"(n, a)Na™* and by Na™(n, y)Na™*; a region 
from 1 to 9 Mev in which activation first rises 
and then levels off with increase in bombarding 
energy, ascribable to sodium contamination in 
the aluminum metal since the curve corresponds 
to 0.02 percent Na producing Na*(d, p)Na™; 
and the interval from 9 to 14.5 Mev in which the 
new reaction sets in with a rapidly increasing 
cross section with rise in energy. 


INTERPRETATION OF DATA 


From the appearance “of the experimental 
excitation function it seemed apparent that the 
rapid rise in production of Na™ with deuteron 
energy could not reasonably be ascribed to neu- 
tron (, a) activation, since practically all of the 
product was generated in the first 0.3 mm of the 
aluminum target. The constant amount of Na™ 
generated in each foil beyond the range of the 
deuterons further bore out the hypothesis that 





targets of thicknesses of the order of deuteron 
ranges are thin for neutrons, and that no ap- 
preciable variation in neutron excitation, at least 
in light elements, can be expected in them. 

It remained, therefore, to consider the possible 
deuteron reactions which might be responsible 
for the observed results. In disintegrating from 
13Al?7+,H*—,,Na™, 5 units of mass and 3 of 
charge must be emitted. Table I shows the four 
possible ways, together with the calculated 
threshold energies, in which the transition can 
result in emission of known elementary particles. 
It is evident that the only possible reaction is 
that of emission of a proton and an alpha-particle, 
all other modes requiring considerably more 
energy than can be imparted by 14-Mev deu- 
terons. 

Assuming that the reaction Al?’(d; p, a)Na™ 
takes place through the formation and disinte- 
gration of a compound nucleus (the Oppen- 
heimer-Phillips mechanism has no meaning at 
deuteron energies high compared to that of the 
barrier, since in this region the deuteron break-up 
occurs inside the target nucleus), the shape of 
the excitation curve should, by analogy with 
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(n, p) and (m, a) excitation functions, be gov- 
erned by the penetrabilities of the proton and 
the alpha-particle outward through the nuclear 
barrier. In the case of proton emission followed 
by loss of an alpha-particle, the total average 
penetrability can be expressed : 


Ei 


1 
Py, a(E:) ~ E. P, (Ey) Pa(Ea)dEy, (1) 


i 


where E; represents the total available energy 
for the two particles, E, and E, the energies 
(including recoil), and P, and P, the individual 
penetrabilities of the proton and of the alpha- 
particle, respectively. The energies are related 
by the equation E;=E,+£.,. A similar integral 
can be formed for the case of alpha-emission 
followed by proton loss; it will be slightly 
different because of differences in the masses and 
charges of the emitting nuclei. ' 

In order to predict the absolute cross section 
for the (d; p, a) reaction, it is assumed that in 
the process of decay of the compound nucleus, 
any particle or set of particles energetically 
possible will emerge, and that the relative pro- 
portions of the various products are determined 
mainly by the probability of formation k, of the 
emergent particles in the nucleus and by their 
penetrability P, through the nuclear barrier. 
The cross section og,z, then, may be written 


O4,2= Qk2P2/Zk2Pz, 


where Q, is the probability of formation of a 
compound nucleus, and the summation extends 
over all possible products, including multiple 
disintegrations. In the assumed decay mechan- 
ism, the factor kz=1 for elementary particles, 
less than unity for combinations such as alpha- 


particles. 
Accordingly, 


O4,a=QeRaPa/Dk:Ps (2) 
and | 
G4; p.a=QR kal’ y,a/UkePs; 
G4; a,p=QckakpPa, p/TkzPs- 


The observed cross section should equal the sum 
Of ¢4:p,a+¢4:0,p- Adding Eqs. (3) and substi- 
tuting Q.ka/DkzP.=0c4,2/Pa obtained by rear- 


(3) 


NUCLEAR EXCITATION FUNCTIONS 





‘rangement of Eq. (2), we find 


Paths Pay G4,a(Py,at+Pa,») (4) 








Fobs = ky 


since P, = 1 because the energy region in question 
is considerably above the barrier energy for the 
(d, a) reaction (6.5 Mev). 

The deuteron energy Eg is corrected for recoil 
of the (Al?7+d) compound nucleus and related 
to the total available energy E, by the formula 


E oe 5 
— M i-Q), ( ) 





where M is the mass of Al*’, m the deuteron mass, 
and Q the reaction energy of the reaction 
AP’"(d; p, a)Na™. 

In the determination of the quantities P,, .(E;) 
and P., »(Z;), the only arbitrary factor is that of 
the nuclear radius. The value used here was that 
reported both by Grahame and Seaborg* and by 
Sherr,’ who found R=5.43X<10-" cm for alumi- 
num by measuring both elastic and inelastic 
scattering of fast neutrons from various targets. 
This radius was divided by A'=3.00 to obtain 
the standard nuclear radius ro>=1.8110-" cm 
which was then applied to the penetrability 
formulae given by Bethe.* These penetrabilities 
were inserted in Eq. (1) and integrated graphi- 
cally. 

Since the (d, a) cross section for aluminum has 
not as yet been measured, the previously deter- 
mined excitation function? for the reaction 
Mg?*(d, a)Na™ was used as an approximation. 
The general trend of these reactions at high 
energies is a decrease inversely proportional to 
the deuteron energy, and therefore the cross 
section can be expressed as a function of Ez 


TABLE I. Threshold energies for possible reactions. 








Particles emitted Deuteron threshold 





(Al? +d) —Na™ energy, Mev 
at+p 5.3 
He®+d 25.2 
2Wd+p 30.8 
3p+2n 35.5 








*D. C. Grahame and G. T. Seaborg, Phys. Rev. 53, 
795 (1938). 
7R. Sherr, Phys. Rev. 68, 240 (1945). 
*H. A. Bethe, Rev. Mod. Phys. 9, 161 (1937). 
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Fic. 2. Comparison of predicted excitation of Al*"(d; p, a)Na™ with the observed data. The calculated absolute 
values are within a factor of two of those found experimentally. Two ordinate scales are used to show the correspondence 


at lower energies. 


(in Mev): 
Ga,a= (7.0/Ea) X 10-*5 cm?. (6) 


Using Mattauch’s® values for the neutral 
masses of Al*”(26.98980) and Na?4(23.99779), the 
reaction energy Q was found to be —4.90+0.30 
Mev. However, since only the difference in 
masses between Al?’ and Na*‘ was needed, an in- 
dependent calculation for Q was carried out, 
using the reaction’ Mg”4+,He*—Al?’ +H! — 1.82 
Mev and the recently determined" beta-ray dis- 
integration Na**-+Mg**+5.51 Mev; this yielded 
Q=-—5.10 Mev. ‘An average value of Q= —5.0 
Mev was used here since Mattauch had employed 
‘ an erroneous figure (5.3 Mev) for the Na* dis- 
integration energy in checking his mass deter- 
minations. 


(1943) Fliigge and J. Mattauch, Physik. Zeits. 44, 181 
10 Observed b ee ee and Miller, Proc. Roy. Soc. 
146, 396 (1934};_ and values recalculated by M. S. 


Livingston and H the, Rev. Mod. Phys. 9, 300 
(1937). 


uL.G. Elliott, M. Deutsch, and A. Roberts, Phys. Rev. 
63, 386 (1943). 


The results of these calculations, compared 
with the experimental cross sections, are given in 
Fig. 2, where the predicted curve, multiplied by 
1.48, is seen to fit the observed points. The effects 
of different assumed values for the quantities 1» 
and Q were tested ; it was found that these could 
be changed only by very small factors and still 
retain the shape of the experimental curve. With 
ro held constant, Q is determined to 0.1 Mev by 
this criterion, while r>9 could be varied only by 
about 0.1X10-" cm with Q held within the 
limits of its probable error (0.3 Mev). The 
absolute cross section can be changed by about 
a factor of two by adjusting the variables within 
the limits given above, showing that the observed 
agreement is well within the probable error. 

From the foregoing analysis there can be little 
doubt that the observed reaction is indeed 
Al?"(d; p, a)Na™, and that it proceeds by means 
of a compound nuclear mechanism. The existence 
of this reaction, and the success of the assump- 
tion regarding the competition for decay of the 
compound nucleus, also suggest an explanation 
for the decrease above the nuclear barrier with 
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increasing energy, observed in experimental deu- 
teron-induced excitation curves: competition 
from reactions involving the emission of a second 
particle. The familiar (d, 2”) is a special case of 
such reactions, in that the outgoing particles are 
not limited by the necessity of penetrating the 
Coulomb barrier and can therefore be produced 
with deuterons of moderate energies. The (d; p, m) 
and (d; n, p) reactions which must compete with 
the (d, p) and (d, m) would be difficult to observe, 
since they lead back to the target isotope and 
could be found by these techniques only where 
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excited states of these isotopes exist. For light 
elements in general, where barriers are low for 
deuterons, the (d;m, a) and (d;a,m) reactions 
usually lead to stable isotopes. With the higher 
deuteron energies now becoming available, many 
more of these multiple disintegrations should be 
found. 

The author wishes to thank Professor John 
W. Irvine, Jr., who kindly carried out the chem- 
ical analysis, and Professors V. F. Weisskopf, 
M. Deutsch, and R. D. Evans, with whom he 
discussed the theoretical aspects of the problem. 
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Thermal Ionization of Impurity Levels in Semi-Conductorst 


B. GoopMAN,* A. W. Lawson,** AnD L. I. ScuirF 
Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 
(Received November 6, 1946) 


The probability of ionization of an impurity level by thermal agitation in a semi-conductor 
is estimated by the use of a simple Debye model. It is found that the probability of ionization 
by one phonon decreases exponentially with the depth of the impurity level below the con- 
duction band, approaching zero as the energy separation becomes equal to k@p. The probability 
of ionization by more than one quantum is also estimated by the use of an Einstein model. It 
appears that the probabilities so calculated may play an important role in determining the 
frequency dependence of the rectification efficiency in crystal rectifiers. 


INTRODUCTION 


HE electrical properties of many semi- 

conductors are determined primarily by 
the presence of impurities that introduce discrete 
electronic energy levels which lie between the 
empty conduction band of the semi-conductor 
and the next lowest lying filled electronic band. 
These discrete energy levels may act either as 
acceptor levels or as donator levels. In the first 
case, they lie near the top of the filled band, are 
empty at zero temperature, and can accept 
electrons from the filled band, thus producing 
“hole conduction.” In the second case, the levels 
lie near the bottom of the empty band, are full 
at zero temperature, and can donate electrons 


t This work was carried out on contract OEMsr-388 
between the National Defense Research Committee and 


the University of Pennsylvania and was originally de- 
scribed in a monthly report submitted July 7, 1943. 
* Now Jewett Fellow, University of California, Berkeley. 
' Now at the Institute for the Study of Metals, Uni- 
versity of Chicago. 


to the empty band, thus augmenting the intrinsic 
electronic conduction. The two situations are 
represented schematically in Fig. 1 and Fig. 2, 
respectively. In the following discussion we shall 
consider only donator levels, which produce the 
so-called N-type semi-conductors, with the 
understanding that our conclusions are also ap- 
plicable to acceptor levels forming P-type semi- 
conductors. 

When a semi-conductor is placed in contact 
with a metal, a rectifying junction may be 
formed if the relative work functions of metal 
and semi-conductor are suitable. In the case of 
an N-type semi-conductor, the effective work 
function of the metal must be larger than that of 
the semi-conductor. Then the equilibrium situa- 
tion near the junction is as represented schemati- 
cally by the potential diagram of Fig. 3.! Practi- 


1 F. Seitz and S. Pasternack, “The Principles of Crystal 
Rectifiers," NDRC Report No. 01-102, June 10, 1942. 
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Fic. 1. Energy level diagram showing schematically the 
excitation of electrons from a normally filled band to 
compet levels producing hole conduction. Eo oy 
ae of the acceptor levels above the normally filled 
nd. 
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Fic. 2. Energy level diagram showing schematically the 
ionization of donator impurity levels to produce free 
electrons in the conduction band of a semi-conductor. 
hod - the depth of the levels below the conduction 








cally all of the donators lying more than thermal 
energy (kT) above the thermodynamic potential 
far from the junction (which is about midway 
between the donator levels and the bottom of the 
conduction band) will be ionized. This produces a 
net positive charge in the semi-conductor near 
the metal, and consequently a net negative 
charge in the metal on the other side of the 
potential barrier. Because of this double layer or 
potential step, the electron flow from semi- 
conductor to metal when the metal is made 
anodic by a given potential is greater than the 
reverse flow when the metal is made cathodic by 
the same amount. This results in a rectifying 
characteristic of the type shown in Fig. 4. 

In addition to causing a net current through 
the junction, an applied voltage also changes 
the bound charge in the semi-conductor near the 
interface. Thus the amount of stored charge 
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depends on the voltage, and the contact has sorhe 
capacitative properties. The over-all electrical 
behavior of the rectifier may be represented 
approximately by the simple circuit of Fig. 5, 
It is evident that the capacity associated with 
the bound charge in the barrier layer will tend to 
short out the back resistance when the frequency 
becomes high enough that the impedance of ¢ 
is comparable with R,.? For high frequencies, 
this should cause the rectification efficiency to 
fall off as the inverse square of the frequency. 
Experimentally, however, in many rectifiers, the 
efficiency drops off much more slowly than this 
with increasing frequency. One possible explana- 
tion of this discrepancy lies in assuming that 
the circuit parameters of Fig. 5 are themselves 
functions of the frequency. The capacity C will 
be expected to decrease with increasing frequency 
if there is a time lag in ionization’ (and deioniza- 
tion) of the donator levels that is of the order of 
the period of the applied field. Such a lag can 
occur owing to the finite probability of ionization 
by thermal agitation of the lattice; an estimate 
of this probability is made in the following 
section. 


THEORY 
The transition probability per unit time for 


the excitation of electrons from donator levels 


is given by 
w= (2x/h)| Vox|*px, (1) 


where px, the energy density of states in the 
conduction band, is given by 


px=4rQmp/h. (2) 


m and p are the mass and momentum of the 


,, >b+%, 


e) CH) eo. 


Fic. 3. Schematic diagram showing the relative energy 
levels in a metal and N-type semi-conductor. With no 
external applied field, the top of the Fermi distribution 
lies approximately half-way between the impurity 
and the conduction band of the semi-conductor. 














2H. A. Bethe, ‘Theory of High F ye’ Rectification 
by oy Crystals,” NDRC Report No. 43-11, October 
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electron, respectively, and © is the volume of 
the crystal. Vox is the matrix element 


Yaw f voVox*dr (8) 


calculated in terms of normalized initial and 
final wave functions 


vo= (1/(xa*) ie", 


(4) 
Wx = (1/23) exp(iK-r). 


Yo is taken for simplicity to be of hydrogenic 


form. ’ 
The perturbation V may be approximated as 


V(r) = —R.(r)-VV0(r), (5) 


h 
ee Volt) = —Ze*r/r?, 


(6) 


R, =e, exp(ie-r). 


Equations (5) and (6) may be interpreted as 
follows. The thermal agitation may be resolved 
into acoustic Debye waves which strain the 
lattice and alter the separation of electron and 
donator ion. The total strain will then involve a 
sum over the propagation vectors @ of the Debye 
waves. However, energy must be conserved 
between the phonons (Debye wave quanta) and 
the electrons, so that 


hv, =h?K2/2m+|Eo|, Eo<0, (7) 


where v, is the frequency of the Debye wave and 
E, the energy of the bound electron state. 
Equations (4)-(6) may be substituted into 








; Fic. 4. Typical characteristic curve of a rectifyin 
me poug the dependence of current J on aaliel 
voltage V. 
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RA = 7 











Fic. 5. Simple equivalent circuit of a rectifying contact. 
R, and C, which are both functions of applied voltage, 
represent the resistance and capacitance of the barrier 
layer respectively, R, is the bulk resistance of the semi- 
conductor. 


Eq. (3) and the integrals carried out to give 


T ; =a —— o 
Vox=4iaze*(—~) & cost, (= = ws ), (8) 
3 a*0,? 


where 6, is the angle between e, and Q,=e—k. 
For the situations of practical interest, ¢, K, and 
1/a are all of the same order of magnitude 
(~107 cm). The quantity aQ, will be of the 
order of or less than unity, and a good enough 
approximation to Eq. (8) may be obtained by 
expanding the inverse tangent and keeping the 
first two terms, and replacing aQ, by aK. This 
gives for the square of the matrix element 


167Z*e'a K? 
| Vox | 2 —— cos?6,. (9) 





The amplitude of the Debye wave may be 


written 
hv. 
€,?= (10) 





where M is the mass of the donator atom. Since 
we have to integrate over the Debye spectrum, 
we need the number of modes between », and 


Vetdy, 
S(ve)dve = (3h®v,?/k*®Op*)dv, (11) 
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where 0p is the Debye temperature. Since the 
range of integration is small (|Z)| is close to 
k@p), the result for the total transition probability 
may be estimated as 


1142*Z*e'am?(2mkT)* (Op /T —C) 
Mnh' (e°—1)C? 
| Eo| 
aT. 





Wr 


(12) 





Typical numerical values for substitution in 
Eq. (12) are as follows: @p/T=2, C=3/2, 
Z=0.1, a&10-7 cm, M=27 for aluminum im- 
purity in a silicon semi-conductor. The small 
value of Z includes the effect of the dielectric 
constant (about 10) in reducing the single un- 
balanced charge of the aluminum ion; similarly 
a is larger than would be expected for a hydrogen 
atom by about a factor of 10. These values 
indicate that w is about 10” sec.—', provided that 
| EZo| <k@p.* 

If |EZo| exceeds k@p, which seems unlikely on 
the basis of present information, then ionization 
involves the absorption of more than one phonon 
from the lattice vibrations. The order of magni- 
tude of two or three phonon effects may be 
estimated by using the simpler Einstein model 
for the oscillations of the crystal lattice as 
follows. We assume that the atoms undergo a 
harmonic displacement from fixed lattice points 
given by 

R=R, sinwt, (13) 


which introduces a perturbation of the electron 


* We are indebted to Professor F. Seitz for informing us 
of the results of an independent estimate of the transition 
probability that is based on the Born-Oppenheimer 
molecular approximation. He obtained a numerical result 
of the same order of magnitude. 
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energy given by 


Ze? Ze? Ze’ ao 
a x (R/r)"P,(cosé’), (14) 


where @ is the angle between r and R, and 
|R| < |r|. Using the fact that (Ka)*<1, one finds 
for the matrix element 

4na*Ze? « 2"n! 

(wa®Q)* n—1 


(2n+1)! 
X P,(cos6’) (sinw#)*, (15) 


Vox = 





If we expand Vox in a Fourier series, the first 
harmonic will be the matrix element for transi- 
tions from levels lying between 0 and 6g, the 
second will correspond to transitions from levels 
lying between k@g and 2k6z, i.e., two phonon 
transitions, etc. Using the fact that RK, 
one finds 

, nwa*Ze? 





(tRoK)P;(cos@’), 


Vox= 


3(ma®Q)! 

4nra*Ze? 
coe 
30(xa®Q)* 


etc. 


(16) 


V7 


(tRoK)*P2(cosé’). 


From the foregoing equations, using the numeri- 
cal value suggested above, one may show that 
each additional phonon absorbed reduces the 
transition probability by approximately a factor 
of 10°. 

It may readily be shown that the transition 
probability due to thermal radiation is of the 
order of 10° sec.—!, and hence may be neglected 
in comparison with the probability of single 
phonon transitions. 
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On a General Condition of Heisenberg for the S Matrix 


i S. T. Ma 
Institute for Advanced Study, Princeton, New Jersey 


(Received October 18, 1946) 


It is shown that the S matrix for an attractive exponential potential, which possesses re- 
dundant zeros, does not satisfy a general condition of Heisenberg. To insure the validity of 
Heisenberg’s condition, we introduce the supplementary condition that the interaction poten- 
tial should vanish for large distances from the scattering center. It is shown that the S matrices 
for the attractive exponential and the Coulomb potential cut off at a large distance R give 
correctly both the eigenvalues of energy and the asymptotic behavior of the wave functions 


for the s states in the limit R-. 





1, INTRODUCTION 


ECENTLY Heisenberg! has proposed a new 

formulation of the quantum theory which 
deals exclusively with such physical quantities 
as are directly observable, namely, the discrete 
energy values of stationary states of closed 
systems, and the asymptotic behavior of wave 
functions at infinity in collision processes. These 
physical quantities are determined by Heisen- 
berg’s S matrix, which plays the same funda- 
mental role in Heisenberg’s new scheme as the 
Hamiltonian and wave equation in the current 
scheme of quantum theory. It is possible to 
derive from the current scheme of quantum 
theory certain general properties of the S matrix, 
such as the unitary property’ and the relativistic 
invariance,? which can be taken over into the 
new formulation of the quantum theory. 

The determination of closed stationary states 
from the S matrix has been investigated by 
Heisenberg and Mller. For the following investi- 
gation it is sufficient to state Heisenberg’s result 
for the s states of a particle under the action of a 
central force. Let u,(r) and u(r) be, respectively, 
the radial wave functions of the discrete and the 
continuous energy states, normalized according 
to the equations 


ff sstorrar=t, £: us*(r)uxe(r)dr = 5(k—k’). 


1 W. Heisenberg, Zeits. f. Physik 120, 513 and 673 (1943) 
Heisenberg’s third and fourth papers on the S matrix have 
not yet been published. I am indebted to Professor Pauli 
for communicating to me some results in the fourth paper. 

*W. Pauli, Meson Theory of Nuclear Forces (Interscience 
Publishers, Inc., New York, 1946), p. 49. 

*C. Moller, Danske Vid. Sels. Math.-fys. Medd. 23, 
No. 1 (1945); 22, No. 19 (1946). 
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These wave functions satisfy the completeness 
relation 


Dd a Un*(r)un(r’) +f ux" (r)ux(r’)dk = 5(r—r’). (1) 
0 
The asymptotic expressions of these wave func- 
tions for large distances are 
Un(?) ~Cn(24)-* exp(— | kn Ir), (2) 
ux(r) ~(2/m)* sin[kr+6(k) J, (3) 


where c, is a constant and 4(k) is the phase of 1. 
Making use of the general relations 


S(k) = exp[215(k)] (4) 
S*(k)S(k) =1, (S) 


Heisenberg obtained from (1)-(3) the following 
general relation for a large value of r: 


and 


fsw exp(tkr)dk =>>n|Cn|?exp(—|Ralr). (6) 


The closed discrete states can be determined 
by analytic continuation from the real to the 
complex values of k. If S() decreases sufficiently 
rapidly for large values of |k|, the left-hand side 
of (6) is equal to 277 times the sum of the residues 
of S(k) exp(tkr) at the poles of S(k) in the upper 
half of the complex plane. Identifying the dis- 
crete states with these poles, Heisenberg obtained 


|Cn|*= 2miLResS(k) Jin. (7) 

For the S matrices such that the relation 
S(k)S(—k) =1 (8) 
in the complex plane can be derived by analytic 
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continuation of the unitary condition (5) for real 
values of k, the poles of S(k) in the upper half- 
plane have corresponding zeros in the lower half- 
plane. Hence, when the poles and zeros are on 
the imaginary axis, 

|¢n|?= —2e[(dS/d|k|) -seat (9) 
This relation expressing the absolute magnitude 
of c, in terms of S(k) can also be directly derived 
by a method of Kramers‘ without using (6). We 
shall see below that there are cases in which (9) 
is satisfied but (6) is not satisfied. 

2. REDUNDANT ZEROS 


As recently reported,’ the S matrix for an 
attractive exponential field presents a new fea- 
ture in the eigenvalue problem. For the s states 
of a particle of mass m in the potential field 


V(r) = — Voexp(—r/a), (10) 


where V» and a are positive constants, the wave 
equation takes the form 


— (h?/2m) (d*u/dr?+k*u)+V(r)u=0. (11) 


The solutions of Eq. (11) for a positive value of 
energy are the Bessel functions 


J +i,[.a exp(—r/2a) ], 
p=2ak, a=2a(2mV>)*/h. 


where 


The wave function which vanishes at the origin is 
ux(r) =4(24)-4| P(tp+1)/Jip(a) | 
X { J-ip(a) Ji.[a exp(—r/2a) ] 
—Ji,(a)J_ipLa exp(—r/2a)]}, (12) 


the normalization constant being chosen such 
that the asymptotic expression of Eq. (12) for 
large r, namely, 


ux(r) ~t(2e)-*| Pp +1)/Fip(@) | 








J ~ig(x) (/2)*” 
—ik 
x| Get) exp(—tkr) 
J ip) (a/2)-*” 
- ik 13 
Sap men) 8) 


*H. A. Kramers, Hand- und Jahrbuch d. Chem. Physik 
1, 312 (1938). 
5S. T. Ma, Phys. Rev. 69, 668 (1946). 





is of the form of Eq. (3). Comparing Eq, (12) 
with Eq. (3) and making use of Eq. (4), we see 
that 





Ji(@)P(ip+1) sa\-% 
S(k) = a 
” eee - © 


S(k) vanishes when J;,(a@) vanishes or T( —ip+1) 
is infinitely large. Now J,,(a) has no zeros in the 
lower half-plane of k except those on the imagi- 
nary axis.* Hence the first condition is equivalent 
to 


J\pi(a) =0. (15) 
The wave function corresponding to a zero k, of 
Eq. (15) is’ 
Un(r) =Cn(2r)—*T(| pn| +1)(2/cx)!>m! 
XJlenila exp(—r/2a)]. (16) 
The absolute magnitude of the constant c, js 
given by Eq. (9) with S given by Eq. (14), as 
can be directly verified from the normalization 


condition. 
The second condition 


I'(—ipp+1)= (17) 
is satisfied by values of p such that 
ip=|p|=1, 2, 3---. (18) 


The wave functions Eq. (12) corresponding to 
the eigenvalues of Eq. (18) vanish identically. 
The zeros given by Eq. (18) can be regarded 
as redundant. We shall denote in the following 
the values of p and k corresponding to the 
redundant zeros by p, and k,, respectively. 


3. VALIDITY OF HEISENBERG’S 
GENERAL RELATION 


Let us now study Heisenberg’s general rela- 
tion, Eq. (6), as it applies to the S matrix for the 
exponential potential. The integral 


f “s(k) exp(ikr)dk 


in (6) can be evaluated by contour integration. 


6 Gray, Mathews and MacRobert, Bessel Functions (The 
Macmillan Company, New York, 1931), p. 88. 

7H. Bethe and R. F. Bacher, Rev. Mod. Phys. 8, 111 
(1936). 
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We have 
f ” 5(R) exp(ikr)dk=lim f S(k) exp(ikr)dk 
—o Cc 


+2mni > Res[S(k) exp(ikr)], (19) 


where C denotes a semi-circle of radius K above 
the real axis, having its center at the origin, and 
the summation extends over all the poles in the 
upper half-plane. Equation (14) can be written 


in the form 
S(k) = Mip(a)/A-ip(x) 

if we put 

Jy(a)= (a/2)?A,(a)/T(p+1). 
Now A,(a) approaches the value 1 as p tends to 
infinity through a sequence of numbers not con- 
sisting of negative integers. Hence if K tends to 
infinity through a sequence such that |p| is not 
integral 


lim f S(k) exp(ikr)dk = lim f exp(ikr)dk 
K-00 c Ko | c 


‘ f exp(ikr)dk =2x8(r) =0 


since r*0. Therefore 


S(k) exp(ikr)dk =2mi  Res[.S(k) exp(ikr)] 
= Dal Cn|? exp(— |Ralr) 


aon - (sa). 


—0 


exp(—|k,|7r). (20) 
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It can be easily shown from (14) that 
(dS/d|k|)—ixr>0, (21) 


so that the second sum over the redundant 
eigenvalues is smaller than zero. Comparing Eq. 
(20) with Eq. (6) we see that Heisenberg’s gen- 
eral condition is not satisfied by Eq. (14). 
It has recently been proposed by D. ter Haar® to 
make use of Eq. (21) for discarding the redun- 
dant zeros. 

This example raises a question of the validity 
of Heisenberg’s relation, Eq. (6). Heisenberg’s 
derivation outlined in Section 1 holds good only 
if the difference between the exact wave function 
u,(r) and its asymptotic expression Eq. (3) is 
small in comparison with the wave function 
u,(r). This is not the case with the exponential 
potential, however, as can be seen from Eggs. 
(12), (3), and (2). 

We can see how the difficulty of redundant 
eigenvalues arises from the approximate treat- 
ment by evaluating the integral 


f u,*(r)u,(r’)dk 
0 


for the exact wave function (12). From Eq. (12) 
we have 


un*(r)up(r’)dk = f : f(kir,r’)dk, (22) 


0 —o 


where 


(Ri 1, 1) = (24) (PG p+ 10 (—tp+1)/Fip(a)I—io(a) } {J—ie(a) Jip(a) Jina exp(—1/2a)] 
X J—ieLa exp(—r’/2a) ]—Jip(a)*J_ipLa exp(—1/2a) ]J_i.[a exp(—r’/2a)}}. (23) 


Now the integral of f(k;r, r’) taken round the semicircle C referred to above has the limit 


lim far r)dk= (2x) [ texpLik(r’—r)]—expLib(r’ +1) ah 
c —o 


Hence 


—« 


= 5(r’—r) —8(r' +r) =d5(r’—r). (24) 


f F(R; 1, r’\dk = 5(r' —r)+20i © Resf(k; 1, 1’). (25) 


But at the poles where I'(ip+1) = © the residue vanishes. Hence 
2ni D) Resf(k; 1, 1’) = —i  [ResJ_:,(a)- 0 (ip +1) 1'(—ip+1)Jip(a) 


=— dD Un*(r)un(r’) 


*D. ter Haar, Physica, in the press. 


X J—ipLa exp(—1r/2a) ]J_i[a exp(—r’/2a)] 
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so that \ 
J un*(r)ur(r’)dk = 5(r’ —1) — Din Un*(r)un(r’), 
0 
which is just Eq. (1). In spite of the factor I'(¢p+1) in the integrand of Eq. (22) the SUM Matio, 


in the final result extends only over the true eigenvalues. The above calculation may be taken to 
be a direct verification of the completeness relation. For large distances r and r’, we have, by Eq. (2) 


fw u(rndk =a" 1) — 2x) Daleal? expl = [bal (+79) 09 
0 


In the approximate treatment, however, we have 


f : uz* (r)uz(r’)dk = 8(r’ —r) —(2r)— f : S(k) exp[tk(r+r’) ]dk 
0 —« 


= 8(r’ —r) — (2) Dalen |? expl— | Ral (7 +1’)] 
-1 = Res { S(k) exp[tk(r+r’) }} ian, 


which has an additional term corresponding to the redundant zeros. 

This discussion shows the necessity of introducing a supplementary condition for the validity of 
Heisenberg’s general relation. We can take for the supplementary condition the condition that the 
potential V(r) should vanish for large distances from the scattering center, because under such a 
condition the wave function u(r) is equal exactly to the expression (3) for large distances and 
therefore Eq. (6) holds with certainty. 


4. POTENTIALS CUT OFF AT A LARGE DISTANCE 


Mller has investigated the exponential potential (10) cut off at a large distance R, i.e., the potential 


V(r) =—Voexp(—r/a) (0<r<R) 2N 
=0 ( r>R). 


The wave function u,(7) is now given exactly by (3) when r>R, and, except for a constant factor, 
given by (12) when r<R. The requirement of continuity of u,(r) and its derivative at r=R gives 


J~ip(@) JipsiLa exp(— R/2a)]+Ji,(a) J-+,-1La exp(— R/2a)] 
Jig) Jip-Lae exp(—R/2a)]+-Jip(a) Tips: exp(—R/2a)] 
For large values of || in the upper half-plane of k such that 2a|k| is not an integer, 


S(k) = —exp(—2tkR) 





(28) 


S(k) exp(ikr) ~ —exp[tk(r—2R)] 








JipsiLa exp(—R/ 2a) ] __ Te) (:) 
Jip—Laexp(—R/2a)] (tp +2)\2 


Xexp(— R/a) exp[ik(r—2R)], 
and so 


lim f S(k) exp(ikr)dk 
K-02 c 


vanishes for r>2R.° In accordance with the general discussion of Section 1 the eigenvalues of k 
are the values of k in the lower half-plane satisfying the equation S(k) =0 or by Eq. (28) 


J ~ip(@) J ips La exp(— R/2a)]+Jip(a)J-ip—iLa exp(—R/2a) ]=0, (29) 
® Whittaker and Watson, Modern Analysis (Cambridge University Press, New York, 1927), p. 115. 
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an equation first obtained by ter Haar. Equation (29) reduces to Eq. (15) when R-«. Thus the 
redundant zeros no longer appear if we cut off the potential at a large distance R and subsequently 
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make R tend to infinity after the eigenvalues have been determined from the S matrix. 
Though the potential (27) becomes the potential (10) when R becomes infinitely large, the expres- 
sion (28) goes over to (14) only when the imaginary part of k lies between +1 /2a. Above this region 


S(k) © and below this region S(k)—0. 


This is not the case, however, with the derivative of S with respect to |k| at a point where (29) 


is satisfied. At such a point 





dp d 
et 7 g(a) Spas exp(—R/2a)]+-Jig(a)J—ip-sLoe exp(—R/2a) J} 
dS __ dlk| dp 
— = —exp(—2ikR) 
d|k| J~ip(a) Jip-iLa exp(—R/2a)]+Jip(a) J—ip4ilo exp( — R/2a)] 





But for large values R 





d 
—{ J~i.(a) Jin4iLa exp(—R/2a) }} ~0, 
dp 


d 
—{Ji,(a)J--1La exp(—R/2a) }} ~~ 
dp 


so that 
dS dJ ip(x) 


J~ip(a) J ipsa exp(— R/2a) ]+Jip(a) J—ipyiLee exp(— R/2a)]~ J—ip(a) Jip-sLa exp(— R/2a) J, 


J~ip-iLa exp(—R/2a)] 


dJ (a) 
p 





J—ip-sLa exp(—R/2a) ], 











——~ —exp(—2ikR) 


d|k| d\k| J-ip(a)Jips[o exp(—R/2a) ] 


which is just the derivative of Eq. (14) with 
respect to |k|, whether the imaginary part of k 
is smaller or larger than —1/2a. From Eq. (9) 
we see that the expression c, exp(—|%,|r) repre- 
sents the correct asymptotic behavior of the 
wave function for all the discrete states as R->@. 
Now we have seen above that Eq. (28) tends to 
Eq. (14) as Ro for real values of & corre- 
sponding to positive energy states. It follows 
therefore that the S matrix for the cut-off expo- 
nential potential gives the correct asymptotic 
expression of the wave function for all the s 
states in the limit Ro. 

It is of interest to observe that for large values 
of R Eq. (28) differs very little from the expres- 
sion 


S(k) = —exp(—2ikR) 
Jip(a) J—ipsLx exp(—R/2a)] 
J ~ip(a) J ipsa exp(— R/2a)] 
when k is real; but the behavior of Eq. (30) is 





(30) 


~ae Tip) dJi(a) 
=(3) —— d|k| / a 





quite different from that of Eq. (28) when k is 
complex. This example shows that the analytic 
continuation in the eigenvalue problem has to 
be carried out with great caution. 

A somewhat different situation occurs in the 
case of the Coulomb potential cut off at a large 
distance R. As shown by Mgller, we have for 
the s states 


I'(1+i/ak) 
T'(1—i/ak) 


where a is a negative constant for an attractive 
field. The eigenvalues are determined by the 
poles where 





S(k) = (2kR)—**/*a (31) 


I'(1+i/ak)= 0, 


from which we have 
t/jak=—n (n=1, 2,3:---). 
In the vicinity of the point k=1/ian, 
S(k) ~~ (—1)***(2kR)*(T'(n+1) }*k/(k—1/ant). 
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Hence by Ea. (7) nae re : aa 
n> n n—1 ! 2 ! 3 4 

len|*=2n| hq] (2R| kel) ™[P(n+1)}*. (32) | cate tenbronaies is 

and L,' is an associated Laguerre polynomial.* is 

As R-, S(k) tends to » or 0 according as k I wish to express my sincere gratitude to tr 

is above or below the real axis, and |c,|>«, Professor Pauli for his advice on this work. pt 

The expression * Note added in proof. sa 

The ei lues given by Eq. (18 i 

|¢n| exp(— | Ral] R) = (22| Ral )#(2R| Ral)” quns Gaon: stele le the oo medal " 


theory. They appear only when we determine the ej 
<I'(n+1)-! exp(—|kal|R), (33) values from the expression of S(k) given by Eq. tif) 
ine clei the behevier of the wave function in Tl 

, F . q- . Dr, D. ter Haar has recently inf 

however, gives just the asymptotic value of ua(r) according to a recent investigation of R. or than 
at r=R, as can be seen from the well-known  ‘wfficient to discard the redundant zeros in all cases on 
‘on f a ‘ 2 the basis of the inequality (21). It should be noted that 
expression for the radial wave function of the s__ the conclusion of the present paper, namely that Heisen- 








states berg’s general condition given by Eq. (6) is not sati 
> | =~ ——— of Se) OD by Ea, (18), is indeed sh 
i re) e sign of inequali », an h ce) 
tn(r) = Nar exp(— |Ralr)Ln%(2 | k, |r), when the sign of (21) yh — m4 
va 
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The Invariant Form of Quantum Equations and the Schroedinger-Heisenberg = 
Parallelism = 
val 
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The Institute for Advanced Study, Princeton, New Jersey 
(Received October 30, 1946) 
her 
Some years have passed since the appearance of Dirac’s striking interpretation of the 
Heisenberg-Schroedinger parallelism in the equations of motion, and yet it appears that there 
still remain obscure points which have led to misunderstanding. Although the Schroedinger whi 
method suffices to deal with many practical problems, the flexibility of the Dirac-Jordan ope 
representation theory, aside from its intrinsic unity and beauty, recommends it to general eral 
theoretical investigation. The present paper considers the invariant form of the equations of 
motion in relation to the Schroedinger-Heisenberg pictures of quantum mechanics insofar as the loss 
latter refer to the “pictured” intrinsic motion of states (vectors) and observables (Hermitian con 
operators). The discussion is divided roughly into 5-sections: Operator and Vector Transforms, S 
Contact Transformations, Heisenberg-Schroedinger Operators, Equations of Motion and 
Invariant Properties, and the Density Operator and Operator Spur. Where possible, quasi- ( 
geometrical diagrams are given to illustrate the relation and transition between the several 
“pictures of the motion.” 
Notation :} ont 
a, 3, q, etc. (boldface type) = operators. 
|)=state or coordinate vector in Hilbert space. 
{| =|)*=adjoint vector. whic 
|), =transformed vector. defir 
{|a|)}=(|-a@-|)=general matrix element. able. 
|’), a’ =prototype eigenvector and eigenvalue belonging to operator a. ; 
: VW 
A. OPERATOR AND VECTOR TRANSFORMS then we say that a, is the transform of a under numt 
HEN two dynamical variables a, and @ -the similarity or collineatory transformation a: 
are related by an equation of the form t The notation here is that used in the author's forth- symb 
coming book on Perturbation Calculus and Representation Here 


«,=UaU—, (1) Theory. 
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generated by the operator U. Furthermore, if U 
is unitary, i.e, U*=U-, then the transformation 
is termed unitary or conjunctive. Similarity 
transformations in general have the important 
property that two operators so related have the 
same eigenvalue spectrum, with the transformed 
eigenvectors given by 


|)e=|a,’)=U|a’). (2) 
Thus 


a|a’)=a'|a’), 


«,U | a’) =UeU-'U | a’) =a’ a’), 


(3) 


showing that U| a’) is an eigen-|) of a, belonging 
to the eigenvalue a,’ =a’,! and conversely we can 
show that any eigenvalue of a, is also an eigen- 
value of a. The unitary transformations have the 
additional important property that real ob- 
servables (Hermitian operators) are transformed 
into real observables, and again all algebraic rela- 
tions between operators and vectors are in- 
variant. To verify these properties we note, first, 
that for a and a, real, 


a,* =U-'*aU* = a, =UeU—, (4a) 
hence 


a«U*U —U*Ue=[a, U*U ]_=0 (4b) 


which shows that U*U commutes with any linear 

operator,? and it follows that U*U is at most a 

scalar multiple of the identity operator. Without 

loss of generality we may take U*U=I, which 

condition characterizes unitary operators. 
Similarly, 


(a+6), =U(a+6)U* =a,+6,, 


(a), - U(ag)U* = UaU*UBU* = a,G,, () 


and, in particular, for 


E=E(a, 6); E,=UrU* =€(a,, 6,) 


which may be readily verified using the general 
definition of well-ordered functions of observ- 
ables. 


1With the convention that @ and a, with the same 
number of primes denote the same number, e.g-, a,’ =a’, 
but it is necessary to retain both symbols for this common 
eigenvalue to preserve the properties and meaning of 
symbols such as (a'|a,’), (a’|B|ar’), etc. 

° mposing x 4 linear operator into the sum of a 
ti-Hermitian operator. 


Hermitian and an 
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B. CONTACT TRANSFORMATIONS 


A contact transformation in quantum as well 
as in classical mechanics is defined as a trans- 
formation from one set of canonical coordinates 
and momenta, q,, p, say, to another q,,-, pyr. The 
quantal condition for a set of variables to be 
canonical is the commutator [q,, p, |_=#h, and 
being algebraic is much simpler than its clas- 
sical counterpart, though more restrictive.* The 
quantal contact transformations are in fact 
simply a special case of the unitary transforma- 
tions mentioned above, as is intrinsically evi- 
dent from the required invariance of the com- 
mutator. To show this explicitly we investigate 
the transformation from the canonical variables 
d-, Pr tO Grr, Prr, and to this end we first 
introduce the two representations in which the 
q’s and q,’s, respectively, are diagonal; the 
phases being so chosen that the correspondence 
(‘g!p|“‘g) = —ihé’(‘g—“‘g)* holds for both in their 
respective variables. 

The linear operator U is defined by its mixed 
representative 


(‘g|U|“g) =6(‘¢,—“'g),® (6) 


and is readily seen to be unitary since 


(‘g|U|“g)= f (‘q|‘“g-)d'g-("'ge | 0 | “g) 
=(‘q|“g-) (7) 


so that 


(‘q|UU*|g) = f (‘q|U|“g)d"q("g| U* |g) 


- f (q|“gr)d“ a" ge| 9) (8) 


= 5(‘gq—“g) =(‘g|I|g). 


*True canonical quantum variables must have con- 
tinuous eigenspectra, and hence the class of such variables 
in quantum theory is | limited. See P. A. M. Dirac, 
Quantum Mechanics (Oxford University Press, New York, 
1935). E. M. Corson, Phys. Rev. 70, 728 (1946). 

‘The primes slanted to the left indicate the complete 
representation, i.e., the complete set of commuting ob- 
servables, e.g., 

(q| Pe| “g) = —#h3(qx’ — gu") - - + 8(Qr-1’ —Gr_-1””) 
x ¥G@, Ge") 8(Gr4a’ — Grr’) + + *8(Qn' Gn’). 

5 The 5 function here is meaningful because of the fact 
that the eigenvalues of q, and q range over the same 
infinite spectrum, from — © to «. 
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Further 
(Ge | Gre | “g) = Gre’5 (‘Ge —“'Q), 0) 
(‘g-|Uq-|g)=5(‘¢e—“g) ar”, 
' and the equality of these expressions gives 
q-r = Uq,U-' = Uq,0*. (10) 
Similarly, 
Ge] PreU |g) = —1hd/8qre’5 (‘Ge —“Q), an 
(‘gr| Upr| “g) = 1hd/dq,"5(‘Ge—“g), 
from which it follows that 
Pr = Up,U*. (12) 


The concept of contact transformations in 
quantum mechanics can be extended to include 
dynamical systems for which no canonical coor- 
dinates and momenta exist in the strict sense, 
by defining such transformations as simply 
unitary. This generalization corresponds to the 
fact that the invariant properties associated with 
unitary transformations give them importance 
in their own right, whereas true contact trans- 
formations are an important but special case. 


C. SCHROEDINGER-HEISENBERG 
OPERATORS 


We consider the special contact transformation 
generated by the unitary operator U = exp(tHt/h) 
where ¢ is the time and H is the Hamiltonian 
operator of the dynamical system of interest, 
which at present is assumed not to depend ex- 
plicitly on ¢. 

Then if the q’s and their eigenvectors |‘g) are 
constant,® satisfying the equations 


“\'=0, “q=0 (13) 
sa. a 


their transforms q, = UqU* and |‘g,) =U|‘g) must 


_ *In the following discussion it will be seen that the 
index +r for “transformed” can in many cases be read as 
“time dependent” with equal meaning. However, zero 
subscripts corresponding to initial values are by custom 
usually omitted as the sense will always be clear from the 
equations and context. Obviously, an eigenvalue per se, 
with index merely identifies the observables in question 
in accord with the convention of footnote 1, and does not 


imply time dependence. 


satisfy the equations " 


. d 4 ) 
th—|‘q,)= —H|*q,), 
dt 





d (14) 
+ =q,H—Hq,=[q,, H]_ 
since 
‘5 = ‘a0 and ae =— ‘oy. (15) 
deh dt h 


For algebraic functions’ of non-commuting ob. 
servables we can readily verify that &(q,, D;) 
=, =UEU*; and if & does not explicitly depend 
on the time it follows that £, satisfies the equa. 
tion® 


Pas £H-H 
4 di =C€r gE. (16) 


This result, which is Heisenberg’s equation of 
motion for functions of observables, may also be 
derived (as shown by Dirac) by noting that if 
the basic vectors obey the relation 


a= H}|) 
—" a 


then an observable attached to this “moving” 
system of coordinates may be characterized by 
the equation 

,|a)= |b) (18) 


where |a) are |b) are vectors embedded in the 
moving representation. From (18) we obtain 


d 
ih t-|4)—E-Hle)= —H|b)= — HE,|a), (19) 


from which (16) follows since |a) is arbitrary. 
Both systems (13) and (14) form. proper repre- 


7 Thus far it has been possible to define only algebraic, 
well-ordered functions unambiguously in quantum theory, 
and this limitation has greatly hampered further develop 
ment of theory. See, however, P. A. M. Dirac, Rev. Mod. 
Phys. 17, 195 (1945). 7 

§ When & involves ¢ explicitly, as for example the density 
operator (Section E), the above equation must be 
generalized to P . 

tha ke =iho ke +Lks, H})- 


which again exactly parallels the classical Hamiltonian 
equations of motion. 
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sentations (actually a single infinity of bases) 
and it is convenient to refer to these coordinate 
systems and attached operators as ‘“‘fixed’’ and 
“moving,” respectively. We may note that H 
has at any time the same functional form in 
terms of the dynamical variables of either system 
(13) or (14), but in the latter, as in (16)-—(19), 
these dynamical variables are represented by 
moving linear operators. However, for the 
operator H, this distinction is in general of a 
trivial nature owing to the special form of the 
generator U. 

In the previous discussion we introduced an 
ordinary time variable in an indirect and some- 
what incidental way as a parameter in a par- 
ticular unitary transformation, and it thus 
appears that the time plays a rather subsidiary 
parametric role in quantum theory. This is in 
fact broadly speaking the case, in the sense that 
it is not possible to introduce an observable 
(operator) which would represent the time.” 

In particular the frequently met statement 
that “the time and the negative of the Hamil- 
tonian are canonically conjugate in analogy with 
classical mechanics” is quite incorrect, and par- 
ticularly surprising in view of the well-known 
quantum condition for canonical conjugacy and 
its corollary requirement of continuous eigen- 
value spectra which so sharply restrict this class 
of dynamical variables. 

The final point which we must consider to 
round out the essential structure of the theory 
is the relation between states, or equivalently, 
the vectors which represent them, at successive 
instants in time. In classical theory the equations 
of motion are relations between dynamical 
variables at different times, while in quantum 
theory they must be expected ultimately to be 
relations between the average values of the 
variables at different times, since the symbolic 
scheme is given physical content by virtue of 
the basic postulate that the number (£) =(| &|) is 
the expectation or average value of the observ- 
able & for the physical state |) provided the 
latter is suitably normalized. 


* As we shall explain later, they are representations in 
what are called the Schroedinger and Heisenberg pictures 
of quantum mechanics. 

° This remains true in the relativistic formulation even 
though the space-time variables are essentially on a 
mathematically equal footing. 
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Of course we cannot derive the equations of 
motion on purely epistemological grounds, but 
we can show that the invariance of our symbolic 
equations under unitary transformations leads to 
several equivalent and equally fundamental 
formulations. In particular, we observe that 
states and dynamical variables are not pre- 
scribed by uniquely determined vectors and 
operators since we may pass from any original 
set to the transforms |)-U]|), —E-U£EU*, etc., 
obtaining the same physical results. Hence 
what we must ultimately determine is the value 
of the expression 


d/dt®), (20) 


and we must not impute any absolute significance 
to the equations of motion insofar as this term 
usually refers to the “pictured” intrinsic motion 
of states and observables. 


D. CHANGE OF STATES IN TIME. 
EQUATIONS OF MOTION 


We extend the usual 3-dimensional sense of 
superposition of states by the assumption that 
the principle also holds in 4 dimensions, which 
from the former point of view simply means that 
states, that are linear combinations of other 
states at one instant of time, remain so at all 
times. This implies that all states of a system 
vary with time, if at all, according to a law of 
the form 


| ; +t) =G(t+ dt, t)| 5," (21) 


where G is an infinitesimal linear operator which 
is assumed to be of the form exp (—7Ht/h), so 
that ih(d/dt) | ; ?)=H|; ¢) and its conjugate 


d 
th; t| = —(; |. (22) 
dt 


Either of these equations may be called Schroe- 
dinger’s wave equation, although this term is 
more familiar in reference to their representative 
forms. If we apply the unitary transformation 
generated by U, again assuming H to be constant 
in time, we obtain an alternative (Heisenberg) 
law for the time dependence of physical states, 


“In this equation, which might be termed the 
Schroedinger postulate, the time dependence is indicated 
by the time variable after the semi-colon, to separate it 
from the quantum variables. 























namely ih(d/dt)|),=0 where 
|)}-=U| ; #)=| ; 0). (23) 
The laws (22) and (23), combined with the 





S-picture 


States as vectors: 
Dynamical variables as 
linear operators : 


To these we might properly add the corre- 
sponding equations governing representation in 
these pictures, namely (13) and (14), but it is 
usual to find only JA. and JJB. mentioned as 
equations of motion. We may quite literally say 
that quantities fixed in the S-picture are moving 
in the H-picture, and conversely. The corre- 
spondence being such that if |1), |2) denote 
states” satisfying the equations of motion (in 
either picture) and & denotes any dynamical 
variable (in the corresponding picture) not ex- 
plicitly dependent on #, then 


d 
7 || 2)=(1| {H—He| 2) (24) 
and in particular 
d 
th-X®) =((&, H]_) (25) 


which is the invariant form of the equations of 
motion, being independent of the picture used. 
To these one might add the invariant form of the 
wave equation (27) below. 


2 Here we might have written | 1; #), &, etc., the meaning 
of the indices npending on the picture, but as in most 
discussions, no distinguishing eel ts are used since 
strictly speaking — and &, represent the same observable 
and similarly for the state vectors |). The spirit of quantum 

eory permits no rigid distinction, and it is not uncommon 
to use & (without index) with both meanings in a single 
equation, the distinction being clear from the context. 
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ih el9<>=-Hi3e> 


IB. fixed: th(d/dt)E=0 


Fic. 1. Geometrical picty 
different formulations of oP 
wave equation. 





results of Section C, define equivalent quantum. 
mechanical equations which are termed the 
Schroedinger and Heisenberg pictures. 

In tabular form 


H-picture 


IA. moving: th(d/dt)|;?)=H|;¢) IIA. fixed: ih(d/dt)|),=0 


IIB. moving: th(d/dt)§,=(,, H]_ 





Equations (24)—(25) simply express the fact 
that in consequence of the linear and invariant 
properties of equations in quantum mechanics, 
the development or behavior of a system in time 
may be described as a development of the states 
themselves, the observables remaining constant 
equal to the value of the dynamical variables 
corresponding to the initial time, or conversely. 

The practical way of writing Eqs. (22) and 
(23) is in terms of representations, and as we 
have shown they must lead to identical physical 
results. If, for example, we take a representation 
in which the coordinate observables are diag- 
onal,* which means that we take the representa- 
tion (13) together with the equations of motion 
(22), or (14) with (23), then these lead to the 
familiar form of Schroedinger’s wave equation. 

S-picture :. 


d 
a ;)= f (‘g| | “g)dq(“q| ; 4) 


= H.,(‘q| ; #). (26a) 
H-picture: 


d 
oe [= f (‘gr| | “gx)d“ g(r |) + 
= op( Gr | )r- (26b) 
8 Logically this may be called Schroedinger's repre- 


sentation in either picture, but it is customary to reserve 
this term for the former fixed representation (13). 
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(a) 


Fic. 2. Geometrical picture of 

ible representations when 

the basic states are all eigen- 
states of the Hamiltonian. 





in Hla><0 
FIXED 


(b) 











in dtlce>=-Hlice> | (f) 


—» (@) 





The invariant form would then be 
. d , 
thal) = (al). (27) 


The difference is viewpoint in the above for- 
mulations of the wave equation is perhaps best 
seen in a geometrical picture," Fig. 1. 

It is thus permissible to speak of the wave 
equation in the S- and H-pictures, which are of 
course identical, as differing geometrically in the 
sense that the time dependence enters the 
generalized direction cosine (probability ampli- 
tude) from the right and left, respectively ; i.e., 
in the former it is introduced through the de- 
pendence of the states themselves, while in the 
latter the time dependence is put into the repre- 
sentation. It is important to understand that the 
q-diagonal Schroedinger representation is only a 
special case, though unquestionably the most 
important, of a whole class of representations 
which have the common property of being fixed, 
and lead to wave equations of the form (26a). 
Whether or not the right-hand side of the wave 


equation involves a sum, or integral, or both, | 


depends on the eigenvalue spectrum of the 
observables defining the representation. 

The Schroedinger equation (26a) is perfectly 
general, and we now wish to specialize to the 


case of a dynamical system for which H is 


constant, i.e., in which there are no exterior 
forces varying with time. Then there are solu- 


“The moving vectors in the diagram are indicated as 


though rotating, althoagh this is admittedly a somewhat 
loose description. 





K=H+V 


tions of the form™ 
(‘q| ‘a; t)=(‘g| ‘a)e*#"0/ (27a) 
with 


f (‘q| H|“g)d"q("q|‘a)=H"('g|‘a).. (27b) 


Here a indicates H and such other observables 
which combine with H to form a complete com- 
muting set of which | a’) is a typical eigenvector. 
According to footnote 8, a function of the 
dynamical variables not explicitly dependent on 
the time is a constant if it commutes with H. 
Such variables are the ‘“‘constants of the motion,” 
and since they must commute with H at all 
times, this usually requires that H be a constant 
(in time). The constancy of H in the Heisenberg 
picture means that 0H/d¢=0, or that H itself 
does not involve ¢ explicitly and hence is a con- 
stant in the S-picture too. Since the constants of 
the motion satisfy the condition (d/dt)a=0 in 
either picture (a, =UaU* =a), there is a certain 
inherent redundancy in the notation in this case. 
We retain the “transform” subscript, etc., in the 
present discussion, however, to.conform with the 
general scheme of presentation. 

Equation (27b) is the standard equation for 
determining the eigenvalues of H, and hence 
when H does not involve ¢ explicitly there exist 
solutions, periodic in time of the above form 
with H’ an eigenvalue of H. This special kind of 
state is such that the average value of an ob- 


18 We consider only (26a) to avoid needless repetition, 
since (26b) leads to the same conclusions with but slight 
change in notation. 
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servable is independent of time, since 


@=Ca:tlel'as =f fCastl's 


Xd'q(‘q| E|“g)d“g("g| ‘a; 2), (28) 


which according to (27a) is independent of ?. 
Thus the probability of obtaining any particular 
result when we make an observation is inde- 
pendent of the time, and such states are termed 
“stationary.” 

Let us now consider the representation we get 
when we take the stationary states to be the 
basic states. From the respective equations of 
motion we obtain the geometrical picture (Fig. 2) 
of the various possible representations when the 
basic states are all eigenstates of the Hamiltonian." 

Representation (c) with the constants of the 
motion @ diagonal is moving in the S-picture and 
is therefore fixed in the H-picture (d). Alter- 
natively, we can introduce a representation with 
the a’s diagonal which is fixed in the S-picture 
(e) and is thus moving in the H-picture (f). In 
this special case there is as indicated above a 
duplicity of motation since (c) is related to (e) 
and (d) by the same phase factor, and likewise 
(f) to (e) and (d). We would prefer to keep the 
notation consistent within each picture, but this 
should cause no difficulty. The fixed represen- 
tations in the two pictures, namely (e) and (d), 
differ only in that the basic vectors of the 
Heisenberg representation (d) satisfy the equa- 
tions of motion JJA. Nevertheless, in the equation 


Cal; ) =Car] yearn 


the left- and right-hand representatives are 
sometimes referred to as the Schroedinger and 
Heisenberg representatives of the “‘same’’ state. 
Actually the two representations do not differ in 
phase, except as a figure of speech, since the 
phase difference in the representatives arises 


16 In the H-picture 
(‘a| E,| ‘ap = SS K'a| ‘gr)d ger | Er | “ged ‘gr(“Gr | ‘ae 


which gives the same result; indeed (‘g,|&-|‘‘gr) is con- 
stant, and this condition characterizes a moving basis in 
the H-picture, just as the constancy of (‘g| &|‘‘g) character- 
izes a fixed basis in the S-picture. 

17 Properly, only the stationary state representation in 
the H-picture is called the Heisenberg representation, 
since it was the first discovered in quantum mechanics; 
Heisenberg (1925). 
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from the fact that |;#) and l)e=U| 54) (re: 
ferring now, as in the Fig. 2, to states i 
by the Hamiltonian K) represent the same 
state, but in different pictures. 

It should be noted that the representative of 
an observable & ” in the representations (c) and 
(d) varies according to the Heisenberg law , . 
constant expi(H’—H"’)t/h, the constant being 
the representative of & in (e) and (f), respec. 
tively. The above diagram also illustrates the 
exception to the fixed-moving parallelism be. 
tween the two pictures, which arises when ye 
compare systems with different Hamiltonians, 
e.g., H amd K=H+V above. The equation of 
motion of state (a) when expressed in the fixed 
Schroedinger type representation (e) leads to 
the Schroedinger wave equation. 


d 
ae | )=Lva'a| K|“a)“er| ; t) 
(29) 
= H'('a| 5 t)+Lvaa|V|“a)“al ; 2). 


Assuming a discrete spectrum for definiteness, 
The same result with appropriate change of 
notation is obtained for the transformed state 
vector |),, i.e., (b) relative to the moving Heisen- 
berg axes (f). This equation may be ascribed as 
giving the variation with time of the direction 
cosine between (a) and (e), or equivalently (b) 
and (f), due to K. Whichever picture we adopt 
to obtain (29), the matrix of V (or V,) appearing 
therein is constant except insofar as V itself 
involves ¢ explicitly. In perturbation problems, 
one is usually more interested in the separate 
effect of the perturbing potential V when con- 
sidered from the point of view of the unper- 
turbed system defined by H. Or as we may say, 
we consider only the variation in the direction 
cosine between (a) and (c), or (b) and (d), due 
to the added potential V. 

In the H-picture we have 


d 
ihLa)e=Doal'a|Ve|"a)"al)e, (30 
and in the S-picture 


d 
ida t| 1) = Lovaas t|V|"a5 "a5 #50). 
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These equations are of course identical, and 
the essential point to note is that, in both, the 
matrix of V varies rapidly according to the 
Heisenberg law ~expi(H’—H")t/h, aside from 
any explicit dependence of V on ¢. 

In general theoretical work one may use one 
or other general picture as convenience dictates, 
and it is not ordinarily necessary to use dis- 
tinguishing indices or subscripts as we have done 
here for emphasis. A note of caution may be 
added to the effect that although it is customary 
to give the laws of motion in thepreviously 
mentioned limited form (JA, JIB), and we may 
for reasons of convenience pass back and forth 
from the S- to H-pictures; they do not hold 
simultaneously. 

Thus as a final illustration 


(iteld= f('add'asCael Eel “anddax("ae1) 
(32) 


= f (; | ‘g)d‘gX‘q| E|“g)d gq] ; #) 
=(; t|E|; 2 


which follows immediately from the transforma- 
tion theory. 


E. THE DENSITY OPERATOR 


Quite frequently we have to deal with dy- 
namical systems whose specification is less com- 
plete than heretofore considered, in the sense 
that it may be unnecessary or even impossible 
to know that the dynamical system is in a speci- 
fied state at a given time. Instead, the available 
information may be only that the system is to 
be found in one or another of a number of ac- 
cessible states according to a given probability 
distribution. 

The problem thus stated is essentially a statis- 
tical one, and in classical theory leads to the 
concepts of the Gibbs ensemble, phase space, 
etc., the details of which are of no great interest 
in the present connection. Here we need only 
recall the well-known fact that the concept of 
phase space has, strictly speaking, no meaning 
in quantum mechanics due to the non-com- 
mutation of canonical coordinates and momenta 
which precludes the simultaneous assignment of 
the values required to specify a phase space. 
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Nonetheless, there exists a close quantal analogy 
to the classical ‘“‘density in phase,” and it is 
primarily the mathematics of this “quantum 
density” which we wish to consider here. We 
introduce the density operator here because it 
has in consequence of its definition certain ap- 
parently anomalous properties which can cause 
confusion, and seem best disposed of while we 
are considering the general equations of motion. 

Let us suppose that it is given that P, is the 
probability that a dynamical system under study 
is at a certain time in the mth state |; #) of a 
certain set of states accessible to it. Here we are 
implicitly assuming that the set of states cover 
a discrete range as indicated by the ordinal 
parameter m; the extension to the case of a con- 
tinuous or mixed range being trivial. As indi- 
cated, we are at present considering the S-picture 
(states as moving vectors) and in an arbitrary 
representation determined by the fixed operators 
— whose eigenspectrum may be taken as con- 
tinuous for definiteness, the wave functions will 
appear as (‘| ; 2), etc. 

For a single state the probability that the £’s 
have values in the range d‘é at £', is |(‘E| ; é)|?, 
and hence for the given distribution in the 
present case the corresponding probability will 


be 
Lal CE]; t)|2P a. , (33) 


This expression can be formally interpreted as a 
diagonal element of the quantum density 
operator @, defined directly by the tensor or 
outer product 


o=Laln; )Pa(n;t|=L,U*|n)P.(n|U, (34) 


or by its representative which is actually more 
general!® 


Célel“&)= Da Elm; Palm; t|“€). (35) 


The density operator is a rather exceptional 
if not hybrid construct since its representative 
in a fixed Schroedinger type representation as 
defined by (35) obviously varies with time. This 
is simply due to the fact that » spans a set of 
moving physical states, and hence itself explic- 
itly varies with time (Eq. (38)) in the S-picture 


18 Since in the approximate treatment of many-particle 
problems, @ can always be defined through (35) even when 
the form (34) cannot be written explicitly. 































in contrast to the more familiar observables 
which do not explicitly depend on time. The 
transform of (34) defines the density operator 0, 
in the H-picture as 


or =UpU* => ,,|")P.(n| (36) 
and 


Er] or] “&)= Den Er |m)Pa(m|“E-). (37) 


The Heisenberg density operator will also vary 
unusually with time since (‘£,| o,|‘‘€,) varies with 
time instead of being constant, as it would be 
for a simply moving observable. In fact we shall 
find that 0, itself is constant. The equation of 
motion of the’density operator is most easily 
obtained_by differentiating (34), which gives 


0 
— eH= —[e, H]_, (38) 


P,, being a constant so long as the system is 
undisturbed.” 
It immediately follows that 


af 0 (39) 
M— oe, =V, 
die 


which is also quite obvious from the constancy 
of |) in the definition (35). The formal definition 
of the density operator is especially convenient 
since from (34) it is trivially evident that o is 
Hermitian.” Further it is readily seen that has 
only positive eigenvalues, thus paralleling the 
non-negative character of the classical density in 
phase. Thus if |p’) be an eigen-|) of 9, we have 


Xa| 2; t)P n(n; t| p’)=p’ | p’) (40a) 
and 


Xia Pal (m; t| p’)|?=p'{p'|p’), i-e., p’>0.%4 (40b) 


1#The same result may be obtained somewhat less 
directly from the definition (35), and similarly (39) from 


37). 
2° Noting that |a){b|*= laa and that O=P,=1. 
_# Alternatively: @e—@?= 2n P2(1—P.)|n; t)(n;¢| from 
which 0=p’=1. 
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For a system distributed over the various states 
in this manner, the average value of an obsery. 
‘able a@ will be 


(a) = Don Pa(n; t| a|m; 2) 


-f fx. P,(m; t|E)d&E| «| “Ed EE]: 1) 
=f fceleloeeel el ae 


- f (‘| ag|*é)d*é. (41) 


The average value of an observable is given by 
summing or integrating the product ap along 
the diagonal, over the eigenvalue spectrum of the 


observables defining the representation. This . 


diagonal sum, or integral (or both) is called the 
spur of the operator product, and for any repre. 
sentation satisfies the conditions 


Sp(aBy) =Sp(Gra)=Sp(ya$). (42) 


The spur of a product of operators is invariant 
under cyclic interchange, and also under unitary 
transformation. 

Thus the physical content of the quantum 
equations of motion can be written most con- 
cisely in the form 


(a)=Sp(pa), (43) 
d 
th— a) = Sp(thoa) 
dt 
=Sp{o(aH—He)} =(eH—He). (44) 
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On the Disintegration of Negative Mesons 


M. Converst, E. PANCINI, AND O. Piccioni* 


Centro di Fisica Nucleare del C. N. R. Istituto di 
Fisica dell’'Universita di Roma, Italia 


December 21, 1946 


N a previous Letter to the Editor,’ we gavea first account 

of an investigation of the difference in behavior between 
positive and negative mesons stopped in dense materials. 
Tomonaga and Araki? showed that, becuase of the Coulomb 
field of the nucleus, the capture probability for negative 
mesons at rest would be much greater than their decay 
probability, while for positive mesons the opposite should 
be the case. If this is true, then practically all the decay 
processes which one observes should be owing to positive 
mesons. 

Several workers? have measured the ratio 7 between the 
number of the disintegration electrons and the number of 
mesons stopped in dense materials. Using aluminum, brass, 
and iron, these workers found values of 7 close to 0.5 
which, if one assumes that the primary radiation consists 
of approximately equal numbers of positive and negative 
mesons, support the above theoretical prediction. Auger, 
Maze, and Chaminade,‘ on the contrary, found 7 to be 


. Close to 1.0, using aluminum as absorber. 


Last year we succeeded in obtaining evidence of different 
behavior of positive and negative mesons stopped in 3 cm 
of iron as an absorber by using magnetized iron plates to 
concentrate mesons of the same sign while keeping away 
mesons of the opposite sign (at least for mesons of such 
energy that would be stopped in 3 cm of iron). We obtained 
results in agreement with the prediction of Tomonaga and 
Araki. After some improvements intended to increase the 
counting rate and improve our discrimination against the 
“mesons of the opposite sign,”” we continued the measure- 


TABLE I. Results of measurements on #-decay rates 
for positive and negative mesons. 











Sign Absorber III IV Hours M/100 hours 
a) + 5 cm Fe 213 106 155.00’ 67 +6.5 
b) - 5 cm Fe 172 158 206.00’ 3 

c) = none 71 69 107.45’ =I 

d) + 4cm C 170 101 179.20’ 36+4.5 
e) — 4cm C+5 cm Fe 218 146 243.00’ 2723.5 
f) — 6.2 cm Fe 128 120 240.00’ 0 
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Fic. 1. Disposition of counters, absorber; and magnetized iron plates. 
~All counters “D” are connected in parallel. 


ments using, successively, iron and carbon as absorbers. 
The recording equipment was one which two of us had 
previously used in a measurement of the meson’s mean 
life. It gave threefold (III) and fourfold (IV) delayed 
coincidences. The difference (III) —(IV) (after applying a 
slight correction for the lack of efficiency of the fourfold 
coincidences) was owing to mesons stopped in the absorber 
and ejecting a disintegration electron which produced a 
delayed coincidence. The minimum detected delay was 
about 1 usec. and the maximum about 4.5 ywsec. Our calcu- 
lations of the focusing properties of the magnetized plates 
(20 cm high; 8=15,000 gauss) and including roughly the 
effects of scattering, showed that we should expect almost 
complete cut-off for the “mesons of the opposite sign." 
This is confirmed by our results, since otherwise it would 
be very hard to explain the almost complete dependence 
on the sign of the meson observed in the case of iron. 

The results of our last measurements with two different 
absorbers are given in Table I. In this table “Sign” refers 
to the sign of the meson concentrated by the magnetic 
field. M=(III1)—(IV)—P(IV), the number of decay elec- 
trons, is corrected for the lack of efficiency (p) in our 
fourfold coincidences (~0.046). 

The value M- (5 cm Fe) is but slightly greater than the 
correction for the lack of efficiency in our counting, so 
that we can say that perhaps no negative mesons and, at 
most, only a few (~5) percent undergo 8-decay with the 
accepted half-life. 

The results with carbon as absorber turn out to be quite 
inconsistent with Tomonaga and Araki's prediction. We 
used cylindrical graphite rods having a mean effective 
thickness of 4 cm because we were unable to procure a 
graphite plate. In addition, when concentrating negative 
mesons, we placed above the graphite a 5-cm thick plate 
of iron to guard against the scattering of very low energy 
mesons which might destroy the concentrating effect of 
our magnets. We alternated the following three measure- 
ments: 


A. Negative mesons with 4 cm C and 5 cm Fe, 

B. Negative mesons with 6.2 cm Fe (6.2 cm Fe is ap- 
proximately equivalent to 4 cm C+5 cm Fe as far as 
energy loss is concerned. 

. Positive mesons with 4 cm C. 
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The comparison between A and B gave the difference in 
behavior between Fe and C, once we had established the 
fact that practically no disintegration electrons came from 
negative mesons in the 5-cm iron plate. The comparison 
between A and C gives the difference in behavior between 
negative and positive mesons in carbon. This must be con- 
sidered as a qualitative comparison because of the slightly 
different action of the magnetic field in concentrating 
mesons of different ranges (4 cm C+5 cm Fe in one case 
and 4 cm of C in the other). We could not, of course, add 
5 cm of Fe for the positive mesons too, since positive 
mesons do decay in Fe. 

The great yield of negative decay electrons from carbon 
shows a marked difference between it and iron as absorbers. 
Tomonaga and Araki’s calculation also give for carbon a 
much higher ratio of capture to decay probability for 
negative mesons, so we are forced to doubt their estimation. 
It is possible that a suitable dependence of the capture 
cross section, ¢,, on the nuclear charge, Z, might explain 
these results; however, if the ratio of the capture to decay 
probability also depends on the density as Tomonaga and 
Araki pointed out, then it would require a very irregular 
dependence on Z to also explain the cloud-chamber pictures 
of some authors* showing negative mesons stopped in the 
chamber without any decay electrons coming out. 

Concerning the difference between M, and M_ in carbon, 
we should like to point out that it is not necessary to 


assume that o, for carbon has an appreciable value for . 


negative mesons. A positive excess, (H,—H_)/(H,+H_) 
of 20 percent in the hard component, as it seems to be’ is 
sufficient to explain our results since this gives H,/H.=1.5 
which is greater than M,/M_ for carbon. Impurities in the 
graphite could also explain some preference for M,, with 
a suitable dependence of «, on Z. 

Further experiments on this subject are now in progress, 
in an attempt to calculate the capture cross section, and 
to know how it depends on Z. 


* Now Visiting Research Associate at Massachusetts Institute of 
T » Camb: . Massachusetts. 
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Further Remarks on the Redundant Zeros in 
Heisenberg’s Theory of Characteristic Matrix 


S. T. Ma 
Institute for Advanced Study, Princeton, New Jersey 
December 27, 1946 


ECENTLY Dr. Opechowski! made a remark on my 
previous communication concerning the redundant 
zeros.2 Reasoning in accordance with the current scheme 
of quantum theory, he came to the conclusion that all the 
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my note were given by the first condition, and there 

no eigenvalues given by the second condition. That is of 
course quite correct. In fact it is a result already well 
known, and that is why the term “redundant zeros” was 
introduced for the eigenvalues given by the second condi. 
tion. 

Dr. Opechowski has, -however, apparently 
the fact that I was studying the eigenvalue problem from 
the point of view of the theory of characteristic Matrix, a 
new scheme of quantum theory recently proposed by 
Heisenberg.* In applying this theory to the spherij 
symmetrical states of a particle in a central field of force, 
one determines first S(k), the eigenvalue of the charac. 
teristic matric for the spherically symmetrical states and 
real values of k. The.quantity S(k) is equal to the ratio of 
the amplitude of the outgoing wave to that of the incom. 
ing wave. One then performs an analytic continuation 
into. the complex plane of k and determines the zeros of 
S(k) in the lower half of the complex plane. Pr ing in 
this way, I obtained both sets of eigenvalues mentioned 
in my note. It is clear that my conclusion based on the 
new scheme of quantum theory cannot be invalidated by 
Dr. Opechowski’s considerations based on the current 
scheme of quantum theory. 

The question of the part played by the characteristic 
matrix as a fundamental concept in the future develop- 
ment of quantum theory need not be considered here, as 
that is irrelevant to the point at issue. 


Iw. owski, Phys. Rev. 70, 772 (1946). 
iW. ~ tag = P Payal 120; 513 ( 

- Heisen! . Zeits. f. Physi i » 673 (1943). Heisenberg’ 
third and fourth papers on the same subject are now in press. A ome 
mary of the four pa has recently been given by Heisenberg in a 
paper entitled “‘On the Theory of Elementary Particles,"’ now also in 
press. Cf. also C. Mgller, Det. Kgl. Danske Vid. Sels. Fys.-Math. Medd, 
23, No. 1 (1945); 24, No. 19 (1946); ter Haar, Physica (1946) in press; 
S. T. Ma, Phys. Rev. (1947), in press. 





Production and Isotopic Assignment of a 
90-Day Activity in Element 43* 
E. E. Motta, G. E. Boyp, AND A. R. Bros! 
Clinton Engineer Works, Monsanto Chemical Company, 
Clinton Laboratories, Oak Ridge, Tennessee 
July 31, 1946 
REVIOUS deuteron bombardments of Mo using the 
cyclotron have led to the production of a variety of 
periods associated with the chemistry to be expected from 
element 43.42 However, the mass assignment of many of 
these isotopes has not been reported. Neutron bombard- 
ments of purified samples of Ru(OH); in the Clinton self- 
sustaining chain reacting pile have been found to produce 
a number of these same activities, and, in one case at least, 
an isotopic assignment seems possible. 

Experiments have been carried out wherein growth was 
observed of a daughter activity formed by K-capture decay 
of a previously-discovered three-day Ru’ isotope.? Ex- 
tensive chemical tests in which the known six-hour 43” 
activity was employed as a monitor have shown this 
daughter activity to be an isotope of element 43, thus 
permitting its assignment to 43%. The activity has been 
observed to decay with a half-life of 935 days over three 
months time. The radiations have been found to consist 
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largely of very soft electrons of 1.5 mg/cm? half-thickness 
in Al, along with possibly a small amount of unconverted 

-radiation and probably x-rays. The half-period 
and character of the radiations possessed by this 43 isotope 
suggest that its activity is the same as that observed by 
Cacciapuoti and Segré*‘ from the deuteron bombardment 


of molybdenum. 


is based on the results of research performed under 
ae Wr 3s-058- ng-71 with the Manhattan Project at the Clinton 


i Ridge, Tennessee. 
Laboratories. Os . Phys, Rev. $2. 1252 (1937) 
, Ph " 


_N. Cacciapuoti and E. Seare 
1G. tr Seaborg and E. Segré, Phys. Rev. 55, 808 (19: 
3 W. H. Sullivan, N. R. Sleight and E. M. Glad vate communi- 


row, pri 
i Plutonium Project Report CC-1493, March 8, 1944). 
a {oCScciapuoti, Phys. Rev. $5, 110 (1939). 





Electrical Conductivity of Na-Ammonia 
Solutions at Low Temperatures 
Lu1c1 GIULOTTO AND ALBERTO GIGLI 


Istituto di Fisica “A. Volta’’ dell’Universita di Pavia, Pavia, Italia 
January 6, 1947 


A. OGG! has observed that the electrical resistance 

e of approximately one-molar sodium solutions in 
liquid ammonia decreases when rapidly frozen, so that 
the resistance of solid solutions is very low compared with 
the resistance of liquid solutions. Moreover, experiments 
carried out with an adaptation of the classical “ring 
experiment” in magnetic field, induced Ogg to think that 
these solutions may act as superconductors up to 180°K. 
With successive experiments, Boorse et al.,2 and Daunt 
and co-workers,’ could not confirm the existence of super- 
conductivity. Hodgins,‘ on the contrary, reports that four 
trials, of the numerous ones he carried out, gave positive 


. results. 


Some experiments on superconductivity and measure- 
ments of resistance carried out by us, seem to contribute 
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Fic. 1. Typical resistance curve of sodium-ammonia solution in the 
neighborhood of the solidification temperature. 
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to the clarification of this phenomenon. The solutions we 
used were approximately of one-molar concentration and 
were prepared by dipping into purified ammonia, glass 
ampoules containing sodium which had been just broken. 
The rapid freezing of the solutions was obtained by the 
falling of a small quantity of solution into a small glass 
having a thin plane bottom, previously immersed in the 
liquid air. Immediately on the bottom of the glass a solid 
disk was formed which in a few seconds reached the tem- 
perature of the liquid air. This operation was performed 
between the pole pieces of an electromagnet capable of 
providing a field of about 500 gauss. Afterwards the glass 
was rapidly extracted and brought near a magnetometer 
whose sensitivity was about 10~* gauss. The duration of 
this last operation was of the order of one second. To avoid 
too rapid heating of the solid disk, a plush disk was applied 
on the bottom of the glass which could remain saturated 
with the liquid air. Twenty trials gave negative results. 
It is to be assumed that the disk was often broken during 
the freezing, which evidently, however, could not entirely 
prevent the circulation of eventual persistent currents. 

Measurements of resistance of the solutions at different 
temperatures have been carried out, while they were 
progressively frozen or heated. The cell containing the 
solution was constructed from a small glass on whose 
bottom two strip platinum electrodes connected with a 
Wheatstone bridge, in which a resistance was continuously 
variable, were disposed, and the soldering of a very thin 
copper-constantan couple was protected from the ammonia 
by a little bed of white wax. 

More or less rapid freezing or heating of the solution 
may be obtained by simple artifices. By not too rapid 
freezing or heating (the time necessary for the solidification 
or melting of the disk was of the order of one minute), it 
was noted that the variation of the resistance as a function 
of the temperature is similar to Fig. 1, which refers to 
one particular experiment with a 0.7 molar solution. 

Evidently the resistance of the liquid solution increases 
very rapidly near the solidification temperature. Below 
this temperature a sudden decrease of resistance was 
observed, to a value much smaller than the one attained 
by the liquid solution. This phenomenon is reversible 
with the temperature. In the particular case shown in 
Fig. 1, the ratio of maximum resistance to minimum re- 
sistance is about 150, and it is to be supposed that with 
still slower freezing or heating this ratio might become 
greater. Evidently, as the solution solidifies, the resistance 
decreases very much. The increase of resistance which is 
noted immediately above this temperature is probably 
owing to the fact that, in the proximity of the solidification 
point, two liquid phases can be formed. With very rapid 
freezing, on the contrary, one may not note the increase 
of resistance which precedes the solidification, or only a 
very small increase may be noted. This probably occurs 
because solid and liquid phases are contemporaneously 
present but not in thermal equilibrium, and therefore it 
is to be assumed that the true range of this phenomenon 
does not depend on the velocity of freezing. 

We could not measure very exactly the resistance of the 
solid at the temperature of the liquid air because, in the 
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greatest number of cases, even a little below the solidification 
temperature cracks were formed in the disk, thus falsifying 
the measurements and causing an illusory increase of 
resistance. In the luckiest cases the resistance at the tem- 
perature of the liquid air was found to be almost equal to 
that of the solid just below the solidification point. 

There was still the possibility that the solid resistance 
was much lower than the one measured with this method, 
because of the resistances in the contacts. Therefore 
measurements with a potentiometer method were carried 
out. The resistance of the solid also was always found 
with this method and so in all our experiments the phe- 
nomenon of superconductivity was absent. Jn the other 
hand, the fact which we demonstrated, that the sudden 
decrease of resistance occurs just at the solidification point, 
induces us to think rather of a phenomenon depending 
just on the change of state, than of a phenomenon of super- 
conductivity. 


1R. A. Ogg, Phys. Rev. 69, 243 (1946). 
2H. A. Boorse, D. B. Cook, R. B. Pontius, and M. W. Zemansky, 


Phys. Rev. 70, 92  . 

3J. G. Daunt, M. t, K. Mendelssohn, and A. J. Birch, 
Phys. Rev. 70, 219 (1946). 
_ 4j. W. Hodgins, Phys. Rev. 70, 568 (1946). 





Relative Abundance of the Copper Isotopes and 
the Suitability of the Photometric Method 
for Detecting Small Variations 
in Isotopic Abundance 


Henry E. DucKworTH AND BENJAMIN G. HoGcG 
Scott Laboratory, Wesleyan University, Middletown, Connecticut 
January 6, 1947 


SING a Dempster double-focusing mass spectro- 
graph,' the abundance ratio of the two copper iso- 
topes was determined photometrically with a view to 
ascertaining the value of the method in detecting small 
variations in isotopic abundances. Following the general 
method of Mattauch and Ewald,** the blackening curve 
of Ilford Q11 plates was obtained by using as standards 
the isotopes of zinc, which have been accurately measured 
by Nier.* A brass rod served as one of the electrodes in a 
Dempster spark source® and provided both zinc and copper 
ions. The intensities of the mass spectrum lines were 
measured with a microphotometer, and a blackening curve 
was drawn for each photographic plate. On most plates 
there were more than one exposure and, in these cases, the 
blackening curve was drawn with considerable certainty. 
The values of the ratio Cu*/Cu® obtained by measure- 
ment of twenty-one mass spectra are listed below: 


2.258 2.335 2.242 2.284 2.135 
2.292 2.433 2.160 2.342 2.400 
2.296 2.215 2.010 2.425 2.392 
2.145 2.370 Ler 2.362 2.230 


These data lead to a mean value of the ratio Cu®/Cu® 
=2.277+0.017, the probable error being based on the 
internal consistency of the experimental results. The man- 
ner in which the present blackening curves were obtained, 
namely, in most cases by use of two or more zinc spectra 
of varying intensity, tended to cancel out small errors in 


aa. 
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the presumed abundance of the zinc standard. 
quently the isotopic constitution of copper is Cy®=¢9 48 
+0.16 percent and Cu®=30.52+0.16 percent. Flugge ang 
Mattauch* give a value for the packing fraction of 
Cu® = —8.13+0.1010~*, while Duckworth? gives Pn 
for Cu® and Cu® to be —7,92+0.25x 10-4, ion} 
value of f= —8.03+0.1010~ to both copper i 

and using the ratio given above, namely, 2.277 40.017 
one obtains a chemical atomic weight of 63.542 40.006, 
This is not in very good agreement with the accepted 
chemical atomic weight* of 63.57 but does agree splendidly 
with the recent value of 63.542 obtained by Hénigschmig 
and Johannsen.*® Ewald’s latest value* for the Cu®/Cys 
ratio is 2.330+0.032. 

Some remarks should be made regarding the efficacy of 
the method for detecting small variations from the natural 
abundance. If enough of the sample material is available 
to make, say, ten determinations of the type described 
above, the probable error in the ratio would be of the order 
of 0.020. Since there is associated with the natural ratio 
a probable error of 0.017, variations would “pr 
(used in the probable error sense) be real if they exceeded 
0.037/2.277 =0.016 of the natural ratio. As the measured 
ratios diverged still further from the natural value, the 
certainty of the variation being a real one would increase 
in the statistical manner. With the present accuracy the 
method would not be useful for variations from the natural 
abundance of less than about one-half percent. 

These measurements were made while the authors were 
at the University of Manitoba. The senior author grate. 
fully acknowledges the support of the National Research 
Council of Canada. The authors were ably assisted in the 
experiment by Miss Constance Cox, a senior student in 
physics at the University of Manitoba. 

1A. J. Dempster, Proc. Am. Phil. Soc. 75, 755 (1935). 

2 J. Mattauch and H. Ewald, Zeits. f. Physik 122, 314 (1944), 

*H. Ewald, Zeits. f. Physik 122, 487 (1944). 

+A. O. Nier, Phys. Rev. 50, 1041 (1936). 

5 A. J. Dempster, Nature 135, 542 (1935). 

6S. F lugge and J. Mattauch, Physik. Zeits. 44, 181 (1943). 

7H. E. Duckworth, Phys. Rev. 62, 19 (1942). 


8 J. Am. Chem. Soc. 63, 850 (1941). 
*O. Hénigschmid and T. Johannsen, Naturwiss. 31, 548 (1943). 





Angular Distribution of the Li’(p, a)a Reaction 


SyLvaAN Rustin, W. A. Fow.er, AND C. C. LAURITSEN 


W. K. Ratngg Dette Laboratory, California Institute of 
echnology, Pasadena, California 


December 30, 1946 


EASUREMENTS of the angular distribution of the 
8-cm alpha-particles produced in the reaction 


Li?’+H*~+2 He‘, 


have been made for the range of proton energy from 400 
to 1400 kev. Approximate agreement with previous work" 
up to 900 kev is observed, but there is only a gradual 
decrease in the asymmetry from 800 to 1400 kev. 

The apparatus used was a scattering camera,’ similar in 
principle to that described by Chadwick, ef al.,* using 
photographic detection of the emitted alpha-particles 
Figure 1 shows the apparatus schematically. The 2” x4" 
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Eastman fine-grain alpha-particle plates were located close 
to the target, the center of which was 15.1 mm from the 

of the surface of the emulsion. The paths of the 
alpha-particles which struck the emulsion along a circular 
arc at 50 mm from the target extended over a continuous 
angular range of 25° to 155° to the direction of the incident 
beam. The target surface was at 22° to the proton beam, 
which was 3/32” X 1/16” in cross section, making a target 
area of 3/32’ X5/32”. An aluminum foil separated the 
target from the photographic plate to stop protons scat- 
tered from the target. 

All the data were taken with a single thin target, con- 
sisting of a 20-kev thick coating of LiOH (measured 
parallel to proton beam), on beryllium metal. Beryllium 
was used in order to minimize x-ray emission from the 
target backing.’ Preliminary data on a clean beryllium 
target showed only a very weak production of long range 
particles, probably caused by low energy deuteron con- 
tamination of the proton beam. The yield of these was too 
low to interfere with the measurements on the 8-cm Li 
alphas, and the occasional proton tracks could usually be 
recognized by their longer range and lower grain density, 
and ignored in counting the alpha-tracks. 

The alpha-particle yield as a function of angle was 
obtained by counting individual tracks in the emulsion, 
using a microscope with a mechanical stage with microm- 
eter screw. A total of 3000-10,000 tracks were counted on 
each plate, at 12-24 different angular positions, in equal 
plate areas (usually 0.8 mm), along a circular arc at 
50 mm from the target. The counts were corrected to equal 
solid angles from the target in center-of-mass coordinates, 
and the observed angles corrected to center-of-mass co- 
ordinates. 
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Fic. 1, Semi-schematic drawing of scatt camera. Most of the 
Po + mn were made after changing to a curved foil holder, to make"the 
kness more nearly the same at all angles. 
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Fic. 2. Some of the observed angular distribution functions, in center- 
of-mass coordinates. Total number of tracks counted for each curve 
are given. Fractional P.E. of each point is roughly 0.7("/N)} where n 
is No. of points; N is No. of tracks. 


The proton charge to the target on each run was meas- 
ured, so that the yield per proton could be calculated from 
the track density and the geometry of the arrangement. 

Angular distribution functions at several different ener- 
gies are shown in Fig. 2. The distribution is of the form: 


¥(6, EZ) = Yo(E)[1+A(E) cos*#], 


where Y(6, EZ) is the yield per unit solid angle, at an angle @ 
in the center-of-mass system with respect to the direction 
of the incident protons, and at a proton energy, E, meas- 
ured in the laboratory system. Yo(Z) is the 90° yield 
function. The unit solid angle used for the yield data is a 
sphere, so that the total yield over the sphere is given by 
Y.(1+A/3). 

Figure 3 shows the function A(Z), obtained from 13 
such curves, at various proton energies. Although Swartz, 
et al.,? indicated a fairly sharp maximum in A at about 
700 kev, we have found only a gradual decrease fabove 
900 kev. The value of A reaches the same maximum 
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Fic. 3. Coefficient of cos’ term in the angular distribution as a 
function of proton Solid curve is fitted to our data. Dashed 
curves are previously pu lished results (references 1 and 2). 


observed by Swartz, et al., but does not appear to change 
rapidly between 800-1400 kev. 

The data which we obtained on yield per proton as a 
function of energy are not particularly accurate "but show 
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no indication of any resonance, at 90°, 180°, or in total’ 


yield over the sphere, in the region up to 1400 kev. The 
90° yield is a linearly increasing function up to 1400 key 
in essential agreement with previous results*’ up to 1000 
kev. 

Dr. R. F. Christy has made some calculations to relate 
these data on the angular distribution with the magnitude 
and energy dependence of the total cross section, in terms 
of the dispersion theory. Although the cross section is not 
accurately known above 0.2 Mev, it seems to show no 
pronounced indication of resonance from 0 to 1.4 Mev. A 
reasonable fit of the data is obtained with two resonances 
having !=1 Mev at 0 and 1.8 Mev roughly. The latter is 
taken to have J=2, but the former may be J=0 or J=?2, 
If o is normalized by the measurement of Haworth and 
King,* the proton width agrees well with that derived 
from the y-ray resonance at 0.44 Mev only if the 0 energy 
resonance has also J=2, 

“ & B Swart Young. eons, & oe. on and ‘SJ. Plain, ey Rox 14 oe, 
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